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ABSTRACT 


Passive  multistatic  radar  (PMR)  makes  use  of  transmission  from  opportunistic  radars  to 
detect  targets.  In  this  thesis,  a  maritime  scenario  is  created  with  multiple  merchant  vessels 
transmitting  at  S-band  (3  GHz),  serving  as  opportunistic  radars  and  a  frigate-size  warship 
acting  as  the  PMR  receiver  and  to  detect  a  low-radar  cross-section  (RCS)  target.  The 
simulations  are  carried  out  with  actual  technical  parameters  from  open  sources  to  best 
approximate  practical  performance.  To  further  improve  realism,  the  bistatic  RCS 
simulation  of  the  stealthy  Republic  of  Singapore  Navy  Formidable-class  frigate  was 
included  to  validate  the  results. 

The  simulation  results  show  that  the  multistatic  geometry  of  four  opportunistic 
transmitters  at  the  same  range  from  the  passive  receiver  with  90-degree  separation  offers 
the  best  coverage.  Passive  detection  of  a  target  of  up  to  a  radial  range  of  30  km  with 
detection  coverage  of  85%  or  better  is  possible.  This  range  coverage  is  similar  to  that  of 
the  monostatic  radar  but  lacks  in  the  area  of  detection  coverage.  The  simulations  also 
demonstrated  that  the  detection  accuracy  is  also  the  highest  using  this  same  geometry. 
The  worst-case  uncertainty  ellipse  around  the  low-RCS  target  is  less  than  150  m. 
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EXECUTIVE  SUMMARY 


Using  opportunistic  transmitters  in  the  form  of  civil  marine  radar  (S-band  and  X-band)  to 
detect  maritime  targets  passively  is  both  viable  and  comparable  in  performance  to 
monostatic  radars.  However,  in  the  case  of  detecting  low-RCS  (radar  cross  section) 
targets,  passive  multistatic  radar  (PMR)  outperforms  monostatic  radars.  The  highly 
regulated  requirements  of  the  International  Marine  Organization  (IMO)  ensure  that  most 
maritime  vessels  out-at-sea  have  the  necessary  information  to  support  both  passive 
bistatic  radar  (PER)  and  PMR.  With  the  world’s  trade  routed  via  the  sea,  there  is  no 
shortage  of  suitable  vessels  possessing  at  least  one  civil  marine  radar  (CMR)  and  one 
marine  automatic  identification  system  (MAIS)  on  the  main  shipping  routes.  Among  all 
the  suitable  vessels,  larger  ships  like  oil  tankers,  container  ships,  etc.,  have  their  CMR 
installed  at  the  highest  point,  which  extends  the  radar  horizon  and  detection  range.  The 
narrow  beamwidth  and  high  transmit  power  of  the  CMR  conveniently  lend  themselves  to 
exploitation  by  a  PER  system  in  terms  of  both  detection  accuracy  and  detection  range. 

The  unique  geometry  formed  between  the  opportunistic  transmitter,  passive 
receiver  and  unsuspecting  target  provides  exceptional  military  applications  not  offered  by 
the  current  military  technologies.  It  offers  true  passive  detection  (with  zero  transmission 
from  the  receiving  warship)  while  at  the  same  time  supports  detection  coverage 
comparable  to  a  monostatic  radar.  Passive  detection  means  there  is  no  minimum  range 
detection  limitation  that  is  unavoidable  in  a  monostatic  radar.  In  the  same  vein,  there  is 
no  need  for  ‘blanking’  of  your  electronic  support  measures  (ESM)  receiver  every  time  the 
onboard  radar  transmits.  Passive  detection  allows  for  permanent  monitoring  of  the 
electromagnetic  spectrum  for  intelligence  and  threat.  On  the  vulnerability  front,  this 
removes  the  risk  posed  by  anti-radiation  missiles  and,  at  the  same  time,  enhances  mission 
survivability,  allowing  the  warship  to  continue  tracking  targets  without  a  working  radar 
or  with  a  damaged  radar. 

This  same  unique  geometry  also  improves  the  warship’s  probability  of  detecting  a 
low-RCS  ship  that  has  been  treated  by  shaping  and  radar  absorbent  material  (RAM)  as 

compared  to  a  monostatic  radar.  The  use  of  low  frequency  S-band  CMR  poses  serious 

xvii 


limitations  to  the  use  of  RAM  treatment  on  a  low-RCS  ship.  It,  therefore,  restricts  the 
effectiveness  of  low-RCS  ships  in  the  use  of  shaping.  However,  it  is  shown  in  the  bistatic 
RCS  measurement  of  the  model  based  on  the  RSS  Formidable-class  stealthy  frigate  that 
shaping  alone  is  capable  of  reducing  the  median  RCS  to  values  not  higher  than  200  m  . 
With  all  else  being  equal,  with  the  effect  of  RCS  on  detection,  the  PER  system  has  an 
added  advantage.  The  bistatic  RCS  simulation  using  FEKO,  a  commercial  RCS 
modelling  software,  shows  more  incidence  angles  exist  that  demonstrate  a  spike  in  RCS 
returns.  Simply  put,  there  are  more  degrees-of-freedom  created  by  the  transmitter-target- 
receiver  geometry  that  has  a  higher  probability  of  directing  spikes  in  the  RCS  returns 
toward  the  receiver.  In  the  case  of  a  monostatic  radar,  there  are  limited  ways  to  maximize 
RCS  returns  because  the  shaping  of  a  low-RCS  is  designed  to  reduce  returns  in  the 
direction  of  the  incidence  angle. 

The  simplicity  of  PMR  also  lies  in  the  fact  that  there  are  no  expensive  systems  to 
build.  Instead,  all  the  information  that  the  PMR  system  needs  is  already  available  on  most 
warships.  This  information  can  be  extracted  and  used  to  carry  out  tracking. 

The  passive  detection  range  and  coverage  is  related  to  many  parameters  on  the 
transmitter,  target  and  receiver.  On  the  opportunistic  transmitter  end,  the  transmit  power 
and  antenna  gain  are  the  two  parameters  that  determine  the  maximum  detection  range. 
However,  it  is  the  relative  position  of  the  transmitter  in  the  PER  geometry  that 
determines  the  detection  coverage  or  detection  gaps.  On  the  receiver  end,  the  deciding 
parameters  are  the  sensitivity  and  bandwidth  of  the  ESM  receiver.  These  two  parameters 
affect  the  signal-to-noise  ratio  and,  thereby,  change  the  maximum  detection  range. 
Similar  to  the  opportunistic  transmitter,  it  is  the  geometry  of  the  PER  that  affects  the 
detection  gaps.  Collectively,  the  factor  that  is  used  here  to  quantify  the  geometry  is  the 
bistatic  angle  formed  between  the  transmitter  and  receiver  (see  Figure  1).  Incidentally, 
the  bistatic  angle  decides  the  bistatic  RCS,  which  is  presented  to  the  receiver  depending 
on  the  incidence  angle  of  the  transmitter’s  radar  emission. 


Figure  1.  Bistatic  Angle  -  p. 

From  the  simulation  plots  presented  in  the  thesis,  it  is  shown  that  when  more 
transmitters  are  spaced  farther  apart,  they  provide  better  overall  coverage  with  small 
detection  gaps.  The  minimum  number  of  transmitters  is  four,  placed  90  degrees  apart, 
with  one  in  each  quadrant.  To  achieve  maximum  coverage  with  this  particular 
combination  of  system  parameters  and  transmitter  constellation,  the  baseline  distance  of 
18,000  m  is  recommended.  By  rotating  the  quadrant  placement  of  transmitters,  one  could 
direct  the  detection  gaps  away  from  a  specific  bearing  of  choice  to  improve  detection. 

Using  more  than  four  transmitters  in  the  suggested  constellation  does  not 
significantly  extend  detection  range.  However,  it  is  shown  in  the  eight-transmitter 
simulation  that  it  does  improve  detection  coverage  from  80%  to  95%  through  the 
reduction  of  detection  gaps.  This  translates  to  the  ability  of  the  PMR  to  track  a  target 
continuously,  provided  the  target  offers  a  strong  bistatic  return  to  all  eight  transmitters.  In 
all  likelihood,  a  low-RCS  target  would  have  presented  at  least  two  to  four  strong  RCS 
returns  to  either  the  quadrant-based  constellation  using  four  transmitters  or  the  two¬ 
layered  quadrant-based  constellation  using  eight  transmitters.  If  the  availability  of  such  a 
constellation  is  not  possible  out-at-sea,  one  should  position  oneself  to  create  such  a 
constellation  without  jeopardizing  a  passive  posture. 
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There  are  two  dominant  factors  that  affect  detection  accuracy,  namely,  the 
beamwidth  of  the  transmitter  and  the  direction-finding  (DF)  accuracy  of  the  receiver.  It  is 
shown  that  the  worst-case  measurement  error  is  approximately  ±  120  m  for  a  low-RCS 
114  m  frigate  at  a  range  of  20  km.  In  most  cases,  a  lower  measurement  can  be  achieved 
using  multiple  transmitters  by  overlapping  their  individual  uncertainty  ellipses.  It  is  found 
that  the  four-transmitter  and  eight-transmitter  constellations  recommended  for  maximum 
detection  coverage  are  also  optimized  to  reduce  the  target’s  uncertainty  ellipse.  It  is 
shown  that  with  an  eight-transmitter  constellation,  the  overlapped  area  by  all  the 
uncertainty  ellipses  has  given  rise  to  a  measurement  error  of  20  m  for  the  same  low-RCS 
target  at  a  range  of  6400  m. 

In  summary,  the  use  of  opportunistic  transmitters  in  the  form  of  civil  marine  radar 
(S-band  and  X-band)  to  detect  maritime  targets  passively  is  both  viable  and  comparable 
in  performance  to  monostatic  radars.  The  PMR  with  four  or  more  transmitters  offers  a 
detection  range  of  up  to  30  km  with  detection  coverage  of  85%  or  better.  The  associated 
worst-case  detection  error  is  approximated  to  be  ±  120  m  for  a  low-RCS  114  m  frigate  at 
a  range  of  20  km.  On  top  of  that,  the  PMR  offers  added  capability  to  defeat  low-RCS 
targets  without  the  need  to  build  extra  combat  systems.  All  information  required  for  the 
PMR  to  work  is  available  in  most  modern  warships. 
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I.  INTRODUCTION 


A.  OVERVIEW 

The  ability  to  detect  targets  passively  has  been  greatly  sought  after  since  the 
development  of  modem  radar.  The  term  ‘passive’  detection  of  targets  should  be 
interpreted  as  the  ability  to  detect  targets  without  deliberate  or  active  transmission  by  the 
sensor  on  the  target  of  interest.  In  this  deadly  game  of  cat-and-mouse,  the  one  that  gets  to 
‘see’  first  strikes  first.  In  modem  day  warships,  new  technologies  and  techniques  like  low 
probability-of-intercept  (LPI)  radars,  low-radar  cross-section  (RCS)  design,  radar 
absorbent  material  (RAM)  coating  and  ultra-sensitive  electronic  warfare  (EW)  receivers 
have  changed  the  nature  of  the  game.  On  the  one  hand,  the  risk  of  being  detected  by 
ultra- sensitive  EW  receivers  has  forced  military  ships  to  go  on  a  strict  radar  transmission 
policy  and  leverage  other  means  of  sensing,  such  as  ship-borne  unmanned  aerial  vehicles 
(UAV)  or  passive  infrared  sensor  and  track  (IRST),  for  target  detection.  On  the  other 
hand,  ships  equipped  with  capable  EPI  radars  still  have  to  contend  with  weak  returns 
from  low-RCS  designed  ships.  Eor  both  cases,  the  effective  detection  range  in  most 
maritime  missions  is  generally  reduced  to  slightly  beyond  the  visual  range  with  the 
exception  of  over-the-horizon  (OTH)  radar. 

A  proposed  solution  to  resolve  this  deadlock,  which  is  investigated  in  this  thesis, 
is  the  use  of  multiple  pairs  of  passive  bistatic  radar  (PER)  to  detect  targets  using  common 
radar  transmission  by  opportunistic  transmitters.  This  approach  offers  many  advantages. 
Eirstly,  there  is  no  radar  transmission  by  a  source  eliminating  the  risk  of  being  intercepted 
by  the  target’s  EW  receiver,  or  worse,  by  anti-radiation  missiles  (ARMs).  Detection 
occurs  passively,  thus  advancing  one’s  observe-orient-decide-act  (OODA)  loop  while 
denying  the  target’s  own  OODA  loop.  Secondly,  the  unique  geometry  formed  by 
separating  transmitter  and  receiver  reduces  the  effectiveness  of  low-RCS  design,  which  to 
date  has  only  been  optimized  to  reduce  monostatic  returns.  Simulations  have  shown  that 
the  RCS  of  stealth  ships  changes  with  different  bistatic  angles  p,  and  some  of  these  angles 
have  considerably  larger  RCS  returns.  This  unique  effect  may  support  detection  of  targets 
that  are  not  in  the  line-of-sight  (EOS)  of  one’s  radar.  In  the  maritime  environment,  the 

I 


opportunistic  transmitters  are  civil  marine  radars  (CMR)  operating  in  S-band  (2.9  GHz  - 
3.1  GHz)  and  X-band  (9.2  GHz  -  9.5  GHz)  installed  on  all  seagoing  vessels  regulated  by 
the  International  Maritime  Organization  (IMO)  under  Regulation  19  in  safety  of  life  at 
sea  (SOLAS),  Chapter  V.  The  use  of  S-band  transmission  as  a  source  also  improves  the 
chance  of  detection  as  most  RAM  treatment  is  frequency  specific  and  wavelength 
dependent.  As  the  wavelength  increases,  the  thickness  of  the  corresponding  RAM 
material  increases;  however,  the  RAM  treatment  is  not  effective  at  lower  frequencies.  In 
addition,  the  abundance  of  opportunistic  transmitters  out-at-sea  allows  one  to  choose 
from  different  combinations  of  bistatic  pairs  to  optimize  detection  coverage,  probability 
of  detection  and  range  accuracy.  Next,  all  hardware  required  for  the  PER  to  work  is 
currently  available  on  most  warships.  There  is  no  requirement  for  additional  equipment. 
Lastly,  it  enhances  the  mission  capability  of  warships  by  allowing  target  detection  with  a 
defective  or  damaged  radar. 

In  Figure  I,  a  typical  bistatic  radar  geometry  is  presented  with  the  transmitter  (Tx) 
and  receiver  (Rx)  separated  by  the  baseline  distance,  also  known  as  the  direct  path  or 
line-of-sight  (LOS)  of  transmission,  denoted  by  L.  The  distance  from  the  transmitter  to 
the  target  is  denoted  as  Rj ,  while  the  distance  from  the  target  (Tgt)  to  the  receiver  is 
denoted  as  .  The  angle  formed  between  Rj  and  is  known  as  the  bistatic  angle  p. 
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Figure  1.  Bistatic  geometry  showing  the  Tx,  Tgt  and  Rx.  From  [1]. 
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1.  History 

The  concept  of  bistatic  radar  is  not  new.  In  fact,  the  oldest  form  of  radar  is  of  the 
bistatic  kind,  which  operated  in  the  forward  scattering  configuration  for  aircraft  detection 
in  the  1930s  [1].  The  operation  of  this  forward  scatter  radar  {fi  =  180°)  is  based  on 
Babinet’s  principle  that  states  that  a  target’s  RCS  is  greatly  enhanced  when  observed  by 
the  receiver  in  a  bistatic  configuration  with  the  transmitter,  target  and  receiver  on  the 
same  baseline,  arranged  in  the  order  stated  earlier  [2].  Although  the  forward  scatter  radar 
supports  extended  detection,  it  is  restricted  in  practical  applications,  serving  more  like  a 
‘trip-wire’.  Initially,  little  was  known  about  the  benefits  that  bistatic  radar  brings. 
Concurrently,  the  monostatic  radar  technology  was  gaining  momentum  with  the 
invention  of  the  magnetron  during  World  War  II.  This  technology  gradually  matured  and 
ushered  in  the  birth  of  pulsed  radar,  moving  target  indicator  (MTI)  radar,  phased  array 
radar  and  synthetic  aperture  radar  (SAR)  [3]. 

Thirty  years  later,  in  the  1960s,  a  bistatic  radar  exploiting  a  transmitter  of 
opportunity,  fitting  the  definition  of  a  PBR,  was  mentioned  in  a  paper  by  Rittenbach  and 
Fishbein  [4]  and  demonstrated  a  more  advanced  application  of  the  bistatic  radar.  The 
most  famous  bistatic  system  built  in  the  1960s  was  the  Sugar  Tree  system,  a  high 
frequency  (HF)  OTH  bistatic  receiver  that  exploited  the  Soviets’  international  broadcast 
transmitter  to  detect  missile  launches  [4].  By  then,  the  mainstream  use  of  radar  was 
dominated  by  monostatic  radar,  closely  followed  by  counter-detection  techniques  known 
as  EW.  Concepts  like  LPI  radars,  low-RCS  designs,  and  RAM  were  explored  and 
pursued,  leaving  bistatic  radar  concepts  relatively  untouched. 

In  the  1990s,  a  renewed  interest  in  bistatic  radar  occurred  as  a  possible  solution  to 
defeat  anti-radiation  missiles  and  low-RCS  designed  platforms.  With  better  technologies 
and  computing  capabilities,  many  possible  applications  of  bistatic  radars  were  re¬ 
examined.  The  most  notable  systems  developed  during  the  period  are  summarized  in 
Table  1. 
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2. 


Bistatic  Radar  Research  Areas 


The  current  research  interest  in  bistatic  radar  and  multistatic  radars  can  be  broadly 
classified  into  three  areas:  (1)  target  detection  and  parameter  extraction,  (2)  multistatic 
radar  systems  and  configuration,  and  (3)  target  tracking  and  accuracy. 

Under  category  (1),  the  forward  scattering  configuration  of  bistatic  radar  is 
combined  with  inverse  SAR  (ISAR)  techniques.  The  bistatic  radar,  particularly  in  the 
forward  scatter  region,  is  used  to  detect  crossing  maritime  targets  cutting  through  the 
baseline  formed  by  the  transmitter  and  receiver.  Target  detection  is  accomplished  through 
differentiating  the  target’s  Doppler  frequency  from  the  sea  clutter.  Using  the  rate  of 
change  of  the  Doppler  frequency,  we  can  approximate  the  target’s  speed  and  length  with 
an  error  of  20%  and  an  associated  probability  of  detection  of  80%  based  on  the  70 
maritime  targets  tracked  [5].  This  system  allows  for  autonomous  extraction  of  target 
parameters  in  support  of  target  identification  and  classification.  A  separate  paper  went 
further  to  investigate  the  differences  between  ISAR  returns  from  forward  scatter  and 
backscatter  [6].  In  that  experiment,  an  additional  transmitter,  operating  at  a  slightly 
different  frequency,  is  collocated  near  the  receiver  to  allow  simultaneous  capturing  of 
both  the  forward  scatter  and  backscatter  when  the  ship  crosses  the  baseline.  The  results 
confirmed  that  the  Doppler  frequency  returns  collected  from  the  forward  scatter  and 
backscatter  have  distinctly  different  profiles  in  the  Doppler  versus  time  plot. 

The  forward  scatter  returns  have  a  better  reproduction  of  the  ship’s  silhouette, 
most  likely  due  to  the  diffraction  of  electromagnetic  (EM)  waves  around  the  ship’s 
superstructure.  As  for  the  case  of  the  backscatter  returns,  the  Doppler  frequency  is  almost 
uniformly  spread  out  across  the  ship  structure.  In  summary,  there  is  great  potential  in 
combining  these  two  returns  to  enhance  the  capability  to  detect,  classify  and  identify 
targets. 
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Table  1 .  Parameters  of  significant  passive  bistatic  radar  programs 
designed  and  tested  for  air  surveillance.  From  [4]. 


System 

Silent 

Sentry™ 

TV-Based 
Bistatic 
Radar  (I) 

TV-Based 
Bistatic 
Radar  (I) 

FM  Radio- 
based 

Bistatic 

Radar 

Multistatic 
HDTV- 
based  Radar 

Developer 

IBM,  now 
Lockheed 
Martin 

University 

College 

London 

DERA 

United 

Kingdom 

NATO 

SAIC,  U.S. 
Army 

Decade 

1980  -  2000 

1980 

1980 

2000 

2000 

Configuration 

Multistatic 

Rx  -1,  Tx  - 
Up  to  6 

Bistatic 

Bistatic, 
Near  forward 

scatter 

Bistatic 

Multistatic 

Rx  -4,  Tx  -  1 

Transmitter 

FM  -  realtime 

TV  -  Non- 

TV  -  Near 

FM- 

HDTV  - 

Operation 

TV  -  Non¬ 
realtime 

realtime 
(543  MHz) 

realtime 
(543  MHz) 

realtime 

realtime 

Baseline 

100  km 
typical 

12  km 

150  km 

50  km 

10  km  typical 

Target 

Aircraft, 

Missile 

Launch 

Aircraft 

Aircraft 

Aircraft 

Aircraft 
below  5000ft 

Target  Data 

Range, 

Range, 

Doppler, 

Range, 

Range, 

Doppler 

bearing 

Bearing 

Bearing 

Doppler, 

Bearing 

Doppler 

Measured 

Rm  =  100km- 

Rm  ~  25km 

Rm  ~  160 

Rm  ~  175 

Rm  ~  30  km 

150km 

Occasional 

km 

km. 

Target 

Performance 

2D  track  on 

A/C 

Detections 

Achieved 

location  via 

Aircraft 

detections, 

on  high  and 

with 

range  multi- 

3D  track  on 

but  mostly 

medium 

innovative 

lateration 

missile  launch 

negative 

altitude  A/C 
but  only  1/3 
tracked 

direct  path 
excision 

Ghost 

excision  via 
doppler 

A  PER  system  is  an  attractive  and  low-cost  choice  when  it  comes  to  detection  of 

low-flying  aircraft  traversing  a  terrestrial  region.  Many  terrestrial  transmitters  such  as 

frequency  modulated  (FM)  radio  antennas,  cellphone  base-stations  and  digital  radio 

antennas  provide  the  transmission  source  for  PER  to  detect  these  aircraft.  The  paper  by 

Griffiths  and  Eaker  [7]  evaluates  each  of  the  opportunistic  transmitters  and  their  detection 

ranges.  The  result  shows  the  FM  radio  station  at  Wrotham  in  England  has  a  detection 

range  of  up  to  27  km  with  a  minimum  signal-to-noise  ratio  (SNR)  of  15  dE  using  a 

coherent  integration  period  of  0.1  s.  In  the  case  of  a  low -power  digital  radio  transmitter, 

the  detection  range  is  reduced  to  9  km  with  the  same  parameters.  The  cellphone  base 
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station  reduced  the  detection  range  to  1.2  km,  the  shortest  of  the  three.  However,  the 
extensive  network  of  cellphone  base-stations  makes  up  for  the  short  detection  range  with 
good  detection  coverage.  Delving  into  how  the  transmission  waveforms  affect  detection, 
the  paper  found  that  wideband  signals  like  FM  and  digital  radio  introduced  added 
complexity  in  the  bistatic  and  multistatic  ambiguity  function  and  concluded  that  the 
narrowband  signal  is  more  suitable  for  PER  system.  In  effect,  the  varying  bandwidth  and 
dynamic  amplitude  range  of  the  broadcast  music  makes  detection  and  tracking  accuracy 
extremely  difficult  to  quantify. 

Data  fusion  is  the  combination  of  information  on  the  same  target  from  different 
sources  (multiple  radars)  to  improve  probability  of  detection  and  tracking  accuracy.  Even 
in  network-centric  warfare  concepts,  it  is  applied  almost  exclusively  to  monostatic  radar 
sources  where  system  parameters  such  as  range  bin,  tracking  errors,  and  SNR  are  known 
and  linear.  In  bistatic  and  multistatic  radar  systems,  these  parameters  change  in  a  non¬ 
linear  fashion  with  respect  to  the  geometry  of  the  transmitter,  receiver  and  target.  In  an 
attempt  to  better  understand  whether  data  fusion  yields  a  similar  advantage  in  a 
multistatic  netted  radar  system,  two  sets  of  commercial-off-the-shelf  (COTS)  transmitter 
and  receiver  pairs  operating  at  2.4  GHz  were  set  up  to  form  a  bistatic  angle  of  20°  with 
the  target  [8].  The  results  showed  that  SNR  improved  by  3  dB  with  data  fusion  of  both 
monostatic  and  bistatic  returns  as  compared  to  only  monostatic  returns.  The  experiment 
also  showed  that  the  fusion  of  non-coherent  returns  does  not  improve  the  SNR. 

The  greatest  advantage  offered  by  the  PER  configuration  lies  in  its  ability  to 
leverage  the  EM  emission  of  opportunistic  transmitters  to  detect  targets.  This  allows  the 
receiver  to  detect  and  track  targets  without  the  target  being  aware  of  the  radar’s  presence. 
One  of  the  key  challenges  of  this  approach  is  to  ensure  consistent  detection  coverage  and 
good  tracking  accuracy.  One  of  the  solutions  proposed  is  to  derive  a  method  to  identify 
the  suitable  transmitter-receiver  geometry  pair  to  achieve  good  coverage.  By  applying  the 
Cramer-Rao  lower-bound  (CREB)  to  the  bistatic  ambiguity  function,  we  find  that  the 
results  show  that  regardless  of  the  transmitter-receiver  pair  chosen,  the  range  and  velocity 
error  is  lowest  when  the  target  is  close  to  and  within  the  vicinity  of  the  transmitter  and 
receiver  [9].  The  error  increases  in  a  radially  outwards  manner  as  the  target  moves 
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farther.  In  summary,  to  achieve  good  detection  coverage  and  accuracy,  there  should  be  a 
transmitter  or  receiver  in  each  sector-of-interest  centered  at  the  receiver.  In  one  particular 
configuration,  an  array  of  transmitters  can  be  chosen  to  form  range  rings  around  the 
receiver. 

The  ambiguity  function  of  the  radar  waveform  determines  the  accuracy  of  target 
tracking  performance  in  both  range  and  velocity.  In  the  case  of  the  bistatic  radar 
(multistatic  included),  the  ambiguity  function  incorporates  many  non-linear  variables 
created  by  the  spatial  separation  of  the  transmitter  and  receiver.  Based  on  simulation,  it 
was  found  that  range  resolution  worsens  with  the  increase  of  the  bistatic  angle.  In  the 
Doppler  frequency  domain,  larger  detection  gaps  and  tracking  errors  are  also  created  due 
to  the  geometry  of  transmitter-target-receiver  [10].  In  a  later  paper  [11],  the  authors 
extend  the  investigation  of  the  ambiguity  function  to  cover  coherent  and  non-coherent 
multi  static  radars.  The  results  agreed  with  the  earlier  findings  and  went  further  to 
emphasize  the  importance  of  time  synchronization  and  positional  accuracy  of  the 
transmitter  and  receiver  to  achieve  good  range  accuracy.  It  also  states  that  coherence 
between  the  transmitter  and  receiver  improves  the  range  resolution  and  accuracy. 

From  [12],  improvement  in  detection  and  tracking  accuracy  can  be  achieved  by 
fusion  of  data  from  multiple  domains,  e.g.,  time-of-arrival  (TOA),  phase-on-arrival 
(POA)  and  angle-of-arrival  (AOA).  This  approach  reduces  the  tracking  error  by 
overlapping  the  area  of  uncertainty  of  the  target  and  taking  the  overlap  area  (smallest)  of 
the  three  as  the  final  range  resolution.  It  also  adds  frequency  diversity  by  transmitting  a 
range  of  frequencies  with  multiple  transmitters  to  further  enhance  performance  against 
targets  of  different  dimensions  that  are  moving  at  different  speeds.  The  results  also  show 
that  further  improvement  in  tracking  accuracy  can  be  achieved  with  a  larger  separation 
between  transmitters. 

In  summary,  there  is  a  strong  interest  in  a  PBR  system  working  with  opportunistic 

transmitters.  Experimental  systems  of  this  type  are  almost  always  used  exclusively  to 

detect  aircraft  in  a  terrestrial  environment.  The  choice  of  transmitter  waveform  is  almost 

always  broadband,  continuous  transmission  with  or  without  FM  chirp.  In  the  maritime 

domain,  the  use  of  bistatic  systems  is  almost  always  in  the  forward  scatter  mode  with  a 
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configurable  transmitter.  The  focus  lies  in  target  parameter  extraction  and  classification 
using  ISAR  techniques.  Multistatic  systems  with  good  detection  range  and  accuracy 
almost  always  maintain  coherency  between  transmitter  and  receiver.  It  was  also  found 
that  a  narrowband  transmitter  waveform  performs  better  in  PER.  Good  spatial  placement 
of  the  transmitters  by  sector  with  sufficient  distance  improves  detection  coverage.  In  the 
case  of  tracking  accuracy,  the  importance  of  having  good  transmitter  and  receiver 
position  information  and  good  time  synchronization  cannot  be  understated.  One  could 
improve  range  accuracy  by  choosing  transmitters  near  to  intended  targets  and  keeping  the 
bistatic  angle  between  transmitter-receiver  to  between  30°  and  160°. 

B.  THESIS  OBJECTIVES 

The  goals  of  this  thesis  are  to  1)  evaluate  the  suitability  of  using  passive 
multistatic  radar  to  detect  low-RCS  marine  targets  using  civil  marine  radar,  2)  propose 
the  optimal  multistatic  radar  system  configuration  for  best  detection  coverage,  and  3) 
quantify  the  detection  accuracy  of  the  optimal  configuration  using  MATLAB  simulation. 
The  simulation  results  are  obtained  using  actual  system  parameters  extracted  from  civil 
marine  radar  equipment  manuals  or  open  sources  materials. 

The  proposed  multistatic  radar  system  configuration  is  a  single  receiver  (warship) 
and  multiple  transmitters  (merchant  ships)  configurations.  The  scenario  assumes  a  single 
warship  is  tasked  to  carry  out  comprehensive  surveillance  of  a  sector  out-at-sea  to 
enhance  maritime  domain  awareness.  Presumably,  the  warship  is  positioned  somewhere 
in  the  South  China  Sea  to  monitor  ships  going  into  and  out  of  the  Singapore  Straits.  The 
targets  are  assumed  to  have  low-RCS  design,  operating  with  either  LPI  radars  or  search 
radars  with  infrequent  and  irregular  scans.  The  own-force  warship  is  equipped  with  a 
broadband  EW  receiver  with  direction-finding  (DF)  capability.  The  merchant  ships  and 
the  warship  are  equipped  with  class-A  maritime  automatic  identification  system  (MAIS) 
regulated  under  SOLAS,  which  functions  like  identification  friend  or  foe  (IFF)  on 
aircraft,  sending  out  information  such  as  identity,  global  positioning  system  (GPS) 
position,  course  and  speed  at  short  intervals  on  a  very  high  (VHF)  frequency.  The 
pictorial  representation  of  the  scenario  is  illustrated  in  Figure  2. 
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Figure  2.  Pictorial  representation  of  the  PER  setup. 


Typically,  the  process  of  target  detection  begins  with  the  warship  receiving  all 
merchant  ships’  MAIS  transmissions  on  its  own  transponder.  The  warship  forms  a 
tactical  picture  of  all  known  targets  on  its  combat  management  system  (CMS)  console 
screen.  Using  the  onboard  EW  and  DF  receiver,  each  of  the  CMR’s  transmissions  in  both 
S-band  and  X-band  can  be  correlated  with  the  individual  merchant  ships.  With  these  two 
sets  of  information,  the  EW  DF  receiver  is  now  able  to  discern  the  same  civic  marine 
radar’s  transmission  parameter  coming  from  two  different  bearings  due  to  the  direct  path 
and  reflection  off  an  unknown  target.  The  difference  in  the  time-of-arrival  of  the  direct 
and  indirect  signals  is  used  to  compute  the  range  of  the  unknown  target  from  the  warship. 
This  set  of  range  and  bearing  information  has  inherent  errors  due  to  DF  error,  data  age 
and  uncertainty  in  the  positional  information  of  the  opportunistic  transmitter.  The  error 
can  be  minimized  by  averaging  and  combining  multiple  range  computations  of  the  same 
target  using  another  bistatic  pair. 
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The  first  step  is  to  determine  the  maximum  possible  detection  range  based  the 
radar  horizon  and  the  maximum  detection  range  using  the  bistatic  radar  equation.  The 
latter  takes  into  consideration  system  and  physical  parameters  like  antenna  gains, 
transmission  losses,  frequency,  bistatic  RCS,  antenna  beamwidth,  bistatic  geometry  of 
transmitter-target-receiver,  receiver  sensitivity  for  detection,  etc. 

Next,  evaluation  of  the  combined  detection  range  for  different  bistatic  geometries, 
based  on  the  SNR  at  the  receiver,  is  presented.  Based  on  the  coverage  characteristics  of 
individual  bistatic  geometry,  combinations  of  multiple  bistatic  pairs  at  different  angles 
and  ranges  are  simulated  to  determine  the  optimal  detection  coverage.  The  simulations 
are  run  by  varying  frequency  between  S-band  and  X-band  while  combining  them  with 
different  bistatic  pairs. 

Finally,  the  range  accuracy  is  determined  for  each  of  the  previously  mentioned 
simulation  runs.  The  objective  is  to  evaluate  the  worst-case  position  uncertainty  of  the 
target  and  assess  how  choices  of  different  bistatic  pairs  are  affecting  the  results.  The 
combination  of  these  two  results  provides  the  basis  by  which  the  receiver  can  choose 
different  combinations  of  bistatic  pairs  for  optimal  detection  coverage  and  detection 
accuracy. 

C.  THESIS  OUTLINE 

The  remaining  parts  of  the  thesis  are  organized  as  follows.  The  background 
theoretical  concepts  needed  to  understand  and  carry  out  the  analysis  are  covered  in 
Chapter  II.  Starting  with  Chapter  III,  we  present  simulation  results  in  a  progressive  order 
with  the  detection  coverage  of  a  single  bistatic  pair.  Different  combinations  of  multistatic 
geometry  are  presented  in  Chapter  IV  with  analysis  on  the  optimal  detection  coverage. 
The  expected  detection  accuracy  of  different  multistatic  geometries  and  possible 
improvements  are  presented  in  Chapter  V.  Finally,  the  conclusions  drawn  based  on  the 
results  obtained  and  suggestions  for  future  work  are  found  in  Chapter  VI. 
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II.  BACKGROUND 


A.  CIVIC  MARINE  RADAR  AND  MAIS 

Probability  of  target  detection  by  radar  is  determined  by  the  SNR.  A  higher  SNR 
increases  the  probability  of  detection.  This  is  also  true  in  the  case  of  a  PER;  the  only 
difference  is  that  the  CMR’s  signal  comes  in  many  forms  and  may  not  be  designed  to 
optimize  detection.  Therefore,  the  first  step  is  to  determine  the  availability  of  the  said 
opportunistic  transmitters,  their  performance  specifications  and  operational 
characteristics. 

The  purpose  of  CMR,  as  defined  by  the  IMO,  is  to  “assist  in  safe  navigation  and 
in  avoiding  collision  by  providing  an  indication,  in  relation  to  one’s  ship,  of  the  position 
of  other  surface  craft.”  The  technology  of  the  CMR  is  heavily  influenced  by  the 
regulations  issued  by  the  UN-based  IMO.  CMR  is  restricted  to  operate  in  the  3  GHz 
range  (for  good  range  performance)  and  the  9  GHz  range  (for  enhanced  detection 
performance  in  rain  and  fog).  Since  2004,  IMO  has  encouraged  the  use  of  coherent  radars 
in  place  of  the  simple  pulsed-radar  CMR  for  improved  detection  performance  against 
heavy  sea  clutter  conditions  [13].  The  New  Technology  Radars,  as  they  are  known,  do 
not  have  restrictions  on  the  type  of  waveforms  that  can  be  used.  Some  companies  have 
developed  LPI  CMR  that  is  marketed  for  its  low  power  requirement  and  improved 
reliability  [14]. 

According  to  IMO  regulations,  any  vessel  exceeding  300  gross  metric  tonnage  is 
to  be  fitted  with  a  CMR.  This  translates  to  more  than  50,000  ships  sailing  in  international 
waters,  and  many  of  them  are  fitted  with  two  radars  operating  in  different  frequency 
bands  to  enhance  navigation  safety.  This  same  regulation  also  applies  in  the  carriage  of 
MAIS  but  with  the  coverage  expanded  to  include  ships  exceeding  500  gross  metric 
tonnage  not  on  international  voyage  and  all  passenger  ships  irrespective  of  size  [15].  A 
good  estimate  of  the  number  of  ships  installed  with  CMR  and  MAIS,  passing  through  the 
Singapore  and  Malacca  Straits  is  shown  in  the  Figures  3  and  4.  Notice  that  in  both  scenes. 
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there  are  over  500  vessels  that  are  within  the  reception  range  of  the  shore -based  MAIS 
system  as  captured  by  a  free  online  resource,  www.marinetraffic.com. 


Figure  3.  Snapshot  of  MAIS  targets  along  Malacca  Straits  dated  Feb.  20,  2014. 


Figure  4.  Snapshot  of  MAIS  targets  along  Singapore  Straits  dated  Feb.  20,  2014. 


According  to  a  report  by  the  U.S.  Energy  Information  Administration,  over 
60,000  seagoing  vessels  transit  through  the  Malacca  Straits  (passing  through  the 
Singapore  Straits)  every  year  [16].  Assuming  even  distribution,  one  expects  a  daily 
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average  of  approximately  150  vessels  transiting  though  the  Singapore  Straits  on  the 
Indian  Ocean-South  China  Sea  shipping  route.  To  further  narrow  this  down,  the 
assumption  made  in  the  simulation  is  the  use  of  frigate-class  warship  as  own  ship  and  the 
largest  merchant  vessels  that  conform  to  the  MALACCAMAX  to  determine  the  radar 
horizon.  MALACCAMAX  is  the  term  used  in  Naval  Architecture  to  denote  the  largest 
seagoing  vessel  with  a  maximum  draft  of  25  m  that  can  fit  through  the  Malacca  Straits. 
According  to  [17],  the  numbers  of  conventional  cargo  ship,  tanker,  container  ships  and 
roll-on  roll-off  (RORO)  ships  passing  through  the  Malacca  Straits  in  1994  was  11620, 
9688,  5244  and  1130,  respectively.  That  translates  to  an  average  of  76  vessels  on  any 
given  day. 

To  compute  the  radar  horizon,  two  parameters  are  needed  according  to  the  radar 
horizon  formula  (standard  atmosphere)  [14]: 

R.^=2.21^,+2.2\^  (1) 

where is  the  range  in  nautical  miles,  h^^is  the  height  of  the  transmitter  antenna  in 
meters  and  is  the  height  of  the  receiver  antenna  in  meters. 

For  seagoing  vessels,  their  dimensions  are  defined  in  length,  breadth,  draught  and 
tonnage.  The  height  above  sea  level  or  air  draught,  as  defined  in  Naval  Architecture 
terms,  is  seldom  specified.  To  approximate  the  height  of  the  CMR  above  sea  level,  the 
new  PANAMAX  standard  states  that  all  ships  passing  through  the  Suez  Canal  must  not 
have  an  air  draught  of  more  than  68  m,  with  the  restriction  set  by  the  Suez  Canal  bridge. 
It  is  reasonable  to  assume  the  likes  of  oil  tankers  and  containers  ship  on  international 
voyage  are  designed  to  conform  to  this  requirement.  For  a  commercial  vessel,  the  CMR  is 
normally  installed  at  the  highest  point.  On  the  other  end,  the  height  of  the  EW  receiver  on 
the  frigate  is  approximated  to  be  at  20  m  above  sea  level.  Therefore,  R„„is  calculated 
from  (1)  to  be  52  km. 
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There  are  many  international  companies  like  the  UK-based  Kelvin  Hughes  and 
Japanese-based  Furuno  that  supply  the  world’s  shipping  industry  with  integrated  bridge 
control  systems  along  with  CMR  and  MAIS.  CMR  such  as  Kelvin  Hughes’s 
MantaDigital  Radar  (See  Appendix  A  for  data  sheet)  and  Furuno’s  FAR28x7  series  (See 
Appendix  B  for  data  sheet)  are  the  most  commonly  installed  onboard  larger  class  vessels. 
Out  of  all  the  technical  specifications,  a  few  parameters  are  needed  for  the  computation 
of  the  SNR  for  detection.  The  specifications  and  comparison  between  the  two  CMRs  are 
presented  in  Table  2. 


Table  2.  Comparison  of  two  of  most  commonly  installed  CMR. 


Technical 

MantaDigital 

FAR28x7  series 

Specifications 

S-Band 

X-band 

S-band 

X-band 

Pulse  Width 

55  ns  (Short), 

250  ns  (Med), 

950  ns  (Long) 

70  ns  (Short), 
250  ns  (Med), 
900  ns  (Long) 

0.5  ps  (Short), 
0.7  ps  (Med), 

1.2  ps  (Long) 

0.5  ps  (Short), 

0.7  ps  (Med), 

1.2  ps  (Long) 

Beamwidth 

(Vertical) 

26° 

- 

25° 

20° 

Beamwidth 

(Horizontal) 

1.9° 

1.25° 

1.8° 

1.23° 

Pulse 

Repetitive 

Frequency 

3000  Hz  (Short), 
750  Hz  (Med), 

750  Hz (Long) 

3000  Hz  (Short), 
1500  Hz  (Med), 
750  Hz  (Long) 

1000  Hz 
(Short), 

600  Hz  (Med), 
600  Hz  (Long) 

1000  Hz  (Short), 
600  Hz  (Med), 
600  Hz  (Long) 

Receiver 
Noise  Figure 

<6dB 

<6dB 

- 

- 

Antenna 

Gain 

28  dB 

30  dB 

- 

- 

Transmitter 

Power 

30  kW 

25  kW 

30  kW 

30  kW 

Frequency 

3050  MHz  +!-  10 
MHz 

9410  MHz -1-/-  30 
MHz 

3050  MHz  +!- 
30  MHz 

9410  MHz  +!-  30 
MHz 

Scan  Rate 

24  rpnV  45  rpm 

24  rpm/  45  rpm 

24  rpm/  42  rpm 

21  rpm/  26  rpm/ 
45  rpm 

In  reference  to  the  technical  specifications  shown  in  Table  2,  both  CMRs  share 
similar  parameters,  but  not  all  parameters  of  the  FAR-28x7  series  are  available  in  its 
product  catalog.  For  the  purpose  of  this  thesis,  the  technical  specifications  of  the 
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MantaDigital  Radar  are  used  to  compute  the  SNR  and  maximum  detection  range  of  the 
PER. 


After  selecting  the  technical  specifications  of  the  opportunistic  radar,  the  next  step 
is  to  obtain  an  accurate  position  of  all  available  opportunistic  radars  so  as  to  calculate  the 
baseline  distance  between  the  transmitter  and  EW  receiver  on  the  warship.  Furthermore, 
the  receiver  needs  to  be  able  to  uniquely  classify  each  CMR  transmission  and  associate  it 
with  the  source  vessel  on  the  CMS.  From  the  operating  specifications  requirement  stated 
by  the  IMO,  the  MAIS’s  function  is  to  support  automated  tracking  and  identification 
among  vessels  and  port  authorities.  It  operates  autonomously  and  continuously  to 
exchange  ship  identity  and  GPS  locations  using  the  VHF  transceiver.  There  are  two  types 
of  MAIS:  Class  A  installed  onboard  a  majority  of  seagoing  commercial  vessel,  and  the 
more  economical  Class  B  used  typically  by  leisure  and  smaller  crafts.  The  MAIS 
communication  network  uses  self-organizing  time  division  multiple  access  (SOTDMA) 
to  allow  for  at  least  4500  communication  timeslots  for  vessels  to  exchange  up  to  27 
different  National  Marine  Electronics  Association  (NMEA)  -  0183  messages,  sorted  into 
different  level  of  priorities.  The  highest  priority  message  is  the  GPS  positional 
information  accompanied  by  a  unique  nine-digit  identification  number  assigned  to  each 
vessel  known  as  the  maritime  mobile  service  identity  (MMSI).  This  information  allows 
the  warship  to  discern  targets  of  interest  from  neutrals  on  its  CMS  display.  Extracted 
from  the  U.S.  Coast  Guard  website,  the  key  difference  in  specifications  between  Class  A 
and  Class  B  MAIS  are  summarized  in  Table  3. 
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Table  3.  Comparison  of  specifications  between  Class  A  and  Class  B 

MAIS. 


Shipboard  AIS 

Class  A 

Class  B/Self- 
Organizing 

Class  B/Carrier 
Sense 

Transmit  Power 
(Watts) 

12.5  W/ 2  W 

5  W/2  W 

2  W 

Primary  Access 
Scheme 

SOTDMA 

SOTDMA  Carrier- 

sense 

TDMA 

non-competing  with 
SOTDMA  units 

Position  Reporting 
Rate 

Either  every  2,  3  Vi,  6 
or  10  s  based  on 
speed  and  course 
change.  Every  3  min 
when  <  3  kts 

Either  every  5,  15  or 

30  s 

based  on  speed  Every 

3  min  when  <  2  kts 

Every  30  s 

Every  3  min.  when  < 

2  kts. 

Static  Data 
Reporting  Rate 

Every  6  min 

Every  6  min 

Every  6  min 

Transmit  Data 

All 

No  Rate  of  Turn, 
Navigation  Status, 
Destination,  ETA, 

No  Rate  of  Turn, 
Navigation  Status, 
Destination,  ETA, 

From  the  U.S.  Coast  Guard  website,  the  worst-case  GPS  positional  accuracy  of 
the  MAIS  is  0.001  minute  in  both  latitude  and  longitude  with  an  update  rate  of  no  longer 
than  10  s  for  a  Class  A  MAIS.  That  translates  to  a  positional  error  of  not  more  than  2  m 
all  round.  With  reference  to  the  MantaDigital  Radar  chosen  earlier,  the  same  class  of 
vessel  is  equipped  with  a  Class  A  MAIS.  For  the  purpose  of  this  thesis,  the  information 
provided  by  Class  A  MAIS  is  used  to  calculate  the  baseline  distance  and  the  detection 
accuracy  of  the  target. 

B.  EW  AND  DF  RECEIVERS 

Onboard  a  warship,  the  EW  receiver  is  responsible  for  the  extraction  of  key 
parameters  from  the  intercepted  transmission  for  follow-on  classification  and 
identification  typically  through  digital  radio  frequency  memory  (DRFM).  The  EW 
receiver  is  capable  of  measuring  the  intercepted  transmission’s  center  frequency,  pulse 
width,  pulse  repetitive  frequency  (PRE),  antenna  scan  rate,  type  of  radar,  transmission 
mode,  etc.  Integrated  with  the  DE  receiver,  the  angle-of-arrival  can  be  determined  and 
used  to  point  to  the  source  of  transmission  with  the  aid  of  a  bearing  line.  Parameters  like 
these  are  saved  in  the  DREM  to  help  with  future  identifications. 
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The  important  parameters  to  quantify  on  the  EW  and  DF  receiver  are  the  lowest 
received  power  that  can  be  detected,  known  as  sensitivity,  bandwidth  of  the  receiver, 
detection  bandwidth  and  DF  accuracy.  Due  to  the  classified  nature  of  these  parameters, 
they  are  seldom  available  in  open  source  literature.  As  such,  the  parameters  are  based  on 
systems  with  published  specifications  that  best  matched  our  scenario  and  as  practical  as 
reasonable. 

The  EW  receiver  specifications  are  based  on  UK-based  Teledyne  Defense  Ftd.’s 
high  performance  receiver  model  QR026  (See  Appendix  C  for  data  sheet).  This  wideband 
receiver  operates  from  0.2  GHz  to  18  GHz  and  is  expandable  to  frequencies  up  to  40 
GHz.  It  has  a  sensitivity  of  approximately  -70  dBm,  a  typical  noise  figure  of  13.5  dB, 
and  a  receiver  gain  of  20  dB  if  operating  with  a  scanning  bandwidth  of  500  MHz.  Most 
EW  receivers  support  tuning  of  bandwidth  to  enhance  detection  of  specific  frequency 
band  and  the  use  of  a  notch  or  Yttrium  iron  garnet  (YIG)  filter  to  carry  out  narrowband 
detection  and  classification.  For  our  case,  the  bandwidth  of  the  receiver  can  tuned  to 
coincide  with  the  frequency  band  of  interest,  thus  reducing  additive  white  gaussian  noise 
(AWGN),  thereby,  improving  the  sensitivity  and  SNR. 

The  DF  receiver  specifications  are  based  on  South  Africa’s  Poynting  Defense’s 
land-based  DF  Antenna  model  DF-A0062  (See  Appendix  D  for  data  sheet),  which  has  an 
effective  DF  range  from  20  MHz  to  6  GHz  and  U.S. -based  DRS  Defense  Solutions 
shipboard  DF  Antenna  array  model  ZA-4501.  Both  DF  antenna  systems  work  on  a  five- 
element  phase  interferometer  concept  with  an  azimuth  DF  accuracy  of  between  1°  to  2°. 
For  the  purpose  of  this  thesis,  the  average  DF  accuracy  of  1.5°  is  used. 

C.  BISTATIC  RCS 

Bistatic  RCS  of  a  target  is  dependent  on  the  angle  formed  between  the  direction  of 
incidence  and  the  direction  of  reception,  known  as  the  bistatic  angle  p.  Unlike  the 
monostatic  RCS,  which  is  a  measure  of  the  amount  of  EM  energy  re-radiated  back  in  the 
direction  of  incidence,  the  bistatic  RCS  is  the  measure  of  the  amount  of  EM  energy 
returned  in  the  direction  of  the  receiver  that  is  not  co-located  with  the  transmitter.  Putting 
everything  in  perspective,  there  is  only  one  set  of  360  values  in  monostatic  RCS 
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corresponding  to  360  degrees  in  azimuth  for  each  target.  However,  in  the  case  of  the 
bistatic  RCS,  each  direction  of  incidence  has  its  own  set  of  RCS  values  for  360  degrees. 
So  to  completely  map  out  the  bistatic  RCS  of  a  target,  there  are  360  sets  of  values  per 
target. 


The  definition  of  RCS  can  be  stated  [2]  as 


a  =  lim  AnR 

/?— ^00 


2  * 


(2) 


For  simple  geometrical  shapes,  the  RCS  can  be  easily  calculated  using 

4;rAi 


cr  = 


2' 


(3) 


where  is  the  equivalent  area  in  m^  and  2  is  the  wavelength  in  meters.  In  the  case  of  a 


complex- shaped  low-RCS  warship,  the  approximation  of  the  RCS  has  to  take  into 
consideration  the  ‘electrical  size’  of  the  target,  which  is  defined  by  the  largest  relative 
dimension  of  the  target  in  terms  of  the  incident  wavelength.  There  are  three  frequency 
regions  in  which  the  RCS  of  a  target  can  vary  significantly,  namely  the  Rayleigh  region 
(low  frequency),  Mie  region  (resonance)  and  Optical  region  (high  frequency).  The  three 
regions  are  classified  based  on  the  magnitude  of  the  ‘electrical  size.’  The  graph  (Figure  5) 
presents  the  effect  on  the  monostatic  RCS  by  varying  the  wavelength  on  the  same  sphere. 
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Sphere  circumierence  in  wavelengths 


Figure  5.  Effects  on  monostatic  RCS  of  a  sphere  with  changes  to  the  wavelength. 

From  [1]. 

Due  to  the  profound  effect  of  the  ratio  of  the  target  length  versus  the  input 
transmission  wavelength,  different  methods  of  RCS  prediction  are  used  in  response  to  the 
dominant  scattering  mechanism  that  is  present.  The  RCS  prediction  methods  most 
commonly  used  are  physical  optics,  microwave  optics,  method  of  moments,  the  finite 
element  method  and  the  finite  difference  method.  For  our  case,  the  low-RCS  warship  is 
chosen  to  be  the  Republic  of  Singapore  Navy’s  (RSN)  Formidable-class  frigates.  A  three- 
dimensional  model  (Figure  6  and  7)  of  the  Formidable-class  frigate  is  created  in  CST 
Studio  and  ported  over  to  FEKO,  a  commercial  RCS  modeling  software,  to  calculate  the 
free  space  monostatic  and  bistatic  RCS  at  S-band  and  X-band  using  physical  optics.  The 
physical  optics  method  is  suited  for  use  to  predict  the  RCS  for  targets  that  are  in  the 
optical  frequency  region. 
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Figure  6.  Actual  photo  of  RSN’s  Formidable-class  stealth  frigate  (top,  from  [18])  versus 

its  model  (side  profile). 

A  quick  comparison  between  the  monostatic  RCS  and  bistatic  RCS  shows  that  the 
RCS  differences  can  fluctuate  up  to  a  factor  of  30  dB  over  an  aspect  change  of  a  few 
degrees  as  illustrated  in  Figure  8.  A  sample  set  of  bistatic  RCS  plots  are  attached  as 
Appendix  E,  and  the  values  are  added  to  the  MATLAB  simulation  in  the  area  of  detection 
coverage. 
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Figure  7.  Actual  photo  of  RSN’s  Formidable-class  stealth  frigate  (top,  from  [19])  versus 

its  model  (front  profile). 


Figure  8.  Comparison  of  monsotatic  RCS  (dBsm)  of  the  model  measured  at  3  GHz  (left) 
and  bistatic  RCS  with  incidence  angle  at  3  degrees  (right). 
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D.  BISTATIC  RADAR  RANGE  EQUATION 


In  reference  [1],  the  bistatic  radar  range  equation  is  derived.  The  single  range  in 
the  monostatic  case  is  now  split  into  two  range  variables  Rj  and  ,  and  the  monostatic 


RCS  is  replaced  with  the  bistatic  RCS  denoted  by  cr^  as 


(j gFj 


(4) 


where 


Rj  =  transmitter-to-target  range, 

Rff  =  target-to-receiver  range, 

Pj  =  transmitter  power  output, 

Gj  =  transmitter  antenna  power  gain, 

=  receiving  antenna  power  gain, 

X  =  wavelength, 
crg=  bistatic  RCS, 

Fj.  =  pattern  propagation  factor  for  transmitter-to-target  path, 

=  pattern  propagation  factor  for  target-to-receiver  path, 
k  =  Boltzmann’s  constant, 

=  receiving  system  noise  temperature, 

B^=  noise  bandwidth  of  receiver’s  pre-detection  filter,  sufficient  to  pass  all 
spectral  components  of  the  transmitted  signal, 

(5/A)min  =  signal-to-noise  power  ratio  required  for  detection, 

Lj  =  transmitting  system  losses  (>  1)  not  included  in  other  parameters, 

=  receiving  system  losses  (>  1)  not  included  in  other  parameters,  and 
K  =  bistatic  maximum  range  product. 


By  grouping  some  of  the  terms  listed,  Eq.  (4)  can  be  presented  in  a  more  compact 
and  generalized  form  for  any  SNR  value.  Defining 

^  ^  P,G,G  A^a.F^^F^ 

(4;rfkT^B^L,L, 

we  can  write,  in  general. 
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K 


(6) 


S/N^ 


p2  p2 


The  SNR  value  is  only  sensitive  to  parameters  associated  with  the  CMR  and  EW 
receiver.  It  is  not  grounded  in  any  coordinate  system  or  locations  of  the  transmitter, 
receiver  and  target. 

From  geometry  (refer  to  Figure  1), 

(7) 

2{R^+R^+L^mO^) 

where  L  is  the  baseline  distance  between  the  transmitter  and  receiver,  and  is  the  angle 
between  North  (y-axis)  and  the  target. 

In  the  real  world,  target  detection  begins  with  the  measurement  of  {Rj  +Rg)  using 
one  of  two  methods,  namely  the  direct  method  which  requires  a  line-of-sight  between  the 
transmitter  and  receiver  or  the  indirect  method  which  requires  synchronization  of  stable 
clocks  between  the  transmitter  and  receiver.  Given  the  use  of  opportunistic  transmitters, 
the  direct  timing  method  is  used  here  and 


{Rj+Rj^)-c^T^+L  (8) 

where 

c  =  speed  of  light,  and 

Ar„=  time  interval  between  the  reception  of  the  transmitted  pulse  and 
reception  of  the  target  echo. 

In  this  case,  the  EW  receiver  is  the  one  providing  this  information.  The  EW 
receiver  must  first  characterize  the  transmission  of  each  CMR  for  quick  recognition. 
Upon  receiving  the  direct  transmission  from  the  CMR  as  verified  by  the  DF  receiver,  the 
timer  starts  counting  until  it  receives  the  same  CMR  transmission  from  a  different 
bearing.  This  time  difference  is  .  The  measured  difference  in  bearing  is  then  used  to 

determine  .  With  that,  the  detection  coverage  based  on  SNR  can  be  plotted.  The  plot  is 
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then  normalized  using  the  SNR^^to  display  the  detection  coverage  supportable  by  the 
EW  receiver.  The  SNR^^  (in  dB)  is  calculated  from 

Sensitivity  =  SNR^^  +  101og(B„ )  +  - 174  (9) 


As  shown  in  Eqs.  (6)  through  (8),  the  SNR  value  corresponding  to  a  specific 
combination  of  Tx-Tgt-Rx  can  be  calculated  given  all  the  parameter  values.  Putting  all 
these  individual  SNR  contours  together  on  a  Cartesian  coordinate  plot,  we  find  the 
detection  coverage  for  a  particular  combination  of  PBR.  Given  the  high  availability  of 
opportunistic  transmitters  out-at-sea,  multiple  combinations  of  PBR  can  be  used  for 
optimum  coverage  in  real-time  to  allow  for  the  highest  probability  of  detection. 


E.  MEASUREMENT  AND  LOCATION  ERROR 


Prom  [1],  there  are  three  geometry-dependent  measurement  and  location  errors  in 
a  PBR  system,  namely  the  receiver-to-target  error  due  to  measurement  of(R^-i-R^), 


baseline  measurement  error  of  L  and  receiver  look  angle  error  due  to  DP  error  by  the  DP 
receiver.  The  receiver-to-target  error  is  the  result  of  transmitter  beamwidth  spread  at  the 
target  at  R^  where  the  edge  of  azimuth  beamwidth  is  strong  enough  to  create  returns  back 


to  the  receiver.  The  baseline  measurement  error  is  dependent  on  the  means  by  which  the 
measurement  of  the  baseline  is  obtained.  In  our  case,  the  baseline  L  is  computed  as  the 
difference  between  the  differential  GPS  locations  of  the  transmitter  and  receiver.  As  such, 
the  measurement  error  is  a  function  of  the  accuracy  of  the  GPS  position  fix.  The  root- 
sum-squared  (rss)  of  the  three  errors  in  [1]  is  expressed  as 


dR*  =i 


dR. 


d(Rj  +  Rr  ) 


d  (Rr  -i-  Rr  ) 


dL 


dL 


dOg 


1/2 


(10) 


Depending  on  the  geometry  of  Tx-Tgt-Rx,  the  three  errors  are  still  additive  but 
with  components  on  different  planes,  as  shown  in  Pigures  9  through  11.  With  changing 
geometry,  the  combined  error  forms  an  ellipse  with  a  different  major  and  minor  axis. 
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Tgt 


Rx 


Figure  9.  Illustration  of  measurement  error  due  to  variation  of  L. 


Tgt 


Figure  10.  Illustration  of  measurement  error  due  to  transmitter’s  . 


Tgt 


Figure  11.  Illustration  of  measurement  errors  due  to  receiver’s  DF  error. 

The  three  components  of  error  contributing  to  the  rss  in  Eq.  (10)  are 

dR^  _  (Rj  +Ri^f  +l3  +  2L(Rj.  +  )  sin  6^ 

d(^Rj  +  Rj^ )  2(^Rj  +  +  L  sin  0j^ ) 

_l  +  e^  +  2e  sin  0^. 

2(l  +  esin6>^)^ 


(11) 
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(12) 


dRj^  [L  +  {Rj  +  R/f )  ]  sin  0^^  +  '2.L,{Rj  +  ) 

dL  ^{Rj  +  ^  sin  6*^ 

_  {e^  + 1)  sin  +  le 

2(l  +  esin^^)^ 
and 

dR^  _  U{Rj  +R^f-l3]co^ 6^ 

90^  2,(^Rj  +  +  L  sin  0^ ) 

_  L(l-e^)cos6'^ 

2(l  +  esin6'^)^ 

where  e  =  eccentricity  given  by  L/(Rj,+Rg). 

To  begin  using  this  PBR  system  in  a  real-world  scenario,  one  has  to  first 
determine  what  CMRs  (opportunistic  radars)  are  within  range  and  what  combination  of 
them  offers  the  best  overall  detection  coverage  and,  if  not,  at  a  minimum,  optimize 
coverage  within  their  operating  areas. 

With  the  warship’s  MAIS  in  the  receive-only  mode,  all  merchant  ships  that  are 
within  the  VHP  reception  range  are  displayed  on  the  CMS  screen.  Based  on  the  returns 
from  the  MAIS,  the  merchant  ships  that  operate  with  Class-A  MAIS  are  selected,  and 
their  S-band  and/or  X-band  CMRs’  parameters  are  collected  and  analyzed  using  the  EW 
receiver.  The  parameters  are  then  saved  into  the  DRFM  and  associated  with  the  bearing 
relative  to  the  warship  (receiver)  given  by  the  DF  receiver.  The  Class-A  MAIS  of  all 
selected  merchant  ships  are  constantly  updating  their  GPS  positions  at  an  update  rate  of 
no  greater  than  10  s.  With  the  unique  identifier  provided  by  the  MAIS,  the  specifications 
of  the  CMR  can  be  determined.  Putting  all  this  information  together,  the  crews  onboard 
the  warship  have  all  the  necessary  information  for  Eqs.  (5)  and  (8)  to  plot  the  detection 
coverage  for  a  single  PBR  system.  All  the  involved  parameters  tagged  to  their  individual 
sources  are  summarized  in  Table  4.  A  further  extension  of  this  is  to  establish  the 
relationship  between  different  Tx-Rx  PBR  pairs  and  their  associated  detection  coverage. 
The  warship  can  choose  multiple  PBR  pairs  and  combine  their  individual  detection 
coverage  plots  to  form  a  more  comprehensive  coverage  in  the  area  of  interest.  Different 
constellations  of  Tx-Rx  pairs  offer  different  advantages  and  disadvantages.  Knowing  this 
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in  advance,  along  with  practice,  allows  the  warship  operators  to  pick  the  best 
constellations  from  the  prevailing  availability  of  opportunistic  transmitters  in  a  minimum 
time. 


Table  4.  Parameters  summary  for  Eqs.  (5)  and  (8). 


Systems 

MantaDigital 

(CMR) 

Class-A 

(MAIS) 

QR026 

(EW 

Receiver) 

DF-A0062 

(DF 

Receiver) 

FEKO 

RCS 

Modeling 

Software 

Remarks 

Pt 

30  kW 
(S-band) 

25  kW 
(X-band) 

- 

- 

- 

- 

- 

Gj 

28  dB 
(S-band) 

31  dB 
X-band) 

- 

- 

- 

- 

- 

Gr 

- 

- 

2dB 

- 

- 

- 

X 

- 

- 

0.0984  m 
(S-band) 
0.0319  m 
(X-band) 

- 

- 

- 

See 

Appendix 

E 

Ft 

- 

- 

- 

- 

- 

1 

See 

reference 

[4] 

Fr 

Ts 

- 

- 

3915  K 
(based  on 
noise  figure 
of  13.5  dB) 

- 

- 

- 

B„ 

- 

- 

500  MHz 

- 

- 

- 

Lt 

- 

- 

- 

- 

- 

2 

See 

reference 

[4] 

Lr 

L 

- 

Based  on 
GPS 

location  of 
Tx  and  Rx 

- 

- 

- 

- 

AT„ 

- 

- 

Based  on 
GPS  location 
of  Tx  and  Rx 

- 

- 

- 

- 

- 

- 

Bearing 

Fine 

- 

- 
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To  plot  an  operationally  usable  detection  plot,  the  minimum  SNR  with  reference 
to  the  sensitivity  of  the  EW  receiver  has  to  be  determined  using  Eq.  (9).  This  defines  the 
value  of  SNR  beyond  which  reception  of  the  returns  from  the  target  is  possible.  Based  on 
a  sensitivity  value  of  -70  dBm,  the  SNRmin  required  in  our  scenario  is  more  than  3.5  dB 
for  single  pulse.  To  allow  for  unaccounted  losses  or  unforeseen  variations,  an  additional 
buffer  is  added  to  increase  the  SNRmin  to  10  dB  in  our  simulation. 

In  situations  where  signal  strength  of  the  returned  signal  is  weak,  non-coherent 
pulse  integration  across  multiple  pulses  can  improve  the  chance  of  detection  by 

approximately  yJW  dB,  where  N  is  the  number  of  pulses  used  in  the  integration  [19].  This 
approximation  is  heavily  dependent  on  the  type  of  receiver  architecture  used  and  is  more 
applicable  for  envelope  detector  receiver  types.  The  number  of  pulses  transmitted 
downrange  can  be  calculated  using 

N  = - ^ - (PRF)  (14) 

360(scanrate) 


where 

0^^  =  Azimuth  beamwidth  of  the  CMR  in  degrees, 

scanrate=  CMR  Antenna  scan  rate  in  revolution  per  second  (rps), 
PRF  =  CMR’s  antenna  PRE  in  Hz. 


In  the  case  of  the  S-band  KH  MantaDigital  radar,  the  number  of  pulses  that  are 
available  for  integration  per  second  is  five  pulses  with  ^^2  of  1-9°,  scanrate  of  45  rpm  and 
PRF  of  750  Hz.  This  translates  to  a  gain  of  2.2  dB.  Thus,  by  allowing  for  more 
integration  time,  it  is  possible  to  increase  detection  range.  The  low  SNR  of  each  pulse 
and  the  long  integration  time  also  bring  with  them  increased  measurement  errors, 
specifically,  in  the  TDOA  measurements  of  the  leading  edge  of  the  returned  pulse.  The 
increase  in  error  is  discussed  in  [20]  for  a  monostatic  radar,  but  the  application  to  PBR 
adopting  the  discussed  form  of  TDOA  is  not  well-defined. 
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F.  DETECTION  ACCURACY 

As  the  PER  system  detects  targets  passively  and  displays  them  on  the  CMS 
screen,  a  problem  of  detection  accuracy  arises.  As  shown  in  Eqs.  (13)  through  (15),  there 
are  measurement  and  location  errors  that  form  an  uncertainty  ellipse  around  the  target. 
On  top  of  these  baseline  errors,  there  are  three  sources  of  errors  that  add  to  the  growing 
uncertainty  ellipse.  They  are:  (1)  azimuth  beam  spread  of  the  CMR  with  target  range,  (2) 
accuracy  of  the  DF  receiver  and  (3)  positional  errors  due  to  data  age  and  moving  targets. 
These  sources  of  errors  are  added  in  their  respective  planes  and  determine  the  uncertainty 
ellipse  for  each  Tx-Rx  pair.  By  knowing  the  passively-detected  target  range  and  azimuth 
with  respect  to  the  warship,  one  can  compute  the  worst-case  location  errors  in  Cartesian 
coordinates.  The  location  errors  can  be  reduced  by  using  more  Tx-Rx  pairs  to  narrow 
down  the  overlapped  area. 

In  the  next  chapter,  the  simulation  of  PER  and  PMR  are  presented.  Different 
opportunistic  transmitter  placements  are  simulated  to  understand  their  effects  on 
detection  coverage  and  detection  accuracy. 
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III.  SIMULATION  RESULTS  FOR  DETECTION  COVERAGE 


In  this  chapter,  the  detection  coverage  of  the  PBR  is  investigated  with  different 
combinations  of  Tx-Rx  constellations.  The  chapter  begins  with  a  single  transmitter  and 
receiver  pair  and  demonstrates  how  a  PBR  system  has  a  higher  probability  of  detecting  a 
low-RCS  target.  The  simulation  results  also  show  the  limited  detection  coverage  of  a 
single  PBR  system.  The  investigation  goes  on  to  explore  a  multiple  PBR  system  and 
evaluates  the  detection  performance  of  the  combined  system  using  different  Tx-Rx 
constellations.  All  simulation  results  are  generated  using  MATLAB  with  the  source  codes 
included  in  Appendix  F. 

A.  SINGLE  TX-RX  PAIR 

To  appreciate  the  benefits  brought  about  by  combining  multiple  PBRs  into  a 
multistatic  radar  system,  one  first  needs  to  understand  the  capabilities  and  limitations  of  a 
PBR  system  operating  with  S-band  and  X-band  CMR. 

I.  Monostatic  RCS  with  S-Band  CMR 

According  Skolnik  [14],  at  low  grazing  angle,  the  median  RCS  of  a  ship  can  be 
estimated  using  the  empirical  formula 

cj  =  524fD^'^ 

with  the/in  MHz  and  ship’s  displacement  D  in  kilotons.  Substituting /as  3050  MHz  and 
D  as  3.2  kilotons,  the  displacement  of  RSS  Formidable-class  frigate,  we  find  the  RCS  of 
a  non-RCS  treated  frigate  to  be  5137  m  .  A  separate  simulation  using  FEKO  on  the 
model  based  on  RSN’s  Formidable-class  frigate  yields  a  more  representative  monostatic 
RCS  plot  (Figure  12).  The  monostatic  RCS  fluctuates  from  0  dBsm  to  72  dBsm 
dependent  on  the  incidence  angle.  Using  the  results  shown  in  Figure  12,  we  estimated  the 
simulated  median  RCS  as  18  dBsm  (63  m  )  for  S-band. 
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Figure  12.  Monostatic  RCS  of  RSN’s  Formidable-class  frigate  simulated  at  3050  MHz 

using  FEKO. 

Putting  both  median  monostatic  RCS  values  into  simulations,  we  found  the 
resulting  detection  coverage  of  a  PER  with  a  S-band  CMR  at  (-5000,0)  as  shown  in 
Figures  13  and  14.  These  plots  are  generated  using  Matlab  scatter  plot  function.  Each 
data  point  is  shown  as  a  colored  dot.  The  density  of  the  data  changes  from  plot  to  plot. 
Note  that  gaps  between  the  individual  dots  on  the  low  density  plots  are  not  areas  of  low 
SNR,  which  has  been  set  at  10  dB. 


Eigure  13.  SNR  detection  plot  (dB)  using  empirical  median  RCS  at  S-band  CMR  with  Rx 

at  the  origin  and  single  Tx  at  (-5000,0). 
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Figure  14.  SNR  detection  plot  (dB)  using  simulated  median  RCS  at  S-band  CMR  with 

Rx  at  the  origin  and  single  Tx  at  ( -5000,0). 

With  the  minimum  SNR  set  at  10  dB  and  the  EW  receiver  bandwidth  of  500  MHz 
over  an  integration  time  of  1  s,  the  radial  detection  coverage  under  such  conditions 
extends  to  approximately  7500  m  from  the  receiver  with  near  100%  coverage.  From 
Figure  14,  the  detection  coverage  with  the  use  of  simulated  median  RCS  shows  a  drop  in 
detection  coverage  range  to  that  of  1500  m  from  the  receiver.  The  observable  gaps  in  the 
plot  of  Figure  14  are  the  result  of  plotting  resolution  and  not  actually  detection  gaps.  The 
detection  void  between  the  transmitter  and  receiver  is  due  to  the  limitation  offered  by  the 
PBR  geometry.  As  the  target  approaches  the  baseline  formed  by  the  transmitter  and 
receiver  for  the  TDOA  calculation  expressed  in  Eq.  (8),  the  difference  between  the  direct 
path  L  and  the  reflected  path  (Rj,+Rj^)  converges  to  zero.  All  equations  expressed  in  this 

thesis  break  down  both  in  theory  and  in  practice  since  the  target  is  now  directly  between 
the  transmitter  and  receiver  line-of-sight.  In  such  cases,  the  PBR  functions  as  a  forward 
scatter  radar  instead  of  a  PBR. 
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2. 


Monostatic  RCS  with  X-Band  CMR 


Using  the  Eq.  (15),  we  can  show  the  monostatic  RCS  of  a  typical  frigate  at  X- 
band  (9410  MHz)  to  be  28,875  m  .  A  separate  simulation  using  FEKO  on  the  model 
based  on  RSN’s  Eormidable-class  frigate  yields  a  more  representative  monostatic  RCS 
plot  (see  Figure  15).  The  monostatic  RCS  fluctuates  from  -10  dBsm  to  80  dBsm 
dependent  on  the  incidence  angle.  Using  the  results  shown  in  Figure  15,  we  find  the 
estimated  simulated  median  RCS  to  be  15  dBsm  (32  m  )  for  X-band. 


Figure  15.  Monostatic  RCS  of  RSN’s  Formidable-class  frigate’s  model  simulated  at  9410 

MHz  using  FEKO. 

For  the  purpose  of  comparison,  we  put  both  median  monostatic  RCS  values 
(median  from  formula  and  median  from  FEKO)  into  the  simulation.  The  resulting 
detection  coverage  of  a  PBR  with  an  X-band  CMR  at  (-5000,0)  is  shown  in  Figures  16 
and  17.  The  magnitude  of  the  SNR  behaves  in  the  same  way  for  all  PBR  systems  with 
high  SNR  around  the  vicinity  of  the  transmitter  and  receiver  and  reduces  in  a  linearly 
radial  manner.  This  simulation  is  consistent  with  findings  by  [10]. 
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Figure  16.  SNR  detection  plot  (dB)  using  empirical  median  monstatic  RCS  from  Eq.  (15) 
at  X-band  CMR  with  Rx  at  the  origin  and  single  Tx  at  (-5000,0). 


Figure  17.  SNR  detection  plot  (dB)  using  simulated  median  monostatic  RCS  from  FEKO 
at  X-band  with  Rx  at  the  origin  and  single  Tx  at  ( -5000,0). 


The  low  FEKO  simulated  monostatic  RCS  has  significantly  reduced  the  detection 
coverage  to  be  restricted  to  regions  around  the  transmitter  and  receiver.  This  indicates 
that  without  using  a  PBR,  the  ability  to  detect  a  low-RCS  target  monostatically  is 
severely  limited. 
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3. 


Bistatic  RCS  with  S-Band  CMR 


With  the  inclusion  of  bistatic  RCS  in  the  simulation,  the  detection  coverage  varies 
with  the  change  in  target  aspect  with  respect  to  the  transmitter  and  receiver.  The  effect  of 
the  bistatic  RCS  is  most  evident  with  changes  made  to  bistatic  angle  by  changing  the 
course  and  heading  of  the  target.  Many  possible  combinations  of  P  with  different  bistatic 
geometry  of  Tx-Tgt-Rx  exist.  The  most  representative  are  those  with  the  target’s  general 
direction  headed  towards  the  receiver  and  the  other  with  the  target’s  broadside  directed 
towards  the  receiver  on  a  crossing  course.  The  former  presents  one  of  the  smallest  bistatic 
RCS  returns,  while  the  latter  presents  a  larger  bistatic  RCS  return.  With  the  ship  forward 
and  aft  line  along  the  0°  and  180°,  it  can  be  seen  from  Figure  18  that  the  bistatic  returns 
are  mostly  small  with  a  narrow  peak  occurring  at  the  negative  of  the  incidence  angle  but 
always  higher  than  the  direct  monostatic  returns. 

The  RCS  data  was  generated  in  1°  increments.  When  used  in  the  PBR 
simulations,  the  incidence  and  observation  angles  were  rounded  to  the  nearest  degree  in 
order  to  obtain  the  RCS. 


Figure  18.  Bistatic  RCS  of  RSN’s  Formidable-class  frigate’s  model  simulated  at  3050 

MHz  using  FEKO. 
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A  comparison  of  Figures  17  and  19  shows  clearly  the  difference  in  detection 
coverage  when  a  different  aspect  of  the  low-RCS  target  is  presented  to  the  transmitter  and 
receiver.  The  detection  range  increases  from  3000  m  to  that  of  beyond  20,000  m  with  a 
change  of  target  aspect.  It  is,  therefore,  important  to  understand  a  single  PER  system  is 
also  limited  in  its  ability  to  ensure  consistent  detection  coverage  and  low  detection  gaps. 
The  resultant  detection  coverage  is  dependent  on  the  Tx-Tgt-Rx  geometry  and  the  aspect 
of  the  target. 


Distance  in  m 


Figure  19.  SNR  detection  plot  (dB)  using  bistatic  RCS  at  S-band  having  reflections  off 
the  forward  aspect  of  the  target  ship  by  Txl  (-5000,0). 
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Figure  20.  SNR  detection  plot  (dB)  using  bistatic  RCS  at  S-band  having  reflections  off 
the  broadside  of  the  target  ship  by  Txl  (-5000,0). 


Due  to  the  large  comparative  difference  in  the  ship’s  dimensions  with  respect  to 
the  wavelength  of  the  transmission,  a  higher  resolution  plot  of  the  target  RCS  reveals 
more  RCS  fluctuations  between  the  current  sampling  intervals  of  one  degree. 

4.  Bistatic  RCS  with  X-Band  CMR 

Similar  conclusions  can  be  drawn  from  the  detection  coverage  plot,  as  shown  in 
Figures  21  and  22,  with  an  X-band  CMR.  In  general,  the  X-band  CMR  as  an 
opportunistic  transmitter  has  a  smaller  detection  coverage  as  compared  to  the  S-band 
CMR.  The  detection  gaps  are  also  more  pronounced.  In  a  real-world  scenario,  the  use  of 
RAM  becomes  more  practicable  at  small  wavelengths.  It  follows  that  a  low-RCS  target  is 
likely  to  have  even  lower  RCS  returns  at  higher  frequencies,  reducing  the  usefulness  of 
X-band  CMR  as  an  opportunistic  transmitter. 
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Figure  21.  SNR  detection  plot  (dB)  using  bistatic  RCS  at  X-band  having  reflections  off 
the  forward  aspect  of  the  target  ship  by  Txl  (-5000,0). 


Figure  22.  SNR  detection  plot  (dB)  using  bistatic  RCS  at  X-band  having  reflections  off 
the  broadside  of  the  target  ship  by  Txl  (-5000,0). 
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B. 


TWO  TX-RX  PAIR 


Equipped  with  the  knowledge  of  how  a  single  PBR  system  is  limited  in  its  ability 
to  provide  consistent  detection  coverage,  it  is  important  to  explore  passive  multistatic 
radar  (PMR)  and  how  a  different  transmitter  constellation  affects  the  detection  coverage. 
The  three  parameters  for  transmitter  constellation  selection  are  the  quantity  of 
transmitters  and  their  position  relative  to  the  receiver. 

1.  Bistatic  RCS  with  Two  S-Band  CMRs 

In  this  two-transmitter  PMR  system  setup,  the  two  transmitters  (Txl  and  Tx2)  are 
90°  degrees  apart  and  share  the  same  baseline  distance  of  5000  m.  Although  the 
subsequent  detection  plots  bear  semblance  to  the  detection  plots  of  a  single  Tx  (Figure 
23),  it  is  important  to  note  that  detection  gaps  have  been  reduced  in  all  cases  (Figure  24), 
and  the  two-transmitter  PMR  system  is  now  more  tolerant  of  target  RCS  changes  due  to 
aspect.  From  Figures  25  and  26,  the  detection  plots  show  the  detection  coverage  of  the 
Tx-Tgt-Rx  geometry  has  at  least  one  broadside  reflection  return  from  the  low-RCS  target. 


Figure  23.  SNR  detection  plot  (dB)  with  2Tx  at  S-band  (90°  apart)  having  reflections  off 

the  forward  aspect  of  the  target  ship. 
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Figure  24.  SNR  detection  plot  (dB)  with  2Tx  at  S-band  having  reflections  off  the 

broadside  of  the  target  ship. 
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Figure  25.  SNR  detection  plot  (dB)  at  S-band  with  Txl  with  reflections  off  the  forward 
aspect  of  the  target  ship  and  Tx2  off  the  broadside  of  the  target  ship. 
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Figure  26.  SNR  detection  plot  (dB)  at  S-band  with  Txl  with  reflections  off  the  broadside 
of  the  target  ship  and  Tx2  off  the  forward  aspect  of  the  target  ship. 


2.  Bistatic  RCS  with  X-Band  CMR 

Extending  the  simulation  using  X-band  CMR,  we  show  the  detection  plots  in 
Figures  27  through  30.  They  are  consistent  with  that  using  S-band  CMR  as  the 
opportunistic  transmitter. 


Figure  27.  SNR  detection  plot  (dB)  with  2Tx  at  X-band  (90°  apart)  with  reflections  off 

the  forward  aspect  of  the  target  ship. 
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From  Figures  29  and  30,  it  can  be  observed  that  there  two  distinct  shades  of  blue 
denoting  the  detection  area  about  5000  m  further  away  from  the  receiver.  The  dark 
shades  denote  regions  where  the  SNR  is  above  the  detection  threshold  of  10  dB  but  less 
than  40  dB.  This  dark  shade  of  blue  exists  only  between  ranges  above  5000  m  and  12,000 
m,  while  there  are  repeated  patterns  of  coverage  and  non-coverage  when  the  target  range 
increases. 

To  explain  this  phenomenon,  recall  the  variables  in  Eqs.  (5)  and  (6)  that  affect  the 
SNR.  The  obvious  factor  is  the  combined  value  of  that  is  non-linear  in  nature 

acting  as  the  denominator  in  Eq.  (6).  With  fi  constant,  the  change  in  R^R^  affects  the 
SNR.  The  more  asymmetric  the  Rj  and  are,  the  larger  the  SNR.  This  explains  why  the 

SNR  increases  at  the  ends  of  the  constant  y5-contours;  however,  if  the  bistatic  RCS  returns 
corresponding  to  a  are  small,  there  are  detection  gaps  for  that  value  of  p.  In  summary, 
there  are  two  variables  involved  in  determining  the  SNR,  namely,  the  bistatic  RCS 
returns  which  are  sensitive  to  P  and  RjR^,  where  the  latter  loses  its  dominance  as  the 
target  moves  further  away  from  the  receiver. 


Eigure  28.  SNR  detection  plot  (dB)  with  2Tx  at  X-band  with  reflections  off  the  broadside 

of  the  target  ship. 
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Figure  29.  SNR  detection  plot  (dB)  at  X-band  with  Txl  with  reflections  off  the  forward 
aspect  of  the  target  ship  and  Tx2  off  the  broadside  of  the  target  ship. 
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Figure  30.  SNR  detection  plot  (dB)  at  X-band  with  Txl  with  reflections  off  the  broadside 
of  the  target  ship  and  Tx2  off  the  forward  aspect  of  the  target  ship. 
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3. 


Combination  and  Constellation 


The  approach  taken  to  evaluate  a  transmitter  constellation  is  to  compare  the 
detection  coverage  with  a  monostatic  radar.  That  requirement  is  to  achieve  all-round 
coverage  with  a  minimum  number  of  detection  gaps.  A  series  of  different  constellations 
are  explored  to  establish  a  decision  matrix  to  choose  the  optimal  number  of  opportunistic 
transmitters  to  maximize  coverage  and  minimize  gaps.  To  simplify  the  simulation 
without  sacrificing  validity,  all  subsequent  detection  plots  are  shown  with  strong  RCS 
returns  using  S-band  CMR.  This  is  fair  as  the  low-RCS  target  presents  a  different  aspect 
to  each  transmitter,  and  at  least  one  of  these  combinations  invokes  a  strong  RCS  return 
back  to  the  receiver. 

As  shown  in  Figures  31  through  36,  both  detection  coverage  and  range  can  be 
improved  by  choosing  an  opportunistic  transmitter  that  has  a  larger  baseline  distance 
from  the  receiver.  The  upper  limit  of  the  choice  of  baseline  distance  lies  in  the 
transmitter’s  power  and  the  receiver’s  sensitivity.  Complemented  by  another  transmitter 
90°  or  180°  apart,  the  detection  gap  is  significantly  reduced.  It  is  also  shown  that 
detection  gaps  in  a  particular  bearing  can  be  reduced  by  a  transmitter  that  lies 
orthogonally  to  it. 


Figure  31.  SNR  detection  plot  with  Txl  at  (-5000,0)  and  Tx2  at  (5000,0). 
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Comparing  Figures  31  and  32,  we  see  that  detection  gaps  to  the  right  of  the 
receiver  are  reduced  with  Tx2  at  90°  as  compared  to  180°.  The  converse  is  also  true  for 
detection  gaps  that  are  reduced  forward  of  the  receiver  if  both  Txl  and  Tx2  lie  on  the 
same  baseline  as  the  receiver. 


Figure  32.  SNR  detection  plot  with  Txl  at  (-5000,0)  and  Tx2  at  (0,8000). 


It  is  noteworthy  that  the  detection  gaps  can  be  classified  into  two  types,  interval 
gaps  and  through  gaps.  The  interval  gaps  represent  detection  gaps  that  appear  in  a  lateral 
manner  and  create  loss  in  tracking  as  the  target  moves  towards  or  away  from  the  receiver. 
The  through  gaps  are  detection  gaps  that  allow  the  target  to  evade  detection  and  reach  the 
receiver.  Refer  to  Figure  32  for  a  depiction  of  how  the  interval  gaps  are  detection  gaps 
that  form  concentric  ovals  connecting  the  transmitter  and  receiver.  The  through  gaps  are 
detection  gaps  connecting  with  the  transmitter  (Txl)  on  the  left  or  the  receiver  on  the 
right. 
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Figure  33.  SNR  detection  plot  with  Txl  at  (-5000,0)  and  Tx2  at  (0,10000). 


Figure  34.  SNR  detection  plot  with  Txl  at  (-8000,0)  and  Tx2  at  (0,10000). 
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Figure  35.  SNR  detection  plot  with  Txl  at  ( - 10000,0)  and  Tx2  at  (0,10000). 


Figure  36.  SNR  detection  plot  with  Txl  at  ( - 10000,0)  and  Tx2  at  (10000,0). 
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To  minimize  detection  gaps  as  a  result  of  the  bistatic  angle  formed  between  the 
Tx-Tgt-Rx,  the  optimal  transmitter  constellation  should  have  transmitters  at  every 
quadrant  and  of  varying  baseline  distance.  This  is  consistent  with  the  findings  by  [10]. 

C.  TRANSMITTER  CONSTELLATIONS 

How  different  constellations  formed  using  4Tx-Rx  and  8Tx-Rx  affects  detection 
coverage  is  explored  in  this  section. 

1.  Four  Tx-Rx  Pair 

In  this  simulation,  four  transmitters  are  placed  10000  m  away  from  the  receiver, 
each  occupying  a  quadrant  at  90°  apart  (Figure  37).  It  is  worth  noting  that  100%  detection 
coverage  is  achieved  within  a  20  km  radius  of  the  receiver.  The  second  detection  tier 
between  20  km  to  28  km  has  about  60%  detection  coverage.  An  increase  in  the  baseline 
distance  from  10  km  to  15  km  (Figure  38)  also  saw  an  improvement  in  detection  range  in 
the  second  detection  tier  with  the  reduction  of  through  gaps.  The  improvement  in 
detection  coverage  in  the  second  tier  increased  from  60%  to  approximately  80%. 


Figure  37.  SNR  detection  coverage  plot  with  4Tx  (0°,  90°,  180°  and  270°)  at  10000  m 

from  the  receiver. 
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Figure  38.  SNR  detection  coverage  plot  with  4Tx  (0°,  90°,  180°  and  270°)  at  15000  m 

from  the  receiver. 


Figure  39.  SNR  detection  coverage  plot  with  4Tx  (45°,  135°,  225°  and  315°)  at  15000  m 

from  the  receiver. 
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As  shown  in  Figure  39,  the  4Tx  constellation  placement  with  transmitters  at  45°, 
135°,  225°  and  315°,  respectively,  has  modified  the  second  tier  detection  coverage 
between  20  km  and  28  km.  By  rotating  the  transmitter  constellation  by  45°,  the  location 
of  through  detection  gaps  in  the  second  tier  has  also  rotated.  Essentially,  the  use  of  a 
different  constellation  of  the  four  transmitters  allows  the  receiver  to  direct  optimal 
detection  coverage  in  its  quadrant  of  choice. 

2.  Eight  Tx-Rx  Pair 

The  use  of  eight  transmitters  offers  the  best  all-round  detection  coverage  with 
near  95%  coverage  up  to  a  range  of  30  km.  From  Figures  40  and  41,  one  can  see  the 
detection  coverage  by  eight  transmitters  placed  at  15000  m  and  18000  m  away  from  the 
receiver,  respectively. 


Figure  40.  SNR  detection  coverage  plot  with  8Tx  (0°,  45°,  90°,  135°,  180°,  225°,  270°  and 

315°)  at  15000  m  from  the  receiver. 
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Figure  41.  SNR  detection  coverage  plot  with  8Tx  (0°,  45°,  90°,  135°,  180°,  225°,  270°  and 

315°)  at  18000  m  from  the  receiver. 


In  this  chapter  it  was  shown  how  the  placement  of  the  opportunistic  transmitters 
affects  detection  range  and  detection  coverage.  In  the  next  chapter  the  behavior  of  the 
measurement  errors  and  how  to  minimize  errors  through  the  use  of  opportunistic 
transmitter  placement  is  investigated. 
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IV.  SIMULATION  RESULTS  FOR  MEASUREMENT  ERRORS 


Although  an  eight-transmitter  constellation  offers  detection  coverage  nearly 
identical  to  that  of  a  monostatic  radar,  there  are  measurement  errors  that  have  to  factor 
into  consideration.  A  target  detection  with  an  unacceptably  large  error  is  not  useful. 
Depending  on  the  geometry  of  the  PBR,  we  found  that  the  worst-case  measurement  error 
varies  across  range  and  bistatic  angle.  As  shown  in  Eqs.  (11)  through  (13),  the  RCS  of  the 
target  does  not  have  a  role  to  play  in  measurement  error;  measurement  error  calculations 
come  after  target  detection  with  a  SNR  of  more  than  10  dB.  Despite  the  varying  nature  of 
the  measurement  errors,  it  is  possible  to  use  a  careful  choice  of  transmitter  constellation 
to  minimize  the  errors. 

A.  UNCERTAINTY  ELLIPSE 

As  explained  earlier  with  the  aid  of  Figure  9,  every  target  position  with  a  Tx-Rx 
geometry  pair  has  its  own  measurement  error.  The  measurement  errors  typically  span 
three  different  axes  and  form  an  uncertainty  ellipse  around  the  target  position.  It  should 
be  interpreted  as  the  target’s  actual  location  is  anywhere  within  the  ellipse.  If  the  target  is 
also  tracked  by  another  Tx-Rx  pair,  the  resultant  uncertainty  ellipse  is  likely  to  have  the 
dominant  measurement  error  on  a  different  axis.  In  such  cases,  the  uncertainty  of  the 
target’s  position  is  further  reduced  to  the  overlapped  area.  This  improvement  extends  to 
multiple  pairs  of  Tx-Rx  on  the  same  target. 

To  be  commensurate  with  the  detection  plots  shown  in  Chapter  III,  the 
uncertainty  ellipse  of  a  different  target  location  is  plotted  using  the  same  transmitter 
constellations.  The  MATLAB  source  code  for  the  following  plots  is  included  as 
Appendix  G. 

I.  Two  Tx-Rx  Pairs 

In  this  case,  the  target  is  20.6  km  away  from  the  receiver  with  two  transmitters  90° 
apart,  as  shown  in  Figure  42.  The  inputs  to  this  simulation  plot  assume  transmitter 
positioning  update  by  the  MAIS  at  5  s  and  a  net  displacement  of  the  target  at  10  knots. 
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The  uncertainty  ellipse  around  the  target  (-5000,  20000)  is  the  overlap  area  between  the 
two  uncertainty  ellipses  formed.  At  this  extended  distance,  the  target’s  measurement  error 
area  is  approximately  120  m  all  round.  Comparing  with  the  target  length  of  114  m,  we 
assess  the  error  as  reasonable. 


Figure  42.  Uncertainty  ellipse  of  Target  (-5000,20000)  with  2Tx  at  90°  apart. 


Figure  43.  Uncertainty  ellipse  of  Target  (5000,20000)  with  2Tx  at  90°  apart. 


A  change  of  target  position  to  the  right  of  the  Tx-Tx  constellation  shows  a  slightly 
larger  uncertainty  ellipse  due  to  geometry  (Figure  43).  As  the  target  approaches  the  Tx- 
Rx  geometry,  one  can  see  the  uncertainty  ellipse  reduces  with  an  all-round  error  of  40  m 


54 


as  shown  in  Figure  42.  In  certain  geometries  (Figure  45),  the  uncertainty  ellipse  is  the 
same  from  two  different  transmitter  locations. 


Figure  44.  Uncertainty  ellipse  of  Target  (5000,5000)  with  2Tx  at  90°  apart. 


Figure  45.  Uncertainty  ellipse  of  Target  ( -2000,10000)  with  2Tx  at  180°  apart. 


2.  Four  Tx-Rx  Pairs 

In  a  four-transmitter  constellation,  the  effect  on  the  uncertainty  ellipse  is  similar 

(Figure  46).  The  major  factors  affecting  the  size  of  the  uncertainty  ellipse  are  the 

beamwidth  of  the  transmitter  and  DF  accuracy  of  the  receiver.  As  the  target  approaches 

the  transmitter  or  the  receiver,  the  dominant  factor  changes.  The  important  observation 

here  is  that  the  movement  of  the  transmitter,  receiver  and  target  has  very  little  influence 
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in  determining  the  uncertainty  ellipse.  With  the  target  at  (4500,4500)  and  using  the  S- 
band  CMR  with  a  beamwidth  of  1.8°,  one  can  see  from  Figure  47  that  the  shape  of  the 
uncertainty  ellipse  has  changed.  Extending  this  further,  we  see  the  uncertainty  ellipse  of 
the  same  constellation  with  the  transmitter  changed  to  X-band  CMR  with  a  beamwidth  of 
0.95°  as  illustrated  in  Figure  48.  For  Figure  49,  the  transmitter,  receiver  and  target  speeds 
have  been  reduced  to  zero  from  10  knots.  It  can  be  seen  that  the  data  age  and  target  speed 
have  negligible  effects  on  the  measurement  error  in  this  setup. 
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Figure  46.  Uncertainty  ellipse  of  Target  (8000,13000)  with  4Tx  at  90°  apart  from  each 

other. 
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Figure  47.  Uncertainty  ellipse  of  Target  (4500,4500)  with  4Tx  at  90°  apart  from  each 

other  using  S-band  CMR. 
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Figure  48.  Uncertainty  ellipse  of  Target  (4500,4500)  with  4Tx  at  90°  apart  from  each 

other  using  X-band  CMR. 


Figure  49.  Uncertainty  ellipse  of  Target  (4500,4500)  with  4Tx  at  90°  apart  from  each 
other  using  X-band  CMR  and  zero  net  speed  for  transmitter,  receiver  and  target. 

3.  Eight  Tx-Rx  Pairs 

One  can  expect  similar  results  from  an  eight-transmitter  constellation  when 
compared  to  the  four- transmitter  constellation;  however,  it  is  worthwhile  to  note  that  the 
higher  the  number  of  transmitters  and  the  more  they  are  spatially  diversified,  the  higher 
the  number  of  target  positions  that  the  PMR  can  achieve  and  the  smaller  the  uncertainty 
ellipse  that  results.  Therefore,  a  higher  accuracy  of  detection  can  be  obtained.  In  Figure 
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50,  an  eight-transmitter  constellation  is  set  up  similar  to  that  shown  in  Chapter  III.  Letting 
the  target  position  be  (4500,4500),  it  can  be  seen  that  a  very  small  measurement  error  of 
approximately  20  m  is  achievable  with  only  four  transmitters  out  of  the  eight  transmitters 
used.  This  should  be  interpreted  as  the  PMR  system  now  has  eight  data  sets  to  compare 
and  choose  from,  so  the  best  sets  can  be  used  to  represent  the  uncertainty  ellipse.  To 
achieve  the  same  low  measurement  error,  one  needs  to  actively  seek  transmitter  spatial 
diversity  if  one  wishes  to  focus  on  obtaining  high  accuracy  on  one  target. 
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Figure  50.  Uncertainty  ellipse  of  Target  (4500,4500)  with  8Tx  at  45°  apart  from  each 
other  using  S-band  CMR  while  the  target  is  moving  at  the  speed  of  10  knots. 
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V.  CONCLUSION 


A.  SUMMARY  OF  THESIS  FINDINGS 

Using  opportunistic  transmitters  in  the  form  of  civil  marine  radar  (S-band  and  X- 
band)  to  detect  maritime  targets  passively  is  both  viable  and  has  comparable  performance 
to  monostatic  radars.  However,  in  the  case  of  detecting  low-RCS  targets,  the  PMR 
outperforms  the  former.  The  highly  regulated  requirements  instilled  by  the  IMO  have 
ensured  that  most  of  maritime  vessels  out-at-sea  have  the  necessary  information  to 
support  a  PER  and  PMR.  With  the  world’s  trade  routed  via  the  sea,  there  is  no  shortage 
of  suitable  vessels  possessing  at  least  one  CMR  and  one  MAIS  on  the  main  shipping 
routes.  Among  all  the  suitable  vessels,  larger  ships  like  oil  tankers,  container  ships,  etc., 
have  their  CMR  installed  at  the  highest  point,  which  extends  the  radar  horizon  and 
detection  range.  The  narrow  beamwidth  and  high  transmit  power  of  the  CMR 
conveniently  lend  themselves  to  exploitation  by  a  PER  system  in  terms  of  both  detection 
accuracy  and  detection  range. 

The  unique  geometry  formed  among  the  opportunistic  transmitter,  passive 
receiver  and  unsuspecting  target  provides  exceptional  military  applications  not  achieved 
by  current  military  technologies.  It  offers  true  passive  detection  (with  zero  transmission 
from  the  receiving  warship)  while  at  the  same  time  supports  detection  coverage 
comparable  to  a  monostatic  radar.  Passive  detection  means  there  is  no  minimum  range 
detection  that  is  unavoidable  in  a  monostatic  radar.  In  the  same  vein,  there  is  no  need  for 
‘blanking’  of  the  ESM  receiver  every  time  the  onboard  radar  transmits.  It  allows  for 
continuous  monitoring  of  the  EW  spectrum  for  intelligence  and  threat.  On  the 
vulnerability  front,  this  has  removed  the  risk  posed  by  ARM  and,  at  the  same  time,  has 
enhanced  mission  survivability,  allowing  the  warship  to  continue  with  tracking 
information  without  a  working  radar  or  with  a  damaged  radar. 

This  same  unique  geometry  also  improves  the  warship’s  probability  of  detecting  a 
low-RCS  ship  treated  both  in  shaping  and  RAM  as  compared  to  a  monostatic  radar.  The 
use  of  low  frequency  S-band  CMR  poses  serious  limitations  to  the  use  of  RAM  treatment 
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on  a  low-RCS  ship.  It,  therefore,  restricts  the  effectiveness  of  low-RCS  ships  in  the  use  of 
shaping.  It  was  shown  in  the  bistatic  RCS  measurement  of  the  model  based  on  the  RSN 
Formidable-class  stealthy  frigate  that  shaping  alone  is  capable  of  reducing  the  median 
RCS  to  values  not  higher  than  200  m  .  With  all  else  being  equal  with  the  effect  of  RCS 
on  detection,  the  PER  system  has  an  added  advantage.  The  bistatic  RCS  simulation  using 
FEKO  shows  that  more  incidence  angles  exist  that  demonstrate  a  spike  in  RCS  returns. 
Simply  put,  there  are  more  degrees-of-freedom  created  by  the  transmitter-target-receiver 
geometry  that  has  a  higher  probability  of  directing  spikes  in  the  RCS  returns  toward  the 
receiver.  In  the  case  of  a  monostatic  radar,  there  are  limited  ways  to  maximize  RCS 
returns  because  the  shaping  of  a  low-RCS  is  designed  to  reduce  returns  in  the  direction  of 
the  incidence  angle. 

The  simplicity  of  PMR  also  lies  in  the  fact  that  there  are  no  expensive  systems  to 
build.  Instead,  all  the  information  that  the  PMR  system  needs  is  already  available  on  most 
warships.  This  information  can  be  extracted  and  then  used  to  carry  out  tracking. 

1.  Conclusion  for  Passive  Detection  Coverage 

The  passive  detection  range  and  coverage  are  related  to  many  parameters  on  the 
transmitter,  target  and  receiver.  On  the  opportunistic  transmitter  end,  the  transmit  power 
and  antenna  gain  are  the  two  parameters  that  determine  the  maximum  detection  range. 
However,  it  is  the  relative  position  of  the  transmitter  in  the  PER  geometry  that 
determines  the  detection  coverage  or  detection  gaps.  On  the  receiver  end,  the  deciding 
parameters  are  the  sensitivity  and  bandwidth  of  the  ESM  receiver.  These  two  parameters 
affect  the  SNR  and,  thereby,  change  maximum  detection  range.  Similar  to  the 
opportunistic  transmitter,  it  is  the  geometry  of  the  PER  that  affects  the  detection  gaps. 
Collectively,  the  factor  that  is  used  here  to  quantify  the  geometry  is  the  bistatic  angle  P, 
defined  earlier  in  this  thesis.  Incidentally,  the  bistatic  angle  determines  the  bistatic  RCS 
of  what  is  presented  to  the  receiver  depending  on  the  incidence  angle  of  the  transmitter’s 
radar  emission. 

Prom  the  simulation  plots  presented  in  Chapter  III,  it  was  shown  that  when  more 
transmitters  are  spaced  further  apart,  they  provide  better  overall  coverage  with  small 
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detection  gaps.  The  minimum  number  of  transmitters  is  four,  placed  90  degrees  apart, 
with  one  in  each  quadrant.  To  achieve  maximum  coverage  with  this  particular 
combination  of  system  parameters  and  transmitter  constellation,  a  baseline  distance  of 
18000  m  is  recommended.  By  rotating  the  quadrant  placement  of  transmitters,  one  can 
direct  the  detection  gaps  away  from  a  specific  bearing  of  choice  to  improve  detection. 

Using  more  than  four  transmitters  in  the  suggested  constellation  does  not 
significantly  extend  the  detection  range.  As  shown  in  the  eight-transmitter  simulation,  it 
does  improve  detection  coverage  from  approximately  80%  to  95%  through  the  reduction 
of  detection  gaps.  This  translates  to  the  ability  of  the  PMR  to  track  a  target  continuously 
provided  the  target  offers  a  strong  bistatic  return  to  all  eight  transmitters.  In  all  likelihood, 
a  low-RCS  target  would  have  presented  at  least  two  to  four  strong  RCS  returns  to  either 
the  quadrant-based  constellation  using  four  transmitters  or  the  two-layered  quadrant- 
based  constellation  using  eight  transmitters.  If  the  availability  of  such  a  constellation  is 
not  possible  out  at  sea,  one  should  position  oneself  to  create  such  a  constellation  without 
jeopardizing  a  passive  posture. 

2.  Conclusion  for  Detection  Accuracy 

There  are  two  dominant  factors  that  affect  detection  accuracy,  namely  the 
beamwidth  of  the  transmitter  and  the  DF  accuracy  of  the  receiver.  It  was  shown  in 
multiple  figures  that  the  worst-case  measurement  error  comes  up  to  approximately  ±120 
m  for  a  low-RCS  114  m  frigate  at  a  range  of  20  km.  In  most  cases,  a  lower  measurement 
can  be  achieved  using  multiple  transmitters  by  overlapping  their  individual  uncertainty 
ellipses.  It  was  found  that  the  four- transmitter  and  eight-transmitter  constellations 
recommended  for  maximum  detection  coverage  are  also  optimized  to  reduce  the  target’s 
uncertainty  ellipse.  It  was  shown  that  with  an  eight-transmitter  constellation,  the  area 
overlapped  by  all  the  uncertainty  ellipses  have  given  rise  to  a  measurement  error  of  20  m 
for  the  same  low-RCS  target  at  a  range  of  6400  m. 
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B.  FUTURE  WORK 


Improvements  are  needed  in  the  simulation  model.  Higher  resolution  RCS  data 
should  be  provided.  This  is  possible  by  calling  FEKO  directly  from  the  Matlab  simulation 
with  the  exact  incidence  and  observation  angles.  Also,  the  effect  of  multipath  from  the 
sea  surface  should  be  included. 

A  natural  extension  of  this  research  is  to  look  into  tracking  algorithms  for  such  a 
PMR  given  that  there  is  only  range  and  no  Doppler  information.  There  are  a  couple  of 
proposed  tracking  methods  that  might  be  suitable  for  the  PMR  system.  There  are 
possibilities  for  extracting  target  Doppler  information  from  opportunistic  targets  using  the 
coherent-on-receive  method  [21].  Such  receivers  support  a  pseudo-coherent  reception  of 
radar  signals  to  extract  Doppler  information.  With  both  range  and  Doppler  information  in 
a  pseudo-coherent  form,  the  improved  PMR  system  can  generate  ambiguity  function  for 
better  detection  accuracy  and  also  adopt  typical  tracking  methods  used  in  monostatic 
radars. 

With  detection  and  tracking  assured,  it  is  possible  to  take  a  step  further  into 
classification  and  identification  of  targets.  From  [6],  forward  scatter  and  backscatter 
based  ISAR  are  currently  being  investigated  for  their  ability  to  detect  target  outline.  It 
might  be  interesting  to  consider  adding  another  dimension  to  this  by  using  bistatic  ISAR. 
In  combination,  the  three  ISAR  plots  might  offer  added  information  about  the  target  to 
help  in  classification  and  identification. 


62 


APPENDIX  A.  SPECIFICATION  SHEET  OF  MANTADIGITAL 
RADAR  BY  KELVIN  HUGHES 
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APPENDIX  B.  SPECIFICATION  SHEET  OF  FAR28X7  SERIES  BY 
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440  VAC  1  o.  SOSO  Hz  wY)  ooUO(^4  tranafprmor  RI>1801 
X  Antenna  UnH 
rAIW837S/2637SW 

200  VAC  3J)  A  3p  50  Hz:  220  VAC  3.0  A  (XS  A  tor  HSC*). 
30.60  Hr.  380  VAC  IS  A  3p  50  Hz;  440  VAC.  IS  A  (17  A  tor 
HSC*).  3p  60  Hz 
•  V  f-AH-2637S  Ofly 

110  VAC.  3p  60  Hz  with  RLLS693.  220  VAC.  3#  50  Hz  Wtth 
nU-6622. 440  VAC.  3o.  50  Hz  w«i  nLi-5466-1 


EQUIPMENT  LIST 
Standard 

X  OaplayUnitMlWSl 
X  ProeoaaorUnnRPU-013 
X  Fufr^opyboord  Control  Und  nCLM)14 
Trackbal  Control  Unit  RClMIS 
(Spoofy  when  ordamg) 

A  Anaanoa  Unit  with  eoG|a(lS3D'4Q<B0m) 

X  Power  Supply  urn  PSLK)07  tor  hAfWilTS 
X  Standard  Spa*»  Parta  and  Inttalanon  Maiarlala 
Option 

1  Perlormonco  Monitor  P1A31  for  XdraxKS 
PtA41  for  SdMrd 
?«  Homola  Corcml  Unit  RC0016 
X  Gyre  Iraartaeo  QOIO 

A  DVMnalogROBCP'TvoraldiK)tCX*03-ilO(SXQAoutpat) 
X  RGB  Connoctp'  OSUB-BNC-1  (lor  VDR) 

X  Memory  Card  Inipimce  Un8  CU-200 
7.  rr»iatormerRL>180X'546X1'b69l*6d22 
X  Rocti«wrRI>342A17468 
X  Juncten  BcK  RJBdXn 
ia  Antenna  Cabla  RW^ggoo 

11.  Hand  Qrp  FP03-09640 

12.  Bradur  rP03-09620 
IX  Swilchng  Hub  HUB-100 
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APPENDIX  C.  SPECIFICATION  SHEET  OF  QR026  EW  RECEIVER 

BY  TELEYDNE  DEFENSE 


TELEDYNE 

DEFENCE 

A  Teledyne  Technologies  Company 


FeatureG 


Applications 


Product  Description 
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TELEDYNE 

DEFENCE 

A  Teiedyne  Technologies  Company 


0.5  to  IBGHz 

High  Performance  Receiver 

TBchnical  Dou  QR026 


Beccrical  Specifcacion 


•HP  Ir-p-ji  -reausncy  fPJTye: 

CpLlan3(  exLensbns  za  :■  40Gi-lz 

•Maximum  ..-.uju  CW 

•eiaan 

•LD  Re  naiiat.;^T 

-SOdBnkMi 

•Ga^  20CB  nom  [3CGMM2  IFl 

POdB  nom  .irnv-'ir) 

•Gan  war-auPrn  :  BWl 

stana  ir  uhe  IP  aanff^^iotn 

1-^aaB  max  nr  any  antiMquny  na  any  if 

•Gan  vsnau  rn  :ar;v  ICO’.'H/  c'  if  bwi 

^AI.SiB  wm. 

•Gpgup  ne(av  Varooon  [SOOM-;  IF  BWl 

Ita  ow  itH  cnm  OOk  Bw:  IOm  ow  Ub  m 

Qefe-jvtn  3np  uyt)  -ilns  in  ccmirai  ICOK 

Itt 

•inpuL  AtcETuatcr 

BadBMva  acapi 

•IF  AnuenuBtor 

ISBI*!.  liBac^iatB 

•VS'AW  tRf  out] 

U:1 

•Nase  Figure 

lasMCiimria-aB 

•irpjt  P1dB 

•Input  IP3 

-zdanfmn^Qdm 

•Single  Itne  SFOfl 

BOdB  im  tSOOMM  aiM 

•in-tano  ?  inne  SFon 

SQOB  mn  aOOMMK  BM 

•inte(T.ally  Generated  Spirtous 

-TOBn  nMB  BiuiL  apuc  pomt 

•Tune  Speed 

1 

1 

1 

1 

i 

1 

G 

i 

fi 

1 

•Tune  ftcctifacy 

auz 

•lure  PTesoijticr 

IMHzMtMromM 

•Pnase  Naiss 

in-/ 

lOkN/ 

1MH/ 

IDNVt/ 

-USdBQHE 

•SI/P 

19  n  m  X  zu  Hgn  Mnna  obannfll 

•ItiLal  pnn»ep  Corrsunptcn 

<BOW  ai  Mi  rack  oonmunatHi 

•Input  PXiwer  Supp^ 

94  -  aa  WC.  47  - 93  Hz 

•'i^^nr 

<  lao 

•Opera' inp  Te.Tperaujre  FTsr.ga 

-in:  CD  ran: 

•□P'_-j'-a 

•  AM  IMO  RQMr  Odqin 

DHiAMnie  ftapurgradan  n  cuanxiierm 

68 


APPENDIX  D.  SPECIFICATION  SHEET  FOR  DF-A0062  DF 
RECEIVER  BY  POYNTING  DEFENSE 


VERSION; 


Poynting 

Antennas 


2.0 


Direction  Finding  and 
Monitoring  Antenna 

20  -  6000  MHz _ 

Product  Code:  DF-A0062 


PRODUCT  DESCRIPTION: 

This  arecaon  llnang  antenna  covere  a  itequency  range 
or  20  MHz  to  6  GHz.  A  separate,  out  nay  ntegrateo 
active  monitoring  arcefYta  gtves  higli  senslOvity  onwi- 
drectlonal  coverage  on  me  same  aas  as  me  3F 
aniema.  eimout  irterTePng  enn  me  Df  array. 

The  tia-slze  etements  on  all  oanos  gtve  exceient  3F 
sensmvtty.  umate  angular  resouoon  tir  snng  signals 
IS  i»ell  under  l*  tor  most  or  me  irequency  range  Gipoie 
eiemerts  provtoe  good  ctoss-pciansation  retecoon.  and 
*air  pertornarKe  tir  signals  ainving  from  up  to  15* 
aoove  or  oeiow  the  honzon. 

The  ircegrated  rmnitonng  antenna  is  m  Mo  hands, 
mounted  aoove  the  Band  C  &  D  DF  antenna  rxxMe 
housed  m  the  radome  Eacn  hand  or  me  monitonng 
anay  is  ampinea  at  the  top  or  me  caUe.  and  passively 
oomewea  to  ghe  oontlnuous  coverage  over  the 
*requercy  range  20  seiz  to  6  GHz.  «im  a  sHgK  output 
connector 


SPECIFICATTONS: 


Product  cooes 

DF-AO062 

DF  sfln  montDTtrtq  up  to  6  GHZ 

E)«cmcai-OF 

DF  tequency  range 

BandA2G-3XMHz: 

Band  5:  300 -1X0  MHz; 

Band  C:  lOX  -  30X  MHz; 
BandD:3OX-60X  MHz 

Normal  iroir.  moedanoe 

SOD 

Artennatype 

5-eiement  DF  irteneiometer 
iFrom  5  to  2-channel  recever 
companoiei 

PoansaPon 

vertical 

Oumut  caCNes 

RG40C  cauee 

DF  oomecmre 

20xTNCmaie 

Monitonng  connector 

t  xN  male 

Eiectncsi  -  momtonna 

Frequency  ranqe 

20-60X  MHZ 

Nomnal  H-olane  oeamwldm 

36C* 

Momnal  E-piane  heanmndm 

60- 

Typical  vswR- 

2:1 

PoansaPon 

vertical 

Connector  type 

N-tvpemaie 

Nomnai  inoedance 

SOD 

inout  ^oitaoe  r^a  coaxr 

13-24VDC 

Input  nsrenr 

« iSOmA 

Poaier  corsumpaon  rnomt' 

«2.2SW 

optdB  ihp.r 

>  ti  dBm 

OIP2  itvo  r 

>  31  dBm 

OIP3  itvo  • 

>21  dBm 

Gensitvty 
W) 
larN- 
Gd6.  BW- 
1  HZ)’ 

20 -too  MHZ 

-20  dByiViW 

IX  -10X 
tJHZ 

-30dBtiVrm 

1-3GHZ 

-35  dBuvm 

3-6GMZ 

-40  dBiivm 

iUftiinirH' 

CrosseecSonal  sand  area 

1  50  m*  rmctupnq  antema  Ewnchi 

Maxrman  lamd  soeed 

t&C  invh  iieimout  ice  loadi 

AssemiMed  neian 

3.60  m 

AssemUeo  Oameter  imaxi 

Z70m 

SHaamo  dmenoons 

220  m  X  0.90  m  X  0.60  m 

weigM  or  antenna 

Inauang  sNppInq  container 

Xkg 

135  iq 

‘n  active  mode 


This  DF  antenna  is  designed  to  he  usatue  vMh  either  a 
5-  or  2-chamei  phase-sensiove  recetier  and  coireiattve 
aigoitmm.  cnaractensatlcn  or  the  antema  car  he 
pained  on  request. 

Related  products:  DF-ADOM.  DF-AD057-03.  DF-AO038 
and  DF-A0098 


ELECTRJCAi.  FEATURES; 

•  Full^izeDF 

•  Wdehand  DF 

•  5-eiemert  imerterometer 

•  Opcmised  ihr  2-cnannei  receivers 

•  High  sensitivity  omrs  antenna  ntegrated 

•  Tactical  rapid  depioyiStoar  design  atm  integrated  dust  caps 
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APPENDIX  E.  SAMPLE  OF  BISTATIC  RCS  RESULTS  OF  MODEL 
OF  RSN’S  FORMIDABLE-CLASS  FRIGATE 


The  following  table  has  a  sample  of  the  bistatic  RCS  values  in  dBsm  for  the 


frequency  of  3050  MHz  and  the  model  shown  in  Figures  6  and  7  as  generated  by  FEKO. 


Reflected 
Angle  in  ° 

Incidence  Angle  in  ° 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

31.53 

31.53 

21.74 

21.74 

14.80 

14.80 

15.80 

15.80 

14.62 

14.62 

2 

28.31 

28.31 

23.84 

23.84 

14.42 

14.42 

12.99 

12.99 

16.08 

16.08 

3 

21.23 

21.23 

18.88 

18.88 

13.18 

13.18 

15.45 

15.45 

5.49 

5.49 

4 

14.15 

14.15 

14.28 

14.28 

18.14 

18.14 

8.34 

8.34 

12.75 

12.75 

5 

19.32 

19.32 

13.93 

13.93 

13.90 

13.90 

17.18 

17.18 

15.36 

15.36 

6 

12.99 

12.99 

13.15 

13.15 

8.37 

8.37 

12.13 

12.13 

13.46 

13.46 

7 

13.99 

13.99 

14.61 

14.61 

16.79 

16.79 

9.87 

9.87 

14.51 

14.51 

8 

1.22 

1.22 

16.13 

16.13 

12.77 

12.77 

13.42 

13.42 

15.77 

15.77 

9 

11.07 

11.07 

16.83 

16.83 

13.94 

13.94 

14.12 

14.12 

17.87 

17.87 

10 

15.99 

15.99 

12.30 

12.30 

16.34 

16.34 

13.45 

13.45 

15.53 

15.53 

11 

17.23 

17.23 

7.59 

7.59 

13.60 

13.60 

15.86 

15.86 

14.21 

14.21 

12 

12.32 

12.32 

17.34 

17.34 

16.82 

16.82 

16.59 

16.59 

16.57 

16.57 

13 

10.85 

10.85 

16.89 

16.89 

18.07 

18.07 

13.87 

13.87 

16.67 

16.67 

14 

16.77 

16.77 

16.29 

16.29 

13.15 

13.15 

17.67 

17.67 

14.66 

14.66 

15 

11.04 

11.04 

15.80 

15.80 

15.05 

15.05 

13.57 

13.57 

17.75 

17.75 

16 

14.73 

14.73 

17.79 

17.79 

17.10 

17.10 

13.81 

13.81 

17.36 

17.36 

17 

15.81 

15.81 

13.67 

13.67 

13.97 

13.97 

19.43 

19.43 

15.44 

15.44 

18 

16.89 

16.89 

14.36 

14.36 

11.99 

11.99 

17.78 

17.78 

14.46 

14.46 

19 

15.63 

15.63 

12.89 

12.89 

16.94 

16.94 

15.27 

15.27 

14.25 

14.25 

20 

14.52 

14.52 

17.17 

17.17 

16.02 

16.02 

14.63 

14.63 

13.51 

13.51 

21 

15.27 

15.27 

18.76 

18.76 

13.62 

13.62 

15.95 

15.95 

15.05 

15.05 

22 

16.27 

16.27 

18.62 

18.62 

16.30 

16.30 

15.81 

15.81 

16.12 

16.12 

23 

18.98 

18.98 

15.81 

15.81 

16.94 

16.94 

14.33 

14.33 

16.22 

16.22 

24 

18.12 

18.12 

14.14 

14.14 

15.03 

15.03 

16.72 

16.72 

16.61 

16.61 

25 

11.94 

11.94 

16.35 

16.35 

16.15 

16.15 

16.83 

16.83 

17.41 

17.41 

26 

13.82 

13.82 

13.11 

13.11 

15.65 

15.65 

16.98 

16.98 

14.61 

14.61 

27 

12.48 

12.48 

17.23 

17.23 

16.91 

16.91 

15.62 

15.62 

15.69 

15.69 

28 

14.67 

14.67 

17.02 

17.02 

17.39 

17.39 

15.53 

15.53 

16.67 

16.67 

29 

15.20 

15.20 

14.33 

14.33 

15.00 

15.00 

16.40 

16.40 

14.78 

14.78 

30 

14.46 

14.46 

13.59 

13.59 

14.24 

14.24 

16.29 

16.29 

13.05 

13.05 

31 

14.06 

14.06 

14.35 

14.35 

14.48 

14.48 

15.07 

15.07 

12.30 

12.30 

32 

14.78 

14.78 

15.32 

15.32 

12.75 

12.75 

14.47 

14.47 

12.93 

12.93 

33 

15.81 

15.81 

13.42 

13.42 

12.58 

12.58 

15.02 

15.02 

13.34 

13.34 

34 

14.25 

14.25 

13.67 

13.67 

15.18 

15.18 

14.85 

14.85 

14.72 

14.72 

35 

11.97 

11.97 

16.41 

16.41 

15.84 

15.84 

13.54 

13.54 

15.30 

15.30 

36 

12.99 

12.99 

14.34 

14.34 

12.07 

12.07 

12.87 

12.87 

13.00 

13.00 

37 

15.63 

15.63 

13.65 

13.65 

13.51 

13.51 

15.26 

15.26 

14.73 

14.73 

38 

10.62 

10.62 

13.39 

13.39 

14.44 

14.44 

11.10 

11.10 

13.19 

13.19 

39 

14.63 

14.63 

13.80 

13.80 

11.52 

11.52 

14.83 

14.83 

13.27 

13.27 

40 

31.53 

31.53 

13.27 

13.27 

12.87 

12.87 

12.94 

12.94 

12.26 

12.26 
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Reflected 
Angle  in  ° 

Incidence  Angle  in  ° 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

141 

25.17 

25.17 

19.12 

19.12 

13.23 

13.23 

14.39 

14.39 

13.48 

13.48 

142 

19.37 

19.37 

15.16 

15.16 

18.71 

18.71 

17.20 

17.20 

15.28 

15.28 

143 

19.92 

19.92 

20.13 

20.13 

16.87 

16.87 

15.18 

15.18 

8.47 

8.47 

144 

9.21 

9.21 

19.72 

19.72 

0.53 

0.53 

17.46 

17.46 

15.95 

15.95 

145 

16.35 

16.35 

17.15 

17.15 

15.94 

15.94 

10.64 

10.64 

16.38 

16.38 

146 

14.20 

14.20 

10.74 

10.74 

12.02 

12.02 

17.08 

17.08 

2.94 

2.94 

147 

16.75 

16.75 

4.29 

4.29 

16.85 

16.85 

13.11 

13.11 

20.67 

20.67 

148 

15.82 

15.82 

14.01 

14.01 

13.23 

13.23 

11.49 

11.49 

15.14 

15.14 

149 

18.16 

18.16 

14.15 

14.15 

12.81 

12.81 

7.26 

7.26 

9.03 

9.03 

150 

11.65 

11.65 

17.56 

17.56 

12.06 

12.06 

13.12 

13.12 

13.29 

13.29 

151 

18.53 

18.53 

20.68 

20.68 

12.32 

12.32 

13.00 

13.00 

15.54 

15.54 

152 

18.15 

18.15 

15.65 

15.65 

16.95 

16.95 

14.03 

14.03 

7.85 

7.85 

153 

18.71 

18.71 

9.39 

9.39 

13.31 

13.31 

16.05 

16.05 

16.40 

16.40 

154 

15.74 

15.74 

12.92 

12.92 

16.24 

16.24 

13.13 

13.13 

12.76 

12.76 

155 

19.69 

19.69 

9.29 

9.29 

15.34 

15.34 

17.32 

17.32 

15.99 

15.99 

156 

5.95 

5.95 

18.21 

18.21 

20.14 

20.14 

13.00 

13.00 

13.18 

13.18 

157 

16.99 

16.99 

16.96 

16.96 

21.29 

21.29 

10.72 

10.72 

18.52 

18.52 

158 

15.73 

15.73 

16.43 

16.43 

18.04 

18.04 

20.36 

20.36 

22.92 

22.92 

159 

18.17 

18.17 

18.27 

18.27 

15.87 

15.87 

23.21 

23.21 

19.38 

19.38 

160 

14.31 

14.31 

22.57 

22.57 

21.28 

21.28 

22.04 

22.04 

25.05 

25.05 

161 

9.92 

9.92 

16.67 

16.67 

16.61 

16.61 

23.11 

23.11 

15.39 

15.39 

162 

20.86 

20.86 

17.48 

17.48 

17.35 

17.35 

16.15 

16.15 

28.07 

28.07 

163 

18.82 

18.82 

22.26 

22.26 

10.20 

10.20 

16.73 

16.73 

29.82 

29.82 

164 

16.47 

16.47 

21.29 

21.29 

20.95 

20.95 

11.73 

11.73 

25.46 

25.46 

165 

9.39 

9.39 

7.17 

7.17 

23.36 

23.36 

10.15 

10.15 

28.05 

28.05 

166 

22.88 

22.88 

23.93 

23.93 

22.85 

22.85 

18.72 

18.72 

31.73 

31.73 

167 

22.27 

22.27 

14.01 

14.01 

20.75 

20.75 

25.96 

25.96 

31.84 

31.84 

168 

16.65 

16.65 

23.70 

23.70 

16.55 

16.55 

11.50 

11.50 

35.03 

35.03 

169 

29.95 

29.95 

20.54 

20.54 

14.62 

14.62 

26.24 

26.24 

37.27 

37.27 

170 

30.67 

30.67 

22.16 

22.16 

22.02 

22.02 

26.95 

26.95 

34.22 

34.22 

171 

24.17 

24.17 

23.35 

23.35 

19.97 

19.97 

27.68 

27.68 

36.87 

36.87 

172 

26.07 

26.07 

24.93 

24.93 

29.18 

29.18 

34.31 

34.31 

49.15 

49.15 

173 

32.45 

32.45 

27.50 

27.50 

35.31 

35.31 

34.74 

34.74 

39.48 

39.48 

174 

34.28 

34.28 

26.45 

26.45 

35.58 

35.58 

48.46 

48.46 

36.34 

36.34 

175 

28.64 

28.64 

27.86 

27.86 

39.40 

39.40 

36.93 

36.93 

39.02 

39.02 

176 

38.80 

38.80 

37.34 

37.34 

44.49 

44.49 

33.15 

33.15 

30.89 

30.89 

177 

37.77 

37.77 

39.45 

39.45 

43.99 

43.99 

24.56 

24.56 

29.71 

29.71 

178 

45.24 

45.24 

47.41 

47.41 

41.08 

41.08 

33.65 

33.65 

33.35 

33.35 

179 

61.50 

61.50 

50.60 

50.60 

43.77 

43.77 

31.57 

31.57 

34.65 

34.65 

180 

78.86 

78.86 

49.50 

49.50 

43.67 

43.67 

37.55 

37.55 

27.72 

27.72 

181 

61.64 

61.64 

61.46 

61.46 

47.32 

47.32 

41.96 

41.96 

36.18 

36.18 

182 

46.01 

46.01 

79.67 

79.67 

48.33 

48.33 

41.59 

41.59 

31.70 

31.70 

183 

37.48 

37.48 

57.93 

57.93 

54.60 

54.60 

32.26 

32.26 

39.70 

39.70 

184 

37.65 

37.65 

48.32 

48.32 

81.19 

81.19 

50.41 

50.41 

44.34 

44.34 

185 

34.87 

34.87 

40.98 

40.98 

49.67 

49.67 

49.87 

49.87 

44.81 

44.81 

186 

33.43 

33.43 

42.10 

42.10 

50.11 

50.11 

82.18 

82.18 

36.68 

36.68 

187 

29.38 

29.38 

34.66 

34.66 

43.76 

43.76 

49.93 

49.93 

51.23 

51.23 

188 

27.37 

27.37 

26.55 

26.55 

41.99 

41.99 

41.30 

41.30 

83.06 

83.06 

189 

14.33 

14.33 

24.06 

24.06 

37.93 

37.93 

42.41 

42.41 

49.20 

49.20 

190 

18.22 

18.22 

32.23 

32.23 

33.01 

33.01 

34.82 

34.82 

49.12 

49.12 

74 


Reflected 
Angle  in  ° 

Incidence  Angle  in  ° 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

191 

24.67 

24.67 

23.26 

23.26 

23.02 

23.02 

36.05 

36.05 

42.32 

42.32 

192 

17.57 

17.57 

26.96 

26.96 

19.98 

19.98 

33.80 

33.80 

43.83 

43.83 

193 

19.73 

19.73 

18.55 

18.55 

26.95 

26.95 

27.81 

27.81 

39.86 

39.86 

194 

15.90 

15.90 

24.28 

24.28 

24.10 

24.10 

30.81 

30.81 

25.96 

25.96 

195 

22.61 

22.61 

22.61 

22.61 

16.25 

16.25 

29.32 

29.32 

26.83 

26.83 

196 

11.71 

11.71 

20.60 

20.60 

21.07 

21.07 

30.66 

30.66 

34.41 

34.41 

197 

16.66 

16.66 

2.21 

2.21 

22.57 

22.57 

31.86 

31.86 

31.80 

31.80 

198 

19.92 

19.92 

17.74 

17.74 

18.11 

18.11 

29.69 

29.69 

22.77 

22.77 

199 

7.48 

7.48 

14.07 

14.07 

16.59 

16.59 

23.17 

23.17 

29.77 

29.77 

200 

8.04 

8.04 

20.62 

20.62 

12.40 

12.40 

11.70 

11.70 

26.32 

26.32 

201 

16.45 

16.45 

9.21 

9.21 

8.05 

8.05 

16.75 

16.75 

16.01 

16.01 

202 

17.43 

17.43 

10.46 

10.46 

17.62 

17.62 

22.39 

22.39 

21.75 

21.75 

203 

14.45 

14.45 

16.91 

16.91 

16.50 

16.50 

20.13 

20.13 

22.61 

22.61 

204 

3.37 

3.37 

18.38 

18.38 

22.88 

22.88 

19.21 

19.21 

14.36 

14.36 

205 

17.14 

17.14 

22.68 

22.68 

24.40 

24.40 

15.20 

15.20 

8.64 

8.64 

206 

18.44 

18.44 

21.96 

21.96 

18.99 

18.99 

19.07 

19.07 

17.71 

17.71 

207 

17.47 

17.47 

11.83 

11.83 

10.29 

10.29 

22.84 

22.84 

18.80 

18.80 

208 

19.00 

19.00 

21.17 

21.17 

17.39 

17.39 

14.10 

14.10 

23.08 

23.08 

209 

15.51 

15.51 

8.20 

8.20 

10.79 

10.79 

18.46 

18.46 

18.48 

18.48 

210 

13.66 

13.66 

20.60 

20.60 

22.07 

22.07 

23.80 

23.80 

27.02 

27.02 

211 

19.13 

19.13 

8.35 

8.35 

16.75 

16.75 

20.60 

20.60 

18.58 

18.58 

212 

16.43 

16.43 

12.76 

12.76 

15.63 

15.63 

27.42 

27.42 

24.24 

24.24 

213 

16.80 

16.80 

16.04 

16.04 

15.93 

15.93 

22.22 

22.22 

17.35 

17.35 

214 

13.54 

13.54 

18.88 

18.88 

18.29 

18.29 

14.52 

14.52 

14.30 

14.30 

215 

15.00 

15.00 

17.92 

17.92 

24.64 

24.64 

2.38 

2.38 

13.60 

13.60 

216 

12.31 

12.31 

22.16 

22.16 

15.44 

15.44 

8.73 

8.73 

13.77 

13.77 

217 

17.63 

17.63 

19.36 

19.36 

13.22 

13.22 

4.40 

4.40 

13.54 

13.54 

218 

17.13 

17.13 

16.26 

16.26 

2.58 

2.58 

7.52 

7.52 

10.42 

10.42 

219 

25.13 

25.13 

12.47 

12.47 

-2.71 

-2.71 

14.27 

14.27 

18.56 

18.56 

220 

11.46 

11.46 

12.77 

12.77 

-0.74 

-0.74 

17.41 

17.41 

18.56 

18.56 

221 

13.62 

13.62 

4.67 

4.67 

11.39 

11.39 

18.06 

18.06 

18.52 

18.52 

222 

14.73 

14.73 

5.61 

5.61 

14.56 

14.56 

19.44 

19.44 

13.22 

13.22 

223 

14.01 

14.01 

13.29 

13.29 

19.75 

19.75 

19.19 

19.19 

13.59 

13.59 

224 

12.99 

12.99 

19.66 

19.66 

20.67 

20.67 

17.67 

17.67 

14.49 

14.49 

225 

12.31 

12.31 

20.20 

20.20 

17.90 

17.90 

10.82 

10.82 

6.52 

6.52 

226 

18.75 

18.75 

19.30 

19.30 

14.68 

14.68 

5.17 

5.17 

14.82 

14.82 

227 

20.83 

20.83 

3.06 

3.06 

15.17 

15.17 

2.18 

2.18 

12.77 

12.77 

228 

15.47 

15.47 

11.32 

11.32 

17.92 

17.92 

1.76 

1.76 

14.52 

14.52 

229 

17.75 

17.75 

14.50 

14.50 

16.50 

16.50 

11.15 

11.15 

11.96 

11.96 

230 

2.76 

2.76 

12.22 

12.22 

10.17 

10.17 

11.95 

11.95 

12.22 

12.22 

231 

12.11 

12.11 

8.65 

8.65 

-1.00 

-1.00 

7.61 

7.61 

15.03 

15.03 

232 

10.97 

10.97 

13.05 

13.05 

7.87 

7.87 

1.35 

1.35 

15.80 

15.80 

233 

14.41 

14.41 

1.19 

1.19 

8.89 

8.89 

11.83 

11.83 

13.59 

13.59 

234 

2.83 

2.83 

4.93 

4.93 

1.01 

1.01 

12.37 

12.37 

10.72 

10.72 

235 

8.00 

8.00 

4.25 

4.25 

4.50 

4.50 

5.05 

5.05 

14.45 

14.45 

236 

3.85 

3.85 

8.27 

8.27 

6.18 

6.18 

7.26 

7.26 

10.84 

10.84 

237 

3.61 

3.61 

-3.69 

-3.69 

-2.64 

-2.64 

-2.80 

-2.80 

8.14 

8.14 

238 

1.33 

1.33 

4.27 

4.27 

8.60 

8.60 

3.87 

3.87 

8.86 

8.86 

239 

-9.57 

-9.57 

-0.57 

-0.57 

-0.70 

-0.70 

8.61 

8.61 

5.38 

5.38 

240 

0.19 

0.19 

9.39 

9.39 

6.38 

6.38 

7.85 

7.85 

4.06 

4.06 

75 


341 

13.75 

13.75 

14.04 

14.04 

16.69 

16.69 

10.31 

10.31 

13.62 

13.62 

342 

16.64 

16.64 

14.48 

14.48 

13.24 

13.24 

13.56 

13.56 

14.19 

14.19 

343 

16.88 

16.88 

18.26 

18.26 

18.09 

18.09 

17.68 

17.68 

24.58 

24.58 

344 

16.87 

16.87 

11.88 

11.88 

20.61 

20.61 

13.88 

13.88 

19.38 

19.38 

345 

8.19 

8.19 

14.83 

14.83 

19.87 

19.87 

25.56 

25.56 

16.03 

16.03 

346 

17.94 

17.94 

24.24 

24.24 

13.24 

13.24 

21.75 

21.75 

16.67 

16.67 

347 

17.19 

17.19 

7.26 

7.26 

24.64 

24.64 

5.61 

5.61 

19.10 

19.10 

348 

24.78 

24.78 

20.85 

20.85 

17.46 

17.46 

14.20 

14.20 

22.48 

22.48 

349 

18.67 

18.67 

25.25 

25.25 

17.22 

17.22 

17.89 

17.89 

24.02 

24.02 

350 

20.67 

20.67 

9.57 

9.57 

16.23 

16.23 

15.73 

15.73 

29.56 

29.56 

351 

24.58 

24.58 

10.26 

10.26 

17.52 

17.52 

24.12 

24.12 

37.39 

37.39 

352 

11.74 

11.74 

9.75 

9.75 

16.80 

16.80 

33.99 

33.99 

64.51 

64.51 

353 

15.73 

15.73 

21.41 

21.41 

25.53 

25.53 

36.90 

36.90 

38.93 

38.93 

354 

15.63 

15.63 

18.43 

18.43 

36.39 

36.39 

64.58 

64.58 

31.57 

31.57 

355 

18.67 

18.67 

26.36 

26.36 

38.28 

38.28 

31.76 

31.76 

22.49 

22.49 

356 

14.94 

14.94 

33.68 

33.68 

64.43 

64.43 

28.69 

28.69 

11.89 

11.89 

357 

26.43 

26.43 

37.85 

37.85 

35.23 

35.23 

22.89 

22.89 

19.38 

19.38 

358 

28.28 

28.28 

64.52 

64.52 

31.68 

31.68 

15.21 

15.21 

6.34 

6.34 

359 

37.85 

37.85 

35.69 

35.69 

23.48 

23.48 

14.92 

14.92 

11.42 

11.42 

360 

64.43 

64.43 

28.41 

28.41 

13.76 

13.76 

12.60 

12.60 

4.39 

4.39 

inciaence 


Angle  in 
16 

16.04 

16.51 
18.11 

13.37 
14.80 
17.74 
16.02 
14.73 

18.93 
15.42 

16.55 

13.21 

15.48 

16.37 
17.13 
14.15 
16.46 
16.87 

16.48 
16.41 

15.61 
13.99 

13.83 
14.77 
15.09 
14.45 
15.00 
12.30 
14.58 
15.82 

14.84 

15.21 

14.64 

12.64 
13.26 
13.89 

11.93 

12.51 

13.57 

12.22 
10.86 
12.70 

12.55 
12.05 

11.55 

9.57 

5.61 
6.90 
2.11 
7.28 


inciaence 


Angle  in 
16 
2.82 

2.24 

4.26 

4.68 
5.07 

10.21 

9.62 
12.28 
12.88 

15.62 

13.84 
14.31 

16.53 

15.64 
15.96 
18.09 

17.85 

14.26 
7.09 

16.21 

23.71 
27.90 
28.33 
25.83 
28.08 

29.72 
25.49 
20.42 
24.00 
18.01 

13.21 
8.77 
15.79 

11.86 

11.53 
8.35 
12.98 
12.06 
14.74 
10.10 
7.80 

10.63 
10.93 

8.64 

6.25 
6.89 

7.69 
11.04 
4.29 
14.60 


inciaence 


Angle  in 
16 
9.91 
13.05 
11.40 
12.71 

12.47 
13.28 
11.81 
12.22 

12.97 
12.03 

10.56 
11.24 
13.32 

13.18 

7.29 

8.97 
-6.16 
17.84 

18.99 

3.99 
-3.29 
5.94 
2.75 

11.19 
6.89 
8.59 

25.48 
11.18 

13.99 
8.74 
10.95 
12.76 
9.26 
5.08 
11.09 

7.42 

9.42 
16.23 
12.11 
9.07 

4.56 
12.86 
12.45 

9.43 
17.31 

7.30 
18.18 
18.92 
18.05 
20.87 


Reflected 
Angle  in  ° 

Incidence  Angle  in  ° 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

151 

9.66 

9.66 

14.74 

14.74 

16.43 

16.43 

17.13 

17.13 

15.55 

15.55 

152 

10.69 

10.69 

11.68 

11.68 

16.60 

16.60 

15.64 

15.64 

13.44 

13.44 

153 

18.65 

18.65 

18.04 

18.04 

18.21 

18.21 

17.20 

17.20 

21.33 

21.33 

154 

18.47 

18.47 

11.43 

11.43 

17.83 

17.83 

11.98 

11.98 

23.36 

23.36 

155 

20.81 

20.81 

8.87 

8.87 

8.40 

8.40 

14.01 

14.01 

27.01 

27.01 

156 

16.17 

16.17 

14.74 

14.74 

19.41 

19.41 

29.56 

29.56 

26.58 

26.58 

157 

14.37 

14.37 

11.35 

11.35 

24.99 

24.99 

32.88 

32.88 

32.84 

32.84 

158 

21.56 

21.56 

19.58 

19.58 

26.64 

26.64 

32.65 

32.65 

36.11 

36.11 

159 

24.22 

24.22 

22.44 

22.44 

30.23 

30.23 

29.60 

29.60 

39.23 

39.23 

160 

18.62 

18.62 

26.54 

26.54 

31.02 

31.02 

29.36 

29.36 

38.50 

38.50 

161 

19.23 

19.23 

24.40 

24.40 

33.88 

33.88 

34.25 

34.25 

41.75 

41.75 

162 

27.19 

27.19 

21.68 

21.68 

31.56 

31.56 

35.21 

35.21 

54.47 

54.47 

163 

30.19 

30.19 

28.26 

28.26 

34.43 

34.43 

40.25 

40.25 

43.28 

43.28 

164 

28.93 

28.93 

33.02 

33.02 

30.40 

30.40 

52.28 

52.28 

36.20 

36.20 

165 

33.03 

33.03 

34.80 

34.80 

36.93 

36.93 

40.23 

40.23 

28.58 

28.58 

166 

28.18 

28.18 

37.36 

37.36 

46.81 

46.81 

36.30 

36.30 

30.53 

30.53 

167 

34.73 

34.73 

37.38 

37.38 

39.51 

39.51 

32.44 

32.44 

34.91 

34.91 

168 

40.28 

40.28 

47.95 

47.95 

41.66 

41.66 

21.85 

21.85 

35.35 

35.35 

169 

39.35 

39.35 

40.01 

40.01 

36.37 

36.37 

32.66 

32.66 

25.63 

25.63 

170 

50.69 

50.69 

40.69 

40.69 

28.49 

28.49 

31.44 

31.44 

31.03 

31.03 

171 

37.33 

37.33 

37.37 

37.37 

35.07 

35.07 

32.70 

32.70 

26.56 

26.56 

172 

37.94 

37.94 

35.89 

35.89 

32.91 

32.91 

27.00 

27.00 

32.41 

32.41 

173 

37.27 

37.27 

36.60 

36.60 

34.32 

34.32 

31.23 

31.23 

31.82 

31.82 

174 

34.86 

34.86 

28.22 

28.22 

30.51 

30.51 

32.62 

32.62 

17.80 

17.80 

175 

32.68 

32.68 

32.99 

32.99 

33.10 

33.10 

34.89 

34.89 

32.87 

32.87 

176 

31.38 

31.38 

29.22 

29.22 

23.01 

23.01 

31.49 

31.49 

22.49 

22.49 

177 

33.01 

33.01 

26.73 

26.73 

35.09 

35.09 

32.50 

32.50 

32.50 

32.50 

178 

22.70 

22.70 

33.06 

33.06 

35.15 

35.15 

11.79 

11.79 

32.61 

32.61 

179 

31.53 

31.53 

37.49 

37.49 

32.11 

32.11 

30.74 

30.74 

23.67 

23.67 

180 

29.61 

29.61 

35.13 

35.13 

30.06 

30.06 

34.77 

34.77 

27.86 

27.86 

181 

32.42 

32.42 

31.90 

31.90 

27.43 

27.43 

37.66 

37.66 

26.77 

26.77 

182 

33.26 

33.26 

31.96 

31.96 

32.34 

32.34 

39.75 

39.75 

28.89 

28.89 

183 

28.02 

28.02 

38.76 

38.76 

38.59 

38.59 

38.83 

38.83 

29.97 

29.97 

184 

42.72 

42.72 

41.91 

41.91 

41.66 

41.66 

26.19 

26.19 

35.39 

35.39 

185 

45.36 

45.36 

35.09 

35.09 

37.95 

37.95 

38.96 

38.96 

37.70 

37.70 

186 

39.30 

39.30 

40.98 

40.98 

38.72 

38.72 

41.36 

41.36 

30.89 

30.89 

187 

40.59 

40.59 

39.56 

39.56 

41.12 

41.12 

38.74 

38.74 

38.75 

38.75 

188 

50.11 

50.11 

46.48 

46.48 

38.01 

38.01 

41.69 

41.69 

27.07 

27.07 

189 

50.90 

50.90 

37.73 

37.73 

42.80 

42.80 

43.06 

43.06 

42.86 

42.86 

190 

83.95 

83.95 

37.84 

37.84 

36.81 

36.81 

41.47 

41.47 

41.85 

41.85 

191 

55.81 

55.81 

58.30 

58.30 

45.54 

45.54 

39.89 

39.89 

41.62 

41.62 

192 

40.31 

40.31 

84.69 

84.69 

52.59 

52.59 

45.87 

45.87 

43.66 

43.66 

193 

44.33 

44.33 

60.30 

60.30 

60.51 

60.51 

43.02 

43.02 

44.69 

44.69 

194 

32.20 

32.20 

51.64 

51.64 

85.36 

85.36 

44.65 

44.65 

42.19 

42.19 

195 

33.14 

33.14 

41.08 

41.08 

55.51 

55.51 

50.19 

50.19 

39.39 

39.39 

196 

39.00 

39.00 

44.49 

44.49 

46.14 

46.14 

85.92 

85.92 

51.24 

51.24 

197 

35.15 

35.15 

43.09 

43.09 

35.30 

35.30 

57.27 

57.27 

60.06 

60.06 

198 

33.91 

33.91 

41.61 

41.61 

37.11 

37.11 

49.93 

49.93 

86.58 

86.58 

199 

28.76 

28.76 

37.96 

37.96 

41.21 

41.21 

43.63 

43.63 

60.17 

60.17 

200 

24.61 

24.61 

18.10 

18.10 

28.91 

28.91 

46.88 

46.88 

40.60 

40.60 

82 


Reflected 
Angle  in  ° 

Incidence  Angle  in  ° 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

201 

26.93 

26.93 

28.86 

28.86 

33.44 

33.44 

37.48 

37.48 

46.33 

46.33 

202 

19.69 

19.69 

23.90 

23.90 

28.74 

28.74 

34.07 

34.07 

37.44 

37.44 

203 

21.25 

21.25 

19.67 

19.67 

32.42 

32.42 

27.21 

27.21 

29.36 

29.36 

204 

18.66 

18.66 

19.73 

19.73 

22.56 

22.56 

30.18 

30.18 

33.55 

33.55 

205 

10.44 

10.44 

22.40 

22.40 

24.29 

24.29 

23.99 

23.99 

33.37 

33.37 

206 

24.83 

24.83 

30.02 

30.02 

23.95 

23.95 

10.09 

10.09 

30.32 

30.32 

207 

22.80 

22.80 

25.21 

25.21 

16.82 

16.82 

18.88 

18.88 

31.67 

31.67 

208 

25.08 

25.08 

25.10 

25.10 

20.90 

20.90 

22.73 

22.73 

17.32 

17.32 

209 

25.88 

25.88 

18.85 

18.85 

16.69 

16.69 

17.23 

17.23 

27.69 

27.69 

210 

4.64 

4.64 

5.57 

5.57 

5.99 

5.99 

27.32 

27.32 

19.74 

19.74 

211 

23.53 

23.53 

16.59 

16.59 

15.15 

15.15 

14.88 

14.88 

19.19 

19.19 

212 

19.98 

19.98 

18.79 

18.79 

12.65 

12.65 

14.87 

14.87 

23.79 

23.79 

213 

20.63 

20.63 

21.88 

21.88 

22.70 

22.70 

21.13 

21.13 

14.20 

14.20 

214 

19.24 

19.24 

23.21 

23.21 

18.60 

18.60 

22.92 

22.92 

12.72 

12.72 

215 

22.14 

22.14 

25.15 

25.15 

22.06 

22.06 

21.86 

21.86 

9.77 

9.77 

216 

22.93 

22.93 

20.15 

20.15 

24.08 

24.08 

13.52 

13.52 

18.52 

18.52 

217 

23.33 

23.33 

22.29 

22.29 

24.36 

24.36 

20.28 

20.28 

23.44 

23.44 

218 

22.53 

22.53 

21.73 

21.73 

18.82 

18.82 

19.74 

19.74 

17.58 

17.58 

219 

17.18 

17.18 

20.49 

20.49 

20.58 

20.58 

20.87 

20.87 

21.98 

21.98 

220 

15.55 

15.55 

18.00 

18.00 

19.70 

19.70 

18.15 

18.15 

21.78 

21.78 

221 

15.22 

15.22 

12.67 

12.67 

14.45 

14.45 

13.46 

13.46 

10.40 

10.40 

222 

10.03 

10.03 

16.70 

16.70 

7.56 

7.56 

16.08 

16.08 

16.38 

16.38 

223 

17.08 

17.08 

9.66 

9.66 

8.92 

8.92 

12.97 

12.97 

15.80 

15.80 

224 

13.40 

13.40 

5.77 

5.77 

9.89 

9.89 

15.47 

15.47 

16.04 

16.04 

225 

-2.65 

-2.65 

8.37 

8.37 

13.50 

13.50 

15.52 

15.52 

14.49 

14.49 

226 

10.79 

10.79 

8.98 

8.98 

7.73 

7.73 

9.05 

9.05 

6.64 

6.64 

227 

13.99 

13.99 

-2.79 

-2.79 

5.46 

5.46 

15.66 

15.66 

15.59 

15.59 

228 

15.07 

15.07 

8.83 

8.83 

8.97 

8.97 

8.75 

8.75 

6.41 

6.41 

229 

13.74 

13.74 

12.52 

12.52 

12.31 

12.31 

-1.11 

-1.11 

-0.51 

-0.51 

230 

13.68 

13.68 

10.72 

10.72 

13.56 

13.56 

-1.10 

-1.10 

4.49 

4.49 

231 

11.18 

11.18 

11.23 

11.23 

16.52 

16.52 

8.99 

8.99 

-9.12 

-9.12 

232 

13.41 

13.41 

15.49 

15.49 

14.13 

14.13 

14.54 

14.54 

4.20 

4.20 

233 

15.07 

15.07 

15.89 

15.89 

8.99 

8.99 

10.30 

10.30 

-6.03 

-6.03 

234 

12.57 

12.57 

15.55 

15.55 

15.78 

15.78 

11.79 

11.79 

10.53 

10.53 

235 

10.26 

10.26 

14.60 

14.60 

15.43 

15.43 

14.90 

14.90 

9.30 

9.30 

236 

14.67 

14.67 

15.61 

15.61 

15.45 

15.45 

13.36 

13.36 

17.97 

17.97 

237 

11.07 

11.07 

13.10 

13.10 

12.80 

12.80 

12.45 

12.45 

13.73 

13.73 

238 

9.84 

9.84 

12.74 

12.74 

13.75 

13.75 

14.33 

14.33 

9.68 

9.68 

239 

8.08 

8.08 

13.99 

13.99 

15.70 

15.70 

18.14 

18.14 

17.22 

17.22 

240 

11.31 

11.31 

10.63 

10.63 

15.28 

15.28 

14.58 

14.58 

16.39 

16.39 

241 

8.86 

8.86 

11.34 

11.34 

11.11 

11.11 

14.97 

14.97 

13.57 

13.57 

242 

6.55 

6.55 

12.75 

12.75 

5.86 

5.86 

14.02 

14.02 

12.58 

12.58 

243 

6.43 

6.43 

3.37 

3.37 

9.60 

9.60 

10.88 

10.88 

15.98 

15.98 

244 

-4.93 

-4.93 

10.53 

10.53 

8.47 

8.47 

14.51 

14.51 

14.03 

14.03 

245 

-3.66 

-3.66 

3.64 

3.64 

5.07 

5.07 

14.81 

14.81 

13.13 

13.13 

246 

14.16 

14.16 

1.68 

1.68 

2.73 

2.73 

11.55 

11.55 

14.45 

14.45 

247 

10.29 

10.29 

12.28 

12.28 

3.31 

3.31 

8.65 

8.65 

9.94 

9.94 

248 

14.57 

14.57 

1.15 

1.15 

-1.60 

-1.60 

9.87 

9.87 

-1.54 

-1.54 

249 

12.25 

12.25 

7.31 

7.31 

2.54 

2.54 

9.24 

9.24 

6.46 

6.46 

250 

11.37 

11.37 

9.89 

9.89 

8.39 

8.39 

5.66 

5.66 

12.62 

12.62 
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Reflected 
Angle  in  ° 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 


Incidence  Angle  in 


13.07  13.07 

17.40  17.40 

16.33  16.33 

0.44  0.44 

15.91  15.91 

14.04  14.04 

16.87  16.87 

8.33  8.33 

13.68  13.68 

5.33  5.33 

4.07  4.07 

12.84  12.84 

14.19  14.19 

12.99  12.99 

10.99  10.99 

0.19  0.19 

11.47  11.47 

6.87  6.87 

15.50  15.50 

7.87  7.87 

10.11  10.11 

12.77  12.77 

14.92  14.92 

13.09  13.09 

15.00  15.00 

12.31  12.31 

12.48  12.48 

14.98  14.98 

13.82  13.82 

14.38  14.38 

13.09  13.09 

9.71  9.71 

12.76  12.76 

10.54  10.54 

8.82  8.82 

10.85  10.85 

11.45  11.45 

11.49  11.49 

11.72  11.72 

12.86  12.86 

11.79  11.79 

10.42  10.42 

10.41  10.41 

7.96  7.96 


inciaence 


Angle  in 
16 
3.72 
6.87 
7.69 
8.10 
5.38 

14.62 

14.25 

8.95 
14.35 
13.52 
14.64 

15.24 

14.77 
14.48 
16.22 
16.12 
14.41 

12.88 

13.37 
18.03 

17.25 
16.04 

15.89 

15.76 

17.31 
15.51 
12.60 

13.63 
14.39 
10.33 
9.49 

15.32 

20.78 
19.10 

19.25 

19.77 

21.37 

11.90 
16.50 
14.43 

14.27 
22.16 
24.57 

36.95 
41.54 
63.71 
35.23 
31.46 

23.28 

16.32 


APPENDIX  F.  MATLAB  SOURCE  CODE  FOR  DETECTION 

COVERAGE  PLOT 


%  Chong  Sze 

Sing  (March  2014) 

%  Bistatic 

RCS  Look-up  Table 

load ( ' 3GHz 

full .mat ' ) 

biRCS3  =  num; 

load ( ' 9GHz 

full .mat ' ) 

biRCS9  =  num; 

sen=10;  %  sensitivity 

%  For  Tx  1 

Xcoordl  =  [ 

] ; 

Ycoordl  =  [ 

] ; 

SNRl  =  [ ]  ; 

SNRdBml  =  [ 

] ; 

Err  thetaRl 

=  []; 

Err  timel  = 

[] ; 

Err  LI  =  [] 

r 

Txl=45; 

Ll=18000; 

%For  Tx2 

Xcoord2  =  [ 

] ; 

Ycoord2  =  [ 

] ; 

SNR2  =  [ ]  ; 

SNRdBm2  =  [ 

] ; 

Err  thetaR2 

=  []; 

Err  time2  = 

[] ; 

Err  L2  =  [] 

r 

Tx2=135; 

L2=18000; 

%For  Tx3 

Xcoord3  =  [ 

] ; 

Ycoord3  =  [ 

] ; 

SNR3  =  [ ] ; 

SNRdBm3  =  [ 

] ; 

Err  thetaR3 

=  []; 

Err  time3  = 

[] ; 

Err  L3  =  [] 

r 

Tx3=225; 

L3=18000; 

%For  Tx4 

Xcoord4  =  [ 

] ; 

Ycoord4  =  [ 

] ; 

SNR4  =  [ ]  ; 

SNRdBm4  =  [ 

] ; 

Err  thetaR4 

=  []; 

Err  time4  = 

[] ; 

Err  L4  =  [] 

r 

Tx4=315; 

L4=18000; 
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%For  Tx5 
XcoordS  =  [ ] ; 

YcoordS  =  [ ] ; 

SNR5  =  [ ] ; 

SNRdBmS  =  [ ] ; 

Err_thetaR5  =  [ ] ; 

Err_time5  =  [ ] ; 

Err_L5  =  [ ] ; 

Tx5=0; 

L5=18000; 

%For  Tx6 
XcoordG  =  [ ]  ; 

YcoordG  =  [ ] ; 

SNR6  =  [ ] ; 

SNRdBme  =  [ ] ; 

Err_thetaR6  =  [ ] ; 

Err_time6  =  [ ] ; 

Err_L6  =  [ ] ; 

Tx6=90; 

L6=18000; 

%For  Tx7 
Xcoord7  =  [ ] ; 

Ycoord7  =  [ ] ; 

SNR7  =  [ ] ; 

SNRdBm7  =  [ ] ; 

Err_thetaR7  =  [ ] ; 

Err_time7  =  [ ] ; 

Err_L7  =  [ ] ; 

Tx7=180; 

L7=18000; 

%For  Tx8 
Xcoord8  =  [ ] ; 

Ycoord8  =  [ ] ; 

SNR8  =  [ ] ; 

SNRdBm8  =  [ ] ; 

Err_thetaR8  =  [ ] ; 

Err_time8  =  [ ] ; 

Err_L8  =  [ ] ; 

Tx8=270; 

L8=18000; 

%  Bi-static  Radar  Equation 
%  Constants 

k  =  1.38e-23;  %  Boltzmann 

c  =  3e8;  %  Speed  of  light 

%  Variables 

%  Parameters  are  based  Furuno  Marine  Radar  FAR-28x7  series  in  both  X-band 
and  S-band 

%  Bi-static  Radar  Equation 
%  Constants 

k  =  1.38e-23;  %  Boltzmann 

c  =  3e8;  %  Speed  of  light 
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Select  =  1;  %chose  which  radar  to  use 
%  Variables 

%  Parameters  are  based  Furuno  Marine  Radar  FAR-28x7  series  in  both  X-band 
%  and  S-band 


%S-band  Radar  (1) 
BeamAz_S  =  1.8; 
BeamEL_S  =  25; 
Scanrate_S  =  45; 
Pt_S  =  30e3; 

Gt_S  =  10^  (28/10) 
Gr_S  =  2; 
f_S  =  3050e6; 
lambda_S  =  c/f_S; 
Ft  S  =  1; 


%  in  degrees 
%  in  degrees 

%  in  rpm  for  range  of  24  nmi  or  higher 
%  Tx  peak  power  output 
%  Tx  Antenna  Gain  given  as  28dB 
%  Rx  Antenna  Gain 
%  frequency 

%  wavelength 

%  Prop  Factor  for  Tx-Tgt  path.  See  Advances  in  Bistatic 


Radar  -  Chp  6. 5. 5. 6 

Fr_S  =1;  %  Prop  Factor  for  Rx-Tgt  path 

Ts  S  =  290+290*  (13 . 5-1) ;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using 

formula  Te=(NF-l)290  and  Ts=Te+290K 

Bn  S  =  500e6;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 

MHz) 


Lt_S  =  2;  %  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp  6. 5. 5. 6 

Lr_S  =  2; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 

PRF_S  =  2000; 

TOT_S  =  BeamAz_S/ (Scanrate_S*360) ; 

N  S  =  TOT  S*PRF  S; 


%X-band  Radar  (2) 
BeamAz_X  =  0.95; 
BeamEL_X  =  25; 
Scanrate_X  =  45; 
Pt_X  =  25e3; 

Gt_X  =  lO'^  (31/10) 
Gr_X  =  2; 
f_X  =  9410e6; 
lambda_X  =  c/f_X; 
Ft  X  =  1; 


%  in  degrees 
%  in  degrees 

%  in  rpm  for  range  of  24  nmi  or  higher 
%  Tx  peak  power  output 
%  Tx  Antenna  Gain  is  given  as  31  dB 
%  Rx  Antenna  Gain 
%  frequency 

%  wavelength 

%  Prop  Factor  for  Tx-Tgt  path.  See  Advances  in  Bistatic 


Radar  -  Chp  6. 5. 5. 6 

Fr_X  =1;  %  Prop  Factor  for  Rx-Tgt  path 

Ts_X  =  290+290*  (6-1) ;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using 

formula  Te=(NF-l)290  and  Ts=Te+290K 

Bn  X  =  500e6;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 

MHz) 


Lt_X  =  2;  %  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp  6. 5. 5. 6 

Lr_X  =  2; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 

PRF_X  =  1500; 

TOT_X  =  BeamAz_X/ (Scanrate_X*360) ; 

N  X  =  TOT  X*PRF  X; 
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if  Select  ==1 

BeamEL  =  BeamEL_S;  %  in  degrees 
BeamAz  =  BeamAz  S;  %  in  degrees 

Scanrate  =  Scanrate_S;  %  in  rpm  for  range  of  24  nmi  or  higher 

Ft  =  Pt_S;  %  Tx  peak  power  output 

Gt  =  Gt_S;  %  Tx  Antenna  Gain  is  given  as  31  dB 

Gr  =  Gr_S;  %  Rx  Antenna  Gain 

f  =  f_S;  %  frequency 

lambda  =  lambda_S;  %  wavelength 

Ft  =  Ft_S;  %  Prop  Factor  for  Tx-Tgt  path.  See  Advances  in 

Bistatic  Radar  -  Chp  6. 5. 5. 6 

Fr  =  Fr_S;  %  Prop  Factor  for  Rx-Tgt  path 

Ts  =  Ts  S;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using  formula 

Te=(NF-l)290  and  Ts=Te+290K 

Bn  =  Bn_S;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 

MHz) 

Lt  =  Lt_S;  %  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp 

6. 5. 5. 6 


Lr  =  Lr_S ; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 
PRF  =  PRF_S; 

TOT  =  TOT_S; 

N  =  N_S; 

else 


BeamEL  = 
BeamAz  = 
Scanrate 
Pt  =  Pt_X 
Gt  =  Gt_X 
Gr  =  Gr_X 
f  =  f  X; 


BeamEL_X;  %  in  degrees 
BeamAz_X;  %  in  degrees 
=  Scanrate_X;  %  in  rpm  for  range 
;  %  Tx  peak  power  output 

;  %  Tx  Antenna  Gain  is  given  as 

;  %  Rx  Antenna  Gain 

%  frequency 


of  24  nmi 
31  dB 


or  higher 


lambda  =  lambda_X;  %  wavelength 

Ft  =  Ft_X;  %  Prop  Factor  for  Tx-Tgt  path.  See  Advances 


in 


Bistatic  Radar  -  Chp  6. 5. 5. 6 

Fr  =  Fr_X;  %  Prop  Factor  for  Rx-Tgt  path 

Ts  =  Ts_X;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using  formula 

Te=(NF-l)290  and  Ts=Te+290K 

Bn  =  Bn  X;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 


MHz) 


Lt  =  Lt_X; 
6. 5. 5. 6 


%  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp 


Lr  =  Lr_X; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 
PRF  =  PRF_X; 

TOT  =  TOT_X; 

N  =  N_X; 

end 
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9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  TKTPTTT  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 
ooooooooooooooooooooo  ilNlrUJ-  oooooooooooooooooooo 

%  We  only  take  in  three  measurments  from  the  ESM  -  L,  Rt+Rr  and  thetaR 

9-9-9-9-9-9-9-9'9-9-9-9-9-9-9-9'9-9-9-9-9-  Tv  1  9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9- 

000000000000000000000  lA  ±  000000000000000000000000000 


for  CombinR=18500 : 200 : 50000 


for  thetaR=0 : 0 . 5 : 360 


RrAct=  (  (CombinR)  ''2  - 

1j1''2)  /  (2*  (CombinR+Ll*sin  (thetaR*pi/180)  )  )  ; 

RtAct=sqrt (RrAct^2+L1^2-2*RrAct*Ll*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round ( (180/pi) *acos ( (RrAct^2+RtAct^2- 
Ll''2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 


%  Organise  the  RCS  values  by 

RCS=10''  (biRCS3  (90-Betal  +  l,  90+Betal)  /lO)  ; 

%RCS=10^ (biRCS3 (Betal  +  1 ,  Betaltl )  /  lO)  ; 

K1  =  Pt*Gt*Gr* (lambda^2) *RCS* (Ft^2) * (Fr^2) ; 
K2  =  (  (4*pi)  ^'3)  *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=Ll/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct^2*RtAct^2) ;  % 

SNRdBm  =  10*log(SNR)+  sqrt (N) ; 

Xcoord  =  RrAct*sin ( (thetaR+Txl) *pi/180)  ; 
Ycoord  =  RrAct*cos ( (thetaR+Txl) *pi/180)  ; 


edge  of  the 
as 


%  Error  to  measurement  of  Rt  +  Rr 
%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 
%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 

%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e''2+2*e*sin  (thetaR*pi/180)  )  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ^2)  ; 

Err_time=Err_factor_time* (BeamAz/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  60nmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e^2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e''2)  *cos  (thetaR*pi/180)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ^2)  ; 

Err_thetaR  =  Err_factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoordl  =  [Xcoordl , Xcoord] ; 

Ycoordl  =  [ Ycoordl , Ycoord] ; 

SNRl  =  [SNR1,SNR]; 

SNRdBml  =  [ SNRdBml , SNRdBm] ; 

Err_timel  =  [Err_timel , Err_time] ; 

Err_Ll  =  [Err_Ll , Err_L] ; 

Err_thetaRl  =  [Err_thetaRl , Err_thetaR] ; 

else 

end 

end 

end 

Err  Totall=zeros (size (Err_timel) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err_timel ) 

Err_Totall (aa)  = 

max (abs (Err_timel (aa) ) , abs (Err_Ll (aa) ) +abs (Err_thetaRl (aa) ) ) ; 
end 

for  aa=l : 1 : length (Err_timel ) 

Err_Totall (aa)  =  max (abs (Err_Totall (aa) ) , abs (Err_thetaRl (aa) ) ) ; 

end 
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9-9-9-9-9-9'9-9-9-9-9'9-9-9'9-9-9-9-9-9-9-  T'-v  9  9-9-9-9'9-9-9-9'9-9-9'9-9-9-9-9'9-9-9-9-9-9-9-9-9-9-9- 

ooooooooooooooooooooo  lA  Z.  ooooooooooooooooooooooooooo 


for  CombinR=18500 : 200 : 50000 
for  thetaR=0 : 0 . 5 : 360 


RrAct=  (  (CombinR)  ''2  - 

L2^2) / (2* (CombinR+L2*sin (thetaR*pi/180)  )  )  ; 

RtAct=sqrt (RrAct^2+L2^2-2*RrAct*L2*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round  (  (180/pi)  *acos  (  (RrAct''2+RtAct''2- 
L2'^2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 


%  Organise  the  RCS  values  by 

RCS=10'^  (biRCS3  (90-Betal  +  l,  90+Betal)  /lO)  ; 

%RCS=10"'  (biRCS3  (Betal  +  1 ,  Betal  +  1 )  /lO)  ; 

K1  =  Pt*Gt*Gr*  (lambda''2)  *RCS*  (Ft^2)  *  (Fr^2)  ; 
K2  =  (  (4*pi)  •^3)  *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L2/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct^2*RtAct^2) ; 

SNRdBm  =  10*log(SNR)+  sqrt (N) ; 

Xcoord  =  RrAct*sin ( (thetaR+Tx2) *pi/180)  ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx2) *pi/180) ; 


%  Error  to  measurement  of  Rt  t  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

edge  of  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 


%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e^2+2*e*sin (thetaR*pi/180) ) / (2* (l+e*sin (thetaR*pi/180) ) ^2) ; 

Err_time=Err_factor_time* (BeamAz/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  de 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  60nmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e''2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e'*2)  *cos  (thetaR*pi/180)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_thetaR  =  Err_factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoord2  =  [Xcoord2 , Xcoord] ; 

Ycoord2  =  [ Ycoord2 , Ycoord] ; 

SNR2  =  [SNR2,SNR]; 

SNRdBm2  =  [ SNRdBm2 , SNRdBm] ; 

Err_time2  =  [Err_time2 , Err_time] ; 

Err_L2  =  [Err_L2 , Err_L] ; 

Err_thetaR2  =  [Err_thetaR2 , Err_thetaR] ; 

else 

end 

end 

end 

Err_Total2=zeros (size (Err_time2) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err_time2 ) 

Err_Total2 (aa)  =  max (abs (Err_time2 (aa) ) , abs (Err_L2 (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err_time2 ) 

Err_Total2 (aa)  =  max (abs (Err_Total2 (aa) ) , abs (Err_thetaR2 (aa) ) ) ; 

end 
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9-9-9-9-9-9-9-9-9'9-9-9'9-9-9-9-9'9-9-9'9-  T'\r  9'9-9-9-9-9'9-9-9'9'9-9'9-9'9-9-9'9-9-9-9-9'9-9-9-9-9- 
ooooooooooooooooooooo  lA  O  ooooooooooooooooooooooooooo 


for  CombinR=18500 : 200 : 50000 
for  thetaR=0 : 0 . 5 : 360 


RrAct=  (  (CombinR)  ''2  - 

L3^2) / (2* (CombinR+L3*sin (thetaR*pi/180) )  )  ; 

RtAct=sqrt (RrAct^2+L3^2-2*RrAct*L3*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round  (  (180/pi)  *acos  (  (RrAct''2+RtAct''2- 
L3'^2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 


%  Organise  the  RCS  values  by 

RCS=10'^  (biRCS3  (90-Betal  +  l,  90tBetal)  /lO)  ; 

%RCS=10''  (biRCS3  (Betaltl ,  Betal  +  1 )  /lO)  ; 

K1  =  Pt*Gt*Gr* (lambda^2) *RCS* (Ft^2) * (Fr^2) ; 
K2  =  (  (4*pi)  ■^3)  *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L3/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct^2*RtAct''2)  ; 

SNRdBm  =  10*log(SNR)+  sqrt (N) ; 

Xcoord  =  RrAct*sin ( (thetaR+Tx3) *pi/180) ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx3) *pi/180) ; 


%  Error  to  measurement  of  Rt  +  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

edge  of  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 


%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e^2+2*e*sin (thetaR*pi/180) ) / (2* (l+e*sin (thetaR*pi/180) ) ^2)  ; 

Err_time=Err_factor_time* (BeamAz/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  60nmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e^2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err  factor_theta  =  (1- 

e''2)  *cos  (thetaR*pi/180)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_thetaR  =  Err_factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoord3  =  [Xcoord3 , Xcoord] ; 

Ycoord3  =  [ Ycoord3 , Ycoord] ; 

SNR3  =  [SNR3,SNR]; 

SNRdBm3  =  [ SNRdBm3 , SNRdBm] ; 

Err_time3  =  [Err_time3, Err_time] ; 

Err_L3  =  [Err_L3 , Err_L] ; 

Err_thetaR3  =  [Err_thetaR3 , Err_thetaR] ; 

else 

end 

end 

end 

Err_Total3=zeros (size (Err_time3) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err_time3 ) 

Err_Total3 (aa)  =  max (abs (Err_time3 (aa) ) , abs (Err_L3 (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err_time3 ) 

Err_Total3 (aa)  =  max (abs (Err_Total3 (aa) ) , abs (Err_thetaR3 (aa) ) ) ; 

end 
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9-9-9-9-9-9-9-9'9-9-9-9-9-9-9-9-9-9-9-9-9-  T'\r  A  9-9-9-9-9-9-9-9-9'9-9-9-9'9-9-9-9-9-9-9-9-9-9-9-9'9-9- 

ooooooooooooooooooooo  lA  4  ooooooooooooooooooooooooooo 


for  CombinR=18500 : 200 : 45000 


for  thetaR=0 : 0 . 5 : 360 


RrAct= ( (CombinR) ^2  - 

L4'^2)  /  (2*  (CombinR+L4*sin  (thetaR*pi/180)  )  )  ; 

RtAct=sqrt (RrAct^2+L4^2-2*RrAct*L4*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round ( (180/pi) *acos ( (RrAct^2+RtAct^2- 
L4'^2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 


%  Organise  the  RCS  values  by 

RCS=10'^  (biRCS3  (90-Betal  +  l,  90+Betal)  /lO)  ; 

%RCS=10"'  (biRCS3  (Betal  +  1 ,  Betal  +  1 )  /lO)  ; 


K1  =  Pt*Gt*Gr*  (lambda''2)  *RCS*  (Ft^2)  *  (Fr^2)  ; 
K2  =  (  (4*pi)  •^3)  *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L4/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct^2*RtAct^2) ; 

SNRdBm  =  10*log(SNR)+  sqrt (N) ; 

Xcoord  =  RrAct*sin ( (thetaR+Tx4) *pi/180) ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx4) *pi/180) ; 


edge  of  the 
as 


%  Error  to  measurement  of  Rt  +  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 

%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e'^2+2*e*sin  (thetaR*pi/180)  )  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_time=Err_factor_time* (BeamAz/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err  time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  GOnmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e''2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e'^2)  *cos  (thetaR*pi/180)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_thetaR  =  Err_factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoord4  =  [Xcoord4 , Xcoord] ; 

Ycoord4  =  [ Ycoord4 , Ycoord] ; 

SNR4  =  [SNR4,SNR]; 

SNRdBm4  =  [ SNRdBm4 , SNRdBm] ; 

Err_time4  =  [Err_time4 , Err_time] ; 

Err_L4  =  [Err_L4 , Err_L] ; 

Err_thetaR4  =  [Err_thetaR4 , Err_thetaR] ; 

else 

end 

end 

end 

Err_Total4=zeros (size (Err_time4) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err_time4 ) 

Err_Total4 (aa)  =  max (abs (Err_time4 (aa) ) , abs (Err_L4 (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err  time!) 

Err_Total4 (aa)  =  max (abs (Err_Total4 (aa) ) , abs (Err_thetaR4 (aa) ) ) ; 

end 
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9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-  T'-v  R  9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9- 
ooooooooooooooooooooo  lA  O  ooooooooooooooooooooooooooo 


for  CombinR=18500 : 200 : 50000 
for  thetaR=0 : 0 . 5 : 360 


RrAct=  (  (CombinR)  ''2  - 

L5^2) / (2* (CombinR+L5*sin (thetaR*pi/180) ) ) ; 

RtAct=sqrt (RrAct^2+L5^2-2*RrAct*L5*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round  (  (180/pi)  *acos  (  (RrAct^2+RtAct''2- 
L5’^2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 

%  Organise  the  RCS  values  by 
RCS=10''  (biRCS3  (90-Betal,  90+Betal)  /lO)  ; 

K1  =  Pt*Gt*Gr* (lambda^2) *RCS* (Ft^2) * (Fr^2) ; 

K2  =  (  (4*pi)  ■^3)  *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L5/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct^2*RtAct^2) ; 

SNRdBm  =  10*log (SNR) +30; 

Xcoord  =  RrAct*sin ( (thetaR+Tx5) *pi/180) ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx5) *pi/180) ; 


edge  of  the 
as 


%  Error  to  measurement  of  Rt  +  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 

%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e'^2+2*e*sin  (thetaR*pi/180)  )  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_time=Err_factor_time* (BeamEL/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  lib) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  60nmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e''2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e'^2)  *cos  (thetaR*pi/180)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_thetaR  =  Err_factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoord5  =  [Xcoord5 , Xcoord] ; 

Ycoord5  =  [ Ycoord5 , Ycoord] ; 

SNR5  =  [SNR5,SNR]; 

SNRdBm5  =  [ SNRdBm5 , SNRdBm] ; 

Err_time5  =  [Err_time5, Err_time] ; 

Err_L5  =  [Err_L5 , Err_L] ; 

Err_thetaR5  =  [Err_thetaR5 , Err_thetaR] ; 

else 

end 

end 

end 

Err_Total5=zeros (size (Err  time5) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err  time5) 

Err_Total5 (aa)  =  max (abs (Err_time5 (aa) ) , abs (Err_L5 (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err  time5) 

Err_Total5 (aa)  =  max (abs (Err_Total5 (aa) ) , abs (Err_thetaR5 (aa) ) ) ; 

end 
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9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-  Tv  (Z,  9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9- 

ooooooooooooooooooooo  lA  O  ooooooooooooooooooooooooooo 


for  CombinR=18500 : 200 : 50000 


for  thetaR=0 : 0 . 5 : 360 


RrAct=  (  (CombinR)  ''2  - 

L6''2)  /  (2*  (CombinR+L6*sin  (thetaR*pi/180)  )  )  ; 

RtAct=sqrt (RrAct^2+L6^2-2*RrAct*L6*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round ( (180/pi) *acos ( (RrAct^2+RtAct^2- 
L6''2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 

%  Organise  the  RCS  values  by 
RCS=10''  (biRCS3  (90-Betal,  90+Betal)  /lO)  ; 

K1  =  Pt*Gt*Gr* (lambda^2) *RCS* (Ft^2) * (Fr^2) ; 

K2  =  ( (4*pi) ^3) *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L6/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct''2*RtAct''2)  ; 

SNRdBm  =  10*log (SNR) +30; 

Xcoord  =  RrAct*sin ( (thetaR+Tx6) *pi/180) ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx6) *pi/180) ; 


%  Error  to  measurement  of  Rt  +  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

edge  of  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 


as 

%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

( l+e''2+2*e*sin  (thetaR*pi/180)  )  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ^2)  ; 

Err_time=Err_factor_time* (BeamEL/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  GOnmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e'^2  +  1)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ^2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e^2) *cos (thetaR*pi/180) / (2* (l+e*sin (thetaR*pi/180) ) ^2) ; 

Err_thetaR  =  Err  factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

XcoordG  =  [Xcoord6 , Xcoord] ; 

YcoordG  =  [ YcoordG , Ycoord] ; 

SNR6  =  [SNR6,SNR]; 

SNRdBm6  =  [ SNRdBm6 , SNRdBm] ; 

Err_time6  =  [Err_time6, Err_time] ; 

Err_L6  =  [Err_L6 , Err_L] ; 

Err_thetaR6  =  [Err_thetaR6 , Err_thetaR] ; 
else 
end 
end 

end 

Err_Total6=zeros (size (Err_time6) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err  time6) 

Err_Total6 (aa)  =  max (abs (Err_time6 (aa) ) , abs (Err_L6 (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err_time6 ) 

Err_Total6 (aa)  =  max (abs (Err_Total6 (aa) ) , abs (Err_thetaR6 (aa) ) ) ; 

end 
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9-9-9-9-9-9'9-9-9-9-9-9-9-9-9-9-9-9'9-9'9-  T'\r  1  9'9-9-9-9'9-9'9-9-9-9'9-9-9-9-9-9-9-9-9-9-9-9'9-9-9-9- 
ooooooooooooooooooooo  lA  /  ooooooooooooooooooooooooooo 


for  CombinR=18500 : 200 : 50000 


for  thetaR=0 : 0 . 5 : 360 


RrAct=  (  (CombinR)  ''2  - 

L7''2)  /  (2*  (CombinR+L7*sin  (thetaR*pi/180)  )  )  ; 

RtAct=sqrt (RrAct^2+L7^2-2*RrAct*L7*cos (pi/2+thetaR*pi/180) ) ; 

Beta=round ( (180/pi) *acos ( (RrAct^2+RtAct^2- 
L7'^2)  /  (2*RrAct*RtAct)  )  )  ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 

%  Organise  the  RCS  values  by 
RCS=10^ (biRCS3 (90-Betal, 90+Betal) /lO) ; 

K1  =  Pt*Gt*Gr* (lambda^2) *RCS* (Ft^2) * (Fr^2) ; 

K2  =  (  (4*pi)  ■^3)  *k*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L7/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct''2*RtAct''2)  ; 

SNRdBm  =  10*log (SNR) +30; 

Xcoord  =  RrAct*sin ( (thetaR+Tx7) *pi/180) ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx7) *pi/180) ; 


edge  of  the 
as 


%  Error  to  measurement  of  Rt  +  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 

%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e'^2+2*e*sin  (thetaR*pi/180)  )  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_time=Err_factor_time* (BeamEL/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  60nmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err  factor_L  =  (- 

(e^2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e''2)  *cos  (thetaR*pi/180)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_thetaR  =  Err  factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoord?  =  [Xcoordl , Xcoord] ; 

Ycoordl  =  [ Ycoordl , Ycoord] ; 

SNR7  =  [SNR7,SNR]; 

SNRdBm?  =  [ SNRdBm? , SNRdBm] ; 

Err_time7  =  [Err_time7 , Err_time] ; 

Err_L7  =  [Err_L7 , Err_L] ; 

Err_thetaR7  =  [Err_thetaR7 , Err_thetaR] ; 

else 

end 

end 

end 

Err_Total7=zeros (size (Err_time7) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err  time?) 

Err_Total? (aa)  =  max (abs (Err_time? (aa) ) , abs (Err_L? (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err_time? ) 

Err_Total? (aa)  =  max (abs (Err_Total? (aa) ) , abs (Err_thetaR? (aa) ) ) ; 

end 
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9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  Tv  ft  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 
000000000000000000000  lA.  O  000000000000000000000000000 


for  CombinR=18500 : 200 : 50000 


for  thetaR=0 : 0 . 5 : 360 


RrAct= ( (CombinR) ^2  - 

L8^2) / (2* (CombinR+L8*sin (thetaR*pi/180) ) ) ; 

RtAct=sqrt  (RrAct''2+L8'^2-2*RrAct*L8*cos  (pi/2+thetaR*pi/180)  )  ; 

Beta=round  (  (180/pi)  *acos  (  (RrAct'^2+RtAct'^2- 
L8^2) / (2*RrAct*RtAct) ) ) ; 

%  Assuming  the  low  RCS  ship  is  heading  south  towards  the  Rx 
%  position.  Then,  Beta/2  is  incident  angle  on  one  side  and 
%  scattered  angle  on  the  other  side. 

Betal=round (Beta/2 ) ; 

%  Organise  the  RCS  values  by 
RCS=10''  (biRCS3  (90-Betal,  90+Betal)  /lO)  ; 

K1  =  Pt*Gt*Gr*  (lambda'^2)  *RCS*  (Ft''2)  *  (Fr'^2)  ; 

K2  =  (  (4*pi)  ^'3)  *k;*Ts*Bn*Lt*Lr; 

K  =  K1/K2;  %  Bi-static  Radar  Constant 


e=L8/ (RrAct+RtAct) ;  %  eccentricity 

SNR  =  K/ (RrAct''2*RtAct''2)  ; 

SNRdBm  =  10*log (SNR) +30; 

Xcoord  =  RrAct*sin ( (thetaR+Tx8) *pi/180)  ; 
Ycoord  =  RrAct*cos ( (thetaR+Tx8) *pi/180) ; 


%  Error  to  measurement  of  Rt  +  Rr 

%  Antenna  scan  rate  of  45  rpm  or  270  deg/s 

%  With  horizontal  beamwidth  of  1.8  deg,  assuming  that  the 

edge  of  the 

%  beamwidth  creates  enough  SNR  to  be  received  by  the  receiver 


as 

%  reflections  from  the  tgt 
%  From  eqn  3.18 
%  rms  error 


Err_f actor_time  = 

(l+e^2+2*e*sin  (thetaR*pi/180)  )  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_time=Err_factor_time* (BeamEL/2*pi/180) *RtAct;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 
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%Err_time3  =  (0 . 55/Bn) *l/sqrt (Bn*SNRl) *c;  %  Ref  to  Alogrithms  for 
Ambiguity  Function  Processing  (Eqn  -  11b) 

%Err_timel  =  Err_time3  +  Err_timel; 

%  Error  to  measurement  of  L 

%  For  a  Class  A  Marine  AIS,  positional  acciracy  is  up  to 
0.0001  mins  where 

%  1  deg  equals  to  GOnmi 

%  Update  rate  is  between  2  to  10  seconds 

%  Assumming  a  opening  speed  of  Tx  at  12  knots  and  Rx  at  12 

knots 

Err_factor_L  =  (- 

(e''2  +  l)  *sin  (thetaR*pi/180)  +2*e)  /  (2*  (l+e*sin  (thetaR*pi/180)  )  ''2)  ; 

Err_L  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ) ;  %  Net 
opening  speed  of  10  knots  +  5s  update  rate 

%  Error  to  measurement  of  thetaR 
%  Error  due  to  thetaR 

%  the  error  will  be  coming  from  the  DF  error  from  the  ESM 

system 

%  which  is  up  to  +/-1.5  deg 

%  I  modified  this  equation  as  I  have  factor  in  the  L  in  my 

error  so  I 

%  normalised  it  with  L  so  as  not  to  double  count 
Err_f actor_theta  =  (1- 

e^2) *cos (thetaR*pi/180) / (2* (l+e*sin (thetaR*pi/180) ) ^2) ; 

Err_thetaR  =  Err_factor_theta* (RtAct*l . 5*pi/180) ; 

if  SNRdBm  >=  sen 

Xcoord8  =  [Xcoord8 , Xcoord] ; 

Ycoord8  =  [ Ycoord8 , Ycoord] ; 

SNR8  =  [SNR8,SNR]; 

SNRdBm8  =  [ SNRdBm8 , SNRdBm] ; 

Err_time8  =  [Err_time8 , Err_time] ; 

Err_L8  =  [Err_L8 , Err_L] ; 

Err_thetaR8  =  [Err_thetaR8 , Err_thetaR] ; 

else 

end 

end 

end 

Err_Total8=zeros (size (Err_time8) ) ; 

%  Find  th  emax  of  the  three  source  of  error  and  use  that  for  that 

%  coordinate  position 

for  aa=l : 1 : length (Err_time8 ) 

Err_Total8 (aa)  =  max (abs (Err_time8 (aa) ) , abs (Err_L8 (aa) ) ) ; 

end 

for  aa=l : 1 : length (Err_time8 ) 

Err_Total68 (aa)  =  max (abs (Err_Total8 (aa) ) , abs (Err_thetaR8 (aa) ) ) ; 

end 
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9-&-9-&-9-9-9-&-9'9-9-&-9-&-9-9-9'9-9-9-  &-9-9-9-9-9-9-9-&-9-9-9-9-9-9-9-&-9-9-9-9-9-&-9- 

00000000000000000000  irJ_CJLo  oooooooooooooooooooooooo 

X=horzcat (Xcoordl , Xcoord2 , Xcoord3 , Xcoord4 , XcoordS , XcoordG , Xcoord7 , XcoordS 

)  ; 

Y=horzcat ( Ycoordl , Ycoord2 , Ycoord3 , Ycoord4 , YcoordS , YcoordG , Ycoord7 , YcoordS 

)  ; 

Z=hor zcat ( SNRdBml , SNRdBm2 , SNRdBm3 , SNRdBm4 , SNRdBmS , SNRdBm6 , SNRdBm7 , SNRdBmS 

)  ; 

XX=horzcat (Err_thetaRl , Err_thetaR2 , Err_thetaR3 , Err_thetaR4 , Err_thetaR5 ,  Er 
r_thetaR6 , Err_thetaR7 , Err_thetaR8 )  ; 

YY=horzcat (Err_timel , Err_time2 , Err_time3 , Err_time4 , Err_time5 , Err  timeG , Er 
r_time7 , Err_time8 ) ; 

ZZ=horzcat (Err_Ll , Err_L2 , Err  L3 , Err_L4 , Err_L5 , Err_L6 , Err  L7 , Err_L8 ) ; 
AA=horzcat (Err_Totall, Err_Total2, Err_Total3, Err_Total4, Err_Total5, Err_Tot 
al6, Err_Total7, Err_Total8) ; 

figure ( 1 ) 
load  seamount 

scatter (X, Y, 4 , Z , 'fill', 'o') 
colorbar ( ' EastOutside ' ) 

%title (' Detection  Coverage  in  Polar  Plot  (4  Tx)  ') 

xlabel (' Distance  in  m' , ' fontsize ' , 14) 

ylabel (' Distance  in  m' ,' fontsize ', 14) 

hold  on 

plot  (0,0,  ' k* ' ) 

text (40,0, ' Rx ' , ' Fontsize ' , 15 ) 

plot  (-12730, 12730,  ' k* ' ) 

text (-12730, 12930, 'Txl', ' FontSize ' , 15) 

plot  (12730, 12730,  ' k* ' ) 

text (12730, 12930, ' Tx2 ' , ' FontSize ' , 15) 

plot (12730, -12730, ' k* ' ) 

text (12730, -12930, ' Tx3 ' , ' Fontsize ', 15) 

plot  (-12730, -12730,  ' k* ' ) 

text (-12730, -12930, ' Tx4 ' , ' FontSize ' , 15) 

plot  (0, -18000,  ' k* ' ) 

text (200, -18000, ' Tx5 ' , ' Fontsize ', 15) 
plot  (0,  18000,  ' k* ' ) 

text (200, 18000, 'Tx6', ' Fontsize ', 15) 
plot  (18000, 0,  ' k* ' ) 

text (18000, 200, ' Tx7 ' , ' Fontsize ', 15) 
plot  (-18000, 0,  ' k* ' ) 

text (-18000, 200, 'Tx8', ' Fontsize ', 15) 
hold  off 
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APPENDIX  G.  MATLAB  SOURCE  CODE  FOR  UNCERTAINTY 

ELLIPSE 


%  Chong  Sze  Sing  (March  2014) 
%  8  Tx  Uncertainty  Ellipse 

TargetX=4500 ; 

TargetY=4500 ; 

%  Initializattion 
%For  Tx  1 
Txl=45; 

Ll=18000; 

%For  Tx2 
Tx2=135; 

L2=18000; 

%For  Tx3 
Tx3=225; 

L3=18000; 

%For  Tx4 
Tx4=315; 

L4=18000; 

%For  Tx5 
Tx5=0; 

L5=18000; 

%For  Tx6 
Tx6=90; 

L6=18000; 

%For  Tx7 
Tx7=179; 

L7=18000; 

%For  Tx8 
Tx8=270; 

L8=18000; 


%  Variables 

%  Parameters  are  based  Furuno  Marine  Radar  FAR-21x7  series  in  both  X- 
band 

%  and  S-band 

%  Bi-static  Radar  Equation 
%  Constants 

k  =  1.38e-23;  %  Boltzmann 

c  =  3e8;  %  Speed  of  light 
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Select  =  1;  %chose  which  radar  to  use 

%S-band  Radar  (1) 

BeamAz_S  =1.8;  % 

BeamEL_S  =25;  % 

Scanrate_S  =  45;  % 

Pt_S  =  25e3; 

Gt_S  =  10^  (28/10)  ; 

Gr_S  =2;  % 

f_S  =3050e6;  % 

lambda_S  =  c/f_S; 

Ft_S  =1;  % 

Radar  -  Chp  6 . 5 . 5 . 6 
Fr_S  =1;  % 

Ts_S  =  290+290*  (6-1) ;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using 

formula  Te=(NF-l)290  and  Ts=Te+290K 

Bn  S  =  100e6;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 

MHz) 

Lt_S  =2;  %  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp  6. 5. 5. 6 

Lr_S  =  2; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 

PRF_S  =  2000; 

TOT_S  =  BeamEL_S/ (Scanrate_S*60) ; 

N  S  =  TOT  S*PRF  S; 


in  degrees 
in  degrees 

in  rpm  for  range  of  24  nmi  or  higher 
%  Tx  peak  power  output 
%  Tx  Antenna  Gain  given  as  28dB 
Rx  Antenna  Gain 
frequency 

%  wavelength 

Prop  Factor  for  Tx-Tgt  path.  See  Advances  in  Bistatic 
Prop  Factor  for  Rx-Tgt  path 


%X-band  Radar  (2) 
BeamAz_X  =  0.95;  % 

BeamEL_X  =25;  % 

Scanrate_X  =45;  % 

Pt_X  =  25e3; 

Gt_X  =  10^  (31/10)  ; 
Gr_X  =2;  % 

f_X  =  9410e6;  % 

lambda_X  =  c/f_X; 
Ft_X  =1;  % 

Radar  -  Chp  6 . 5 . 5 . 6 
Fr_X  =1;  % 

Ts_X  =  290+290*  (6-1) 
formula  Te=(NF-l)290 
Bn_X  =500e6;  % 

MHz) 


in  degrees 
in  degrees 

in  rpm  for  range  of  24  nmi  or  higher 
%  Tx  peak  power  output 
%  Tx  Antenna  Gain  is  given  as  31  dB 
Rx  Antenna  Gain 
frequency 

%  wavelength 

Prop  Factor  for  Tx-Tgt  path.  See  Advances  in  Bistatic 
Prop  Factor  for  Rx-Tgt  path 

;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using 

and  Ts=Te+290K 

Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 


Lt_X  =2;  %  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp  6. 5. 5. 6 

Lr_X  =  2; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 

PRF_X  =  1500; 

TOT_X  =  BeamEL_X/ (Scanrate_X*360) ; 

N  X  =  TOT  X*PRF  X; 
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if  Select  ==1 

BeamEL  =  BeamEL_S;  %  in  degrees 
BeamAz  =  BeamAz_S;  %  in  degrees 

Scanrate  =  Scanrate_S;  %  in  rpm  for  range  of  24  nmi  or  higher 

Pt  =  Pt_S;  %  Tx  peak  power  output 

Gt  =  Gt_S;  %  Tx  Antenna  Gain  is  given  as  31  dB 

Gr  =  Gr_S;  %  Rx  Antenna  Gain 

f  =  f_S;  %  frequency 

lambda  =  lambda  S;  %  wavelength 

Ft  =  Ft_S;  %  Prop  Factor  for  Tx-Tgt  path.  See  Advances  in 

Bistatic  Radar  -  Chp  6. 5. 5. 6 

Fr  =  Fr_S;  %  Prop  Factor  for  Rx-Tgt  path 

Ts  =  Ts_S;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using  formula 

Te=(NF-l)290  and  Ts=Te+290K 

Bn  =  Bn_S;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 

MHz) 

Lt  =  Lt_S;  %  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp 

6. 5. 5. 6 


Lr  =  Lr_S ; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 
PRF  =  PRF_S; 

TOT  =  TOT_S; 

N  =  N_S; 

else 


BeamEL  = 
BeamAz  = 
Scanrate 
Pt  =  Pt_X 
Gt  =  Gt_X 
Gr  =  Gr_X 
f  =  f  X; 


BeamEL_X;  %  in  degrees 
BeamAz  X;  %  in  degrees 
=  Scanrate_X;  %  in  rpm  for  range 
;  %  Tx  peak  power  output 

;  %  Tx  Antenna  Gain  is  given  as 

;  %  Rx  Antenna  Gain 

%  frequency 


of  24  nmi 
31  dB 


or  higher 


lambda  =  lambda_X;  %  wavelength 

Ft  =  Ft_X;  %  Prop  Factor  for  Tx-Tgt  path.  See  Advances  in 

Bistatic  Radar  -  Chp  6. 5. 5. 6 

Fr  =  Fr_X;  %  Prop  Factor  for  Rx-Tgt  path 

Ts  =  Ts_X;  %  System  Noise  Temp  of  Rx  of  NF=6dB  using  formula 

Te=(NF-l)290  and  Ts=Te+290K 

Bn  =  Bn  X;  %  Noise  Bandwidth  of  ESM  Rx  Wideband  Channalised  (500 


MHz) 


Lt  =  Lt_X; 
6. 5. 5. 6 


%  Ref  to  See  Advances  in  Bistatic  Radar  -  Chp 


Lr  =  Lr_X; 

%  this  is  calculation  of  pulse  integration  sqrt (N) 
PRF  =  PRF_X; 

TOT  =  TOT_X; 

N  =  N_X; 

end 
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Rr=sqrt  (TargetX^2+TargetY''2 )  ; 

theta=atan (TargetX/TargetY) ;  %  Angle  made  with  True  North  into 
clockwise  direction 

S- 9- S- 9- S- 9- S- 9- S- 9- S- 9- S-  Tvl  nr^nr-r^o  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 
ooooooooooooo  oULli-Ot::  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve (a^2-2*Rr*a-L1^2+2*Rr*Ll*sin (theta) ) ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 

e=Ll/ (Rr+Rt) ;  %  eccentricity 

beta= (180/pi) *acos ( (Rr^2+Rt^2-L1^2) / (2*Rr*Rt) ) ;  %  I  am  using  beta  to 
add  to  the  angle  of 

Err_f actor_time  =  (l+e'^2+2*e*sin  (theta)  )  /  (2*  (l+e*sin  (theta)  ^2)  )  ; 
Err_timel=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta  corresponding 
to  half  of  the  beamwidth  of  1.8  deg 


Err_factor_L  =  abs ( (- (e^2  +  l) *sin (theta) +2*e) / (2* (l+e*sin (theta) ) ^2) )  ; 
Err_Ll  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  t  5s  update  rate 


Err_f actor_theta  =  (l-e''2)  *cos  (theta)  /  (2*  (l+e*sin  (theta)  )  ^^2)  ; 
Err_thetaRl  =  Err  factor_theta* (Rt*l . 5*pi/180) ; 

%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_timel  >  Err_thetaRl 
Err_Maxl=Err_timel ; 

Err_Minl=Err_thetaRl ; 

%  this  handles  the  four  different  quadrant  and  if  err_time 
dominates 

if  theta>=0  &&  theta<=pi/4 
anglel= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
anglel=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
anglel= (-theta+pi/4) ; 

else 

anglel=pi/2+ (-theta+pi/4) ; 

end 
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else 

Err  Maxl=Err_thetaRl ; 

Err_Minl=Err_timel ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle 1=- theta; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle l=pi/ 2 -theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle 1=- theta; 

else 

angle l=pi/ 2 -theta; 

end 

end 

%  Decide  what  angle  to  put  for  each  Tx  contribuion 

S- 9- S- 9- S- 9- S- 9- S- 9- S- 9- S-  Tv9  nr^n-rr^o  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 

ooooooooooooo  oULlXOt:  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve  (a''2-2*Rr*a-L2''2+2*Rr*L2*sin  (theta)  )  ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 

e=L2/ (Rr+Rt) ;  %  eccentricity 

beta=  (180/pi)  *acos  (  (Rr''2+Rt^2-L2^2)  /  (2*Rr*Rt)  )  ;  %  I  am  using  beta  to 
add  to  the  angle  of 

Err_f actor_time  =  (l+e'^2+2*e*sin  (theta)  )  /  (2*  (l+e*sin  (theta)  ''2)  )  ; 
Err_time2=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta  corresponding 
to  half  of  the  beamwidth  of  1.8  deg 


Err_factor_L  =  abs ( (- (e^2  +  l) *sin (theta) +2*e) / (2* (l+e*sin (theta) ) ^2) )  ; 
Err_L2  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 


Err_f actor_theta  =  (l-e''2)  *cos  (theta)  /  (2*  (l+e*sin  (theta)  )  ^^2)  ; 
Err_thetaR2  =  Err_factor_theta* (Rt*l . 5*pi/180) ; 


113 


%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_time2  >  Err_thetaR2 
Err_Max2=Err_time2 ; 

Err_Min2=Err_thetaR2 ; 

%  this  handles  the  four  different  quadrant  and  if  err_time 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle2= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle2=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
angle2= (-theta+pi/4) ; 

else 

angle2=pi/2+ (-theta+pi/4) ; 

end 


else 

Err_Max2=Err_thetaR2 ; 

Err_Min2=Err_time2 ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle2=- theta; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle2=pi/ 2 -theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle2=- theta; 

else 

angle2=pi/ 2 -theta; 

end 

end 

%  this  function  will  always  draw  the  majoraxis  on  the  horizontal  plane 
and 

%  minoraxis  in  vertical  plane 
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9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 

ooooooooooooo  lAO  OULiXOt::  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve  (a''2-2*Rr*a-L3^2+2*Rr*L3*sin  (theta)  )  ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 

e=L3/ (Rr+Rt) ;  %  eccentricity 

beta=  (180/pi)  *acos  (  (Rr'^2+Rt^2-L3'^2)  /  (2*Rr*Rt)  )  ;  %  I  am  using  beta  to 
add  to  the  angle  of 

Err_f  actor_time  =  (l+e^2+2*e*sin  (theta)  )  /  (2*  (l+e*sin  (theta)  ''2))  ; 
Err_time3=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta  corresponding 
to  half  of  the  beamwidth  of  1.8  deg 


Err_factor_L  =  abs  (  (-  (e''2  +  l)  *sin  (theta)  +2*e)  /  (2*  (l+e*sin  (theta)  )  ^2)  )  ; 
Err_L3  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 


Err_f actor_theta  =  (l-e''2)  *cos  (theta)  /  (2*  (l+e*sin  (theta)  )  ^'2)  ; 
Err_thetaR3  =  Err_factor_theta* (Rt*l . 5*pi/180) ; 

%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_time3  >  Err_thetaR3 
Err_Max3=Err_time3 ; 

Err_Min3=Err_thetaR3 ; 

%  this  handles  the  four  different  quadrant  and  if  err_time 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle3= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle3=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
angle3= (-theta+pi/4) ; 

else 

angle3=pi/2+ (-theta+pi/4)  ; 

end 
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else 

Err  Max3=Err_thetaR3 ; 

Err  Min3=Err_time3 ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle 3=- theta; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle 3=pi/2- theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle3=- theta; 

else 

angle3=pi/2 -theta; 

end 

end 

9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  '^rMTrr’O  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 

ooooooooooooo  oULl-LL-ti  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve (a^2-2*Rr*a-L4^2+2*Rr*L4*sin (theta) )  ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 

e=L4/ (Rr+Rt) ;  %  eccentricity 

beta= (180/pi) *acos ( (Rr^2+Rt^2-L4^2) / (2*Rr*Rt) ) ;  %  I  am  using  beta  to  add 
to  the  angle  of 

Err_f actor_time  =  (l+e''2+2*e*sin  (theta)  )  /  (2*  (l+e*sin  (theta)  ''2)  )  ; 
Err_time4=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta  corresponding 
to  half  of  the  beamwidth  of  1.8  deg 


Err_factor_L  =  abs ( (- (e^2+l) *sin (theta) +2*e) / (2* (l+e*sin (theta) ) ^2) ) ; 
Err_L4  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 


Err_f actor_theta  =  (l-e'^2)  *cos  (theta)  /  (2*  (l+e*sin  (theta)  ) ''2)  ; 
Err_thetaR4  =  Err  factor_theta* (Rt*l . 5*pi/180) ; 
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%  Decide  which  one  is  the  major  and  minor  axis 
if  Err  time4  >  Err  thetaR4 

Err  Max4=Err  time4; 

Err  Min4=Err  thetaR4; 

%  this  handles  the  four  different  quadrant  and  if  err  time 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle4= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle4=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
angle4= (-theta+pi/4) ; 

else 

angle4=pi/2+ (-theta+pi/4) ; 

end 

else 

Err  Max4=Err  thetaR4; 

Err  Min4=Err  time4; 

%  this  handles  the  four  different 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle 4 =- theta; 

quadrant  and  if  err  theta 

elseif  theta>pi/4  &&  theta<=pi/2 
angle 4 =pi/2- theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle4=- theta; 

else 

angle4=pi/2 -theta; 

end 

end 

9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  TvR  ‘^rMTrr’O  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 
0000000000000  oULlXL-ti  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve (a^2-2*Rr*a-L5^2+2*Rr*L5*sin (theta)  )  ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 
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e=L5/ (Rr+Rt) ;  %  eccentricity 

beta= (180/pi) *acos ( (Rr^2+Rt^2-L5^2) / (2*Rr*Rt) ) ;  %  I  am  using  beta  to 
add  to  the  angle  of 

Err_f actor_time  =  (l+e^2+2*e*sin (theta) ) / (2* (l+e*sin (theta) ^2) ) ; 
Err_time5=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 

Err_factor_L  =  abs  (  (-  (e''2  +  l)  *sin  (theta)  t2*e)  /  (2*  (l+e*sin  (theta)  )  '^2)  )  ; 
Err_L5  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 

Err_f actor_theta  =  (l-e''2)  *cos  (theta)  /  (2*  (l+e*sin  (theta)  )  ^^2)  ; 
Err_thetaR5  =  Err  factor_theta* (Rt*l . 5*pi/180) ; 

%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_time5  >  Err_thetaR5 
Err_Max5=Err  time5; 

Err_Min5=Err_thetaR5 ; 

%  this  handles  the  four  different  quadrant  and  if  err_time 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle5= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle5=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
angle5= (-theta+pi/4) ; 

else 

angle5=pi/2+ (-theta+pi/4) ; 

end 

else 

Err_Max5=Err_thetaR5 ; 

Err_Min5=Err_time5 ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle5=- theta; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle 5=pi/2- theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle5=- theta; 

else 

angle5=pi/2 -theta; 

end 

end 


118 


9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  Tvin  ‘^rMTrr’O  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 
0000000000000  XAu  oUL-lJ-L-t:  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve (a^2-2*Rr*a-L6^2+2*Rr*L6*sin (theta) ) ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 

e=L6/ (Rr+Rt) ;  %  eccentricity 

beta= (180/pi) *acos ( (Rr^2+Rt^2-L6^2) / (2*Rr*Rt) ) ;  %  I  am  using  beta  to 
add  to  the  angle  of 

Err_f actor_time  =  (l+e'^2+2*e*sin  (theta)  )  /  (2*  (l+e*sin  (theta)  ^2)  )  ; 
Err_time6=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 


Err_factor_L  =  abs ( (- (e^2  +  l) *sin (theta) +2*e) / (2* (l+e*sin (theta) ) ^2) )  ; 
Err_L6  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 


Err_f actor_theta  =  (l-e''2)  *cos  (theta)  /  (2*  (l+e*sin  (theta)  )  ^2)  ; 
Err_thetaR6  =  Err_factor_theta* (Rt*l . 5*pi/180) ; 

%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_time6  >  Err_thetaR6 
Err  Max6=Err_time6 ; 

Err_Min6=Err_thetaR6 ; 

%  this  handles  the  four  different  quadrant  and  if  err_time 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle6= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle6=pi/2+ (-theta-pi/4) ; 

elseif  thetakO  &&  theta>=-pi/4 
angle6= (-theta+pi/4) ; 

else 

angle6=pi/2+ (-theta+pi/4)  ; 

end 
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else 

Err  Max6=Err_thetaR6 ; 

Err  Min6=Err_time6 ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle 6=- theta; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle6=pi/2- theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle 6=- theta; 

else 

angle6=pi/2 -theta; 

end 

end 

9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9-  Tv7  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 

0000000000000  ±  2i.  /  oULiJ-L-t:  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve  (a''2-2*Rr*a-L7^2+2*Rr*L7*sin  (theta)  )  ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 

e=L7/ (Rr+Rt) ;  %  eccentricity 

beta=  (180/pi)  *acos  (  (Rr'^2+Rt''2-L7'^2)  /  (2*Rr*Rt)  )  ;  %  I  am  using  beta  to 
add  to  the  angle  of 

Err_f actor_time  =  (l+e^2+2*e*sin (theta) ) / (2* (l+e*sin (theta) ^2) ) ; 
Err_time7=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta 
corresponding  to  half  of  the  beamwidth  of  1.8  deg 


Err_factor_L  =  abs  (  (-  (e''2  +  l)  *sin  (theta)  +2*e)  /  (2*  (l+e*sin  (theta)  )  ^2)  )  ; 
Err_L7  =  Err_factor  L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 


Err_f actor_theta  =  (l-e^2) *cos (theta) / (2* (l+e*sin (theta) ) ^2) ; 
Err_thetaR7  =  Err_factor_theta* (Rt*l . 5*pi/180) ; 
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%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_time7  >  Err_thetaR7 
Err_Max7=Err_time7 ; 

Err_Min7=Err_thetaR7 ; 

%  this  handles  the  four  different  quadrant  and  if  err_time 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle7= (-theta-pi/4) ; 

elseif  theta>pi/4  &&  theta<=pi/2 
angle7=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
angle7= (-theta+pi/4) ; 

else 

angle7=pi/2+ (-theta+pi/4) ; 

end 


else 

Err_Max7=Err  thetaR7; 

Err_Min7=Err_time7 ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle 7 =- theta, • 

elseif  theta>pi/4  &&  theta<=pi/2 
angle7=pi/ 2 -theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle 7 =- theta; 

else 

angle7=pi/ 2 -theta; 

end 

end 

9- S- 9- S- 9- S- 9- S- 9- S- 9- S- 9-  Tvfi  c^r^iiT-r’o  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 
ooooooooooooo  ±A.O  oUQ-LL-ti  ooooooooooooooo 

%  Eqn  5 . 1 
syms  a 

syms  A  positive 

A=solve  (a'^2-2*Rr*a-L8'^2+2*Rr*L8*sin  (theta)  )  ; 

RrRt=double (A) ; 

Rt=abs (RrRt (1) ) -Rr; 
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e=L8/ (Rr+Rt) ;  %  eccentricity 

beta=  (180/pi)  *acos  (  (Rr’^2+Rt''2-L8'^2)  /  (2*Rr*Rt)  )  ;  %  I  am  using  beta  to  add 
to  the  angle  of 

Err_f actor_time  =  (l+e^2+2*e*sin (theta) ) / (2* (l+e*sin (theta) ^2) ) ; 
Err_time8=Err_factor_time* (BeamAz/2*pi/180) *Rt;  %  r*theta  corresponding 
to  half  of  the  beamwidth  of  1.8  deg 

Err_factor_L  =  abs ( (- (e^2+l) *sin (theta) +2*e) / (2* (l+e*sin (theta) ) ^2) ) ; 
Err_L8  =  Err_factor_L* ( (0 . 0001*1824) + (10*1824*5/3600) ); %  Net  opening 
speed  of  10  knots  +  5s  update  rate 

Err_f actor_theta  =  (l-e^2) *cos  (theta) /  (2*  (l+e*sin (theta) ) ^2) ; 

Err_thetaR8  =  Err  factor_theta* (Rt*l . 5*pi/180) ; 

%  Decide  which  one  is  the  major  and  minor  axis 
if  Err_time8  >  Err_thetaR8 
Err_Max8=Err_time8 ; 

Err_Min8=Err_thetaR8 ; 

%  this  handles  the  four  different  quadrant  and  if  err_time  dominates 
if  theta>=0  &&  theta<=pi/4 
angle8= (-theta-pi/4) ; 
elseif  theta>pi/4  &&  theta<=pi/2 
angle8=pi/2+ (-theta-pi/4) ; 

elseif  theta<0  &&  theta>=-pi/4 
angle8= (-theta+pi/4) ; 

else 

angle8=pi/2+ (-theta+pi/4) ; 

end 

else 

Err  Max8=Err_thetaR8 ; 

Err_Min8=Err_time8 ; 

%  this  handles  the  four  different  quadrant  and  if  err_theta 
dominates 

if  theta>=0  &&  theta<=pi/4 
angle 8 =- theta, • 

elseif  theta>pi/4  &&  theta<=pi/2 
angle8=pi/2- theta; 

elseif  theta<0  &&  theta>=-pi/4 
angle 8 =- theta; 

else 

angle8=pi/2 -theta; 

end 

end 
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9- 9- 9- 9- 9- 9- 9- 9-  P  1  -I- -I-  -i  n  rr  9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 

00000000  lrJ_ULL-Liiy  ooooooooooooooooooo 

pdeellip (TargetX, TargetY, Err  Maxi, Err  Mini, anglel- (Txl*pi/180) , 
pdeellip (TargetX, TargetY, Err  Max2,Err  Min2, angle2- (Tx2*pi/180) , 
pdeellip (TargetX, TargetY, Err  Max3,Err  Min3, angle3- (Tx3*pi/180) , 
pdeellip (TargetX, TargetY, Err_Max4, Err_Min4, angle4- (Tx4*pi/180) , 
pdeellip (TargetX, TargetY, Err  Max5,Err  Min5, angleS- (Tx5*pi/180) , 
pdeellip (TargetX, TargetY, Err  Max6,Err  Min6, angleG- (Tx6*pi/180) , 
pdeellip (TargetX, TargetY, Err_Max7, Err_Min7, angle7- (Tx7*pi/180)  , 
pdeellip (TargetX, TargetY, Err_Max8, Err_Min8, angle8- (Tx8*pi/180)  , 

hold  on 
axis  square 
grid  on 
plot  (0,0,  ' k* ' ) 

text(40,0, ' Rx ' , ' FontSize ' , 14 ) 

plot (-12730, 12730, ' k* ' ) 

text (-12730, 12930, ' Txl ' , ' FontSize ', 15) 

plot (12730, 12730, ' k* ' ) 

text (12730, 12930, ' Tx2 ' , ' FontSize ', 15) 

plot (12730, -12730, ' k* ' ) 

text (12730, -12930, ' Tx3 ' , ' FontSize ', 15) 

plot (-12730, -12730, ' k* ' ) 

text (-12730, -12930, ' Tx4 ' , ' FontSize ', 15) 

plot  (0, -18000,  ' k* ' ) 

text (200, -18000, ' Tx5 ' , ' FontSize ', 15) 
plot  (0,  18000,  ' k* ' ) 

text (200, 18000, 'Tx6', ' FontSize ', 15) 
plot  (18000, 0,  ' k* ' ) 

text (18000, 200, ' Tx7 ' , ' FontSize ', 15) 
plot  (-18000, 0,  ' k* ' ) 

text (-18000, 200, 'Tx8', ' FontSize ', 15) 
plot (4500, 4500, ' k* ' ) 

text (4600, 4500, 'Target', ' FontSize ', 15) 
hold  off 


'Txl ' ) ; 
' Tx2 ' ) ; 
' Tx3 ' ) ; 
' Tx4 ' ) ; 
' Tx5 ' ) ; 
'Tx6' ) ; 
' Tx7 ' ) ; 
' Tx8 ' ) ; 
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